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ABSTRACT OF THE DISSERTATION 
QCD STRUCTURE OF NUCLEAR INTERACTIONS 

by 

Carlos G. Granados 
Florida International University, 201 1 
Miami, Florida 
Professor Misak Sargsian, Major Professor 

The research presented in this dissertation investigated selected processes involv- 
ing baryons and nuclei in hard scattering reactions. These processes are characterized 
by the production of particles with large energies and transverse momenta. Through 
these processes, this work explored both, the constituent (quark) structure of baryons 
(specifically nucleons and A-Isobars), and the mechanisms through which the interac- 
tions between these constituents ultimately control the selected reactions. 

The first of such reactions is the hard nucleon-nucleon elastic scattering, which was 
studied here considering the quark exchange between the nucleons to be the dominant 
mechanism of interaction in the constituent picture. In particular, it was found that 
an angular asymmetry exhibited by proton-neutron elastic scattering data is explained 
within this framework if a quark-diquark picture dominates the nucleon's structure in- 
stead of a more traditional SU{6) three quarks picture. The latter yields an asymmetry 
around 90° center of mass scattering with a sign opposite to what is experimentally 
observed. 

The second process is the hard breakup by a photon of a nucleon-nucleon system in 
light nuclei. Proton-proton (pp) and proton-neutron (p??) breakup in ^He, and AA-isobars 
production in deuteron breakup were analyzed in the hard rescattering model (HRM), 
which in conjunction with the quark interchange mechanism provides a Quantum Chro- 
modynamics (QCD) description of the reaction. Through the HRM, cross sections for 
both channels in ^'He photodisintegration were computed without the need of a fitting 
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parameter. The results presented here for pp breakup show excellent agreement with 
recent experimental data. 

In AA-isobars production in deuteron breakup, HRM angular distributions for the 
two AA channels were compared to the pn channel and to each other. An important 
prediction from this study is that the A++ A- channel consistently dominates A+A", which 
is in contrast with models that unlike the HRM consider a AA system in the initial state 
of the interaction. For such models both channels should have the same strength. These 
results are important in developing a QCD description of the atomic nucleus. 
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Chapter I 
INTRODUCTION 



This introduction presents an overview of the progress made in nuclear physics towards 
achieving an unified description of the strong nuclear force; a description that can match 
the phenomenological successes of the standard model of electroweak interactions. In 
line with the specific goals of the studies presented in this dissertation, the following 
sections place emphasis on the search for an elementary particle description of nucleon 
nucleon, and more generally hadron hadron interactions. Section I.l addresses the dis- 
covery of the nucleus and of its constituents (nucleons) the proton and the neutron as 
well as the concerns that led to the discovery of the nuclear force. Section 1.2 describes 
the road towards formulating a field theory of the strong nuclear force motivated by 
the quantum electrodynamics (QED) description of electromagnetic interactions. It 
also outlines the experimental and phenomenological work that lead to the discovery 
of quarks and gluons and to the formulation of quantum chromodynamics (QCD) as 
the fundamental field theory of the strong interaction. Section 1.3 focuses on general 
features of hard exclusive processes, and on how the results from QCD phenomenology 
explored in section 1.2 can be used here to investigate the QCD mechanisms through 
which the particles involved in these processes interact. This section introduces the 
framework and some aspects of the methodology in which the processes of interest in 
this dissertation are studied in the chapters that follow. 

I.l The Strong Interaction 

The discovery of subatomic particles in the late XIX and early XX century in conjunc- 
tion with the formulation of special relativity and Quantum mechanics steered major 
efforts and developments in physics towards identifying the fundamental structure of 
matter. These developments resulted in the formulation of the Standard Model (SM) 
of particle physics, which has met outstanding phenomenological success in describ- 
ing electromagnetic and weak interactions. All the fundamental elementary particles 
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(fermions) are subject to either of these two interactions; more generally, they are sub- 
ject to the electroweak interaction which unifies the two interactions under a consis- 
tent mathematical framework. The validation of the SM has been achieved through 
the experimental programs that test phenomenological predictions emanating from this 
framework. Such tests are facilitated by major advances in accelerator technologies 
that allow studies at ever increasing energy regimes and much improved detection and 
particle identification capabilities. These conditions have made possible the production 
and direct detection of leptons and of the otherwise hypothetical weak bosons which 
mediate the weak interaction and that are predicted by the standard model. 

On the other hand, the standard model picture of the strong interaction faces a more 
challenging phenomenological treatment. In the SM, the strong interaction is described 
by quantum chromodynamics (QCD), a quantum field theory in which quarks and glu- 
ons are respectively the elementary fermions and bosons. But these quarks and gluons 
are confined in hadrons (a family of particles of which protons and neutrons are part of) 
and have not been directly observed as free particles. However, as is discussed in the fol- 
lowing sections, the existence of the strong force was postulated to explain the stability 
of the nucleus as a system of protons and neutrons (nucleons); thus in a first approach, 
a theory of the strong interaction would have nucleons as its fundamental elementary 
particles to later on include mesons as the elementary bosons. The origins of this de- 
scription are summarized below and in section 1.2, in which it's also emphasized that its 
validity is limited to separation distances at which the interacting nucleons can still be 
considered elementary particles. At separation distances smaller than Ifm (lO^^m), the 
nucleons' internal structures (quark-gluon distributions) play a more active role in the 
interaction, and approaches assuming elementary (structureless) nucleons are no longer 
suited for describing phenomena at this scale. 

After introducing QCD, the theory of interactions of quarks and gluons, in section 
1.2, section 1.3 comments on general features that allow the use of QCD approaches 
for developing a quantitative description of the strong nuclear force at separation dis- 
tances in this regime (<lfm). In the chapters that follow, the reactions of interest in this 
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dissertation are investigated within the framework of QCD approaches. 
The Atomic Nucleus 

The atomic nucleus concentrates most of the atomic mass in a radius about 10® times 
smaller than the atomic radius. It also carries the net positive charge that balances the 
electronic cloud surrounding it in a neutral atom. This accepted picture of the atom 
was partially completed by Rutherford in 1911 [1] when results from a scattering ex- 
periment of alpha particles incident upon a thin film of gold showed that some small, 
but larger than expected, portion of the beam was deflected at large angles while most 
of the beam passed through the film with almost no deflection [2]. Rutherford showed 
that under a coulomb interaction, the observed deflection of alpha particles at large an- 
gles could only be explained by a large charge and mass concentration in a very small 
volume at the center of each gold atom. The picture of the atom then switched from 
the one in which the positive charge was uniformly distributed throughout the atom 
(J. Thompson's model) and according to which the large angle deflection of the alpha 
beam would had been much more suppressed, to the one in which the mass and the net 
positive charge were concentrated in a small volume (nucleus) surrounded by orbiting 
electrons that expand to the atom's size [1]. 

This 'planetary' model of the atom was still plagued with flaws. The most impor- 
tant being the fact that the atom's stability under this model would contradict classical 
electrodynamics. The orbiting electrons are accelerating and therefore should radiate 
electromagnetic waves losing kinetic energy to their electromagnetic field and therefore 
decaying by spiraling down towards the nucleus under the Coulomb force. This issue 
was then solved by the emerging quantum mechanics picture of the atom, in which the 
orbital structure of the electronic distribution was replaced by discrete quantum states 
of the electronic field surrounding the nucleus. An atom could only radiate or absorb 
energy when transitioning between these states, while the lowest energy state has still a 
definite finite energy. 
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The Strong Nuclear Force 

Concerns about the composition of the nucleus were still to be addressed as well. One 
of the early models of the nucleus in which it consisted of protons and electrons has 
stability issues. The huge repulsing force that positive charges would experience at the 
nuclear distances could not be balanced by the overall charge distribution of the atom. 
The nucleus would quickly disintegrate under the Coulomb repulsion. Therefore, a 
force of a different nature was thought to be responsible for keeping the nucleus' com- 
ponents in place. This nuclear force should be much stronger than the electromagnetic 
force at nuclear distances, but it should have a range limited to the nuclear distances as 
well since the evidence of this force is not found anywhere else outside the nucleus. 

Nuclear Composition 

A lot of progress was made in understanding the nuclear structure of atoms by assum- 
ing the existence of this force, but without a detailed theory of its origin the successes 
were only confined to some general features. This picture was much improved by the 
discovery of the neutron by James Chadwick in 1932 [3]. The nucleus was now thought 
of as a system of nucleons (protons (positive charge carriers) and neutrons (of neutral 
charge)) confined by an effective nuclear potential. Both neutrons and protons are the 
sources of the force field generated by this potential which at the nuclear scale is at- 
tractive. Because neutrons are not meaningful sources of an electromagnetic field, their 
presence in nuclei helped explain the stability trend of the atoms with increasing atomic 
number. As the atomic number (the number of protons in the nucleus) increases, the 
ratio of the number of neutrons to protons has to increase as well in order to balance the 
increasing electromagnetic repulsion felt by the protons with an increasing number of 
attractive nuclear force sources. The short range of the nuclear force will also constrain 
the size of the stable nuclei. 

The existence of neutrons in nuclei also solved the contradiction regarding the ob- 
served spin of nuclei inherent to the early (proton-electron) nuclear models. Under 
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these models the nucleus was made of 14 protons and 7 electrons, both particles 
with total spin 1/2. Their corresponding spins will add up to a system with a total half 
integer spin. However, it was observed that ^^N obeys Bose statistics which implies that 
its total spin was of integer value. Because, like protons and electrons, neutrons are 
fermions (particles with half integer spin ), the nucleus made up of 7 protons and 7 
neutrons would naturally account for expected total integer spin of the "N nucleus. 

Just as for the atom, a quantum mechanics picture of the nucleus emerged in which 
discrete nuclear states correspond to the different configurations of protons and neutrons 
filling different nuclear energy levels. Transitions between nuclear states could involve 
energy absorption or energy emission in the form of radiation. A nucleus with fixed 
atomic and mass numbers would have a unique photon emission spectrum coming from 
the decays of its allowed exited states just as an atom of a given element would. A decay 
to a different nuclear state in which the atomic and/or the mass numbers are altered is 
signaled by the emission of massive charged or neutral particles such as neutrons, alpha 
particles and beta particles. At this stage, it was understood that the two (and more) 
body nucleon-nucleon interactions had to be the source of the nuclear potentials which 
makes a nucleus a bound system. 

Force Carriers 

A picture of nucleon-nucleon (NN) interaction was however yet to be reached to the 
level at which for instance the electromagnetic interaction was then understood. Through 
quantum electrodynamics (QED), the relativistic quantum field theory of electromag- 
netism, charged particles interact by transferring energy and momentum to one another 
through the exchange of field quanta (photons). In 1935, H. Yukawa proposed that just 
as in electromagnetism the Coulomb force was the result of the exchange of virtual 
photons, the strong nuclear force as well results from the exchange of corresponding 
virtual bosons (mesons) [4]. The short range nature of this strong force comes from 
these mesons having nonzero masses. While in the non-relativistic limit the photon 
exchange interaction among charged particles can be effectively described by the clas- 
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sic form of the Coulomb potential (--1/r), the strong interaction mediated by massive 
mesons is in the nonrelativistic limit dominated by the Yukawa potential (e "' /r), with 
n being the mass of the exchanged meson. Therefore, the mass of the exchange meson 
would define the range of the force. A candidate for the particle mediating the nucleon 
nucleon interaction was found in 1947. Known as the n meson or pion, it has zero 
spin (scalar) and could have a charge of +e (n+), -e{n-) or no charge (tt") with a mass 
^ --140MeV. Later on, other heavier mesons of spin 1 (vector mesons) p and w were 
discovered [5] [6]. Contribution from the exchange of these vector mesons to the (NN) 
interaction are thought to be relevant in understanding the behavior of the (NN) force 
at distances far smaller than those were this force is acceptably described by either the 
pion exchange model or the Yukawa potential However, as the distance between the 
interacting nucleons keeps getting smaller, the particle exchange mechanisms required 
to explain the empirical behavior of the strong force grow in number and complexity. 

1.2 From Hadrons to QCD 

The meson field theory for the strong interaction was motivated by the phenomeno- 
logical success of QED. This theory tries to include all possible interactions between 
nucleons in such a way that preserves the increasing number of global symmetries. Evi- 
dence of these symmetries- arose from the discovery not only of new mesons, but also of 
new fermions that seem to interact via the strong interaction. The latter was concluded 
from the observed high rates ( half lives of ^ iO^*s) at which these fermions were decay- 
ing into nucleons that could only be explained by strong interaction couplings. These 
fermions along with the lighter and more stable nucleons came to form a new category 
of particles named baryons. Mesons and baryons form the group of particles that can 
interact through the strong interaction also known as hadrons. While baryons played a 
role analogous to that of fermions in QED, mesons as mentioned earlier assumed a role 
analogous to that of photons, as the carriers of the nuclear force. 

'See e.g. Ref.[7] for an extended overview of the meson theory of nuclear interactions. 
-Each symmetry of the strong interaction was inferred from experimental evidence of certain conservation 
laws, some of which included the conservation of isospin and strangeness. 

6 



Isospin 



The nuclear force for instance contributes with the same strength in both proton proton 
(jrp) and neutron neutron (nn) interaction. This experimental fact suggests that under 
the strong interaction protons and neutrons are different states of the same particle. In 
analogy with the spin formalism, this set of two states is labeled by a quantum number 
called isospin (/). Each state is in this set is labeled by another quantum number (t) 
analogous to the spin projection (m). Each state then can be transformed into the other 
by a rotation in this internal 'isospace'. By convention, for a proton (neutron) / = 1/2 
and T = +i/2( 1/2). Then a pp system is turned into a nn system through a rotation in 
the isospace. As mentioned earlier the strength of the force remains the same after this 
rotation is performed. Therefore, the strong interaction is symmetric under rotations 
in the isospace. For now, the known mesons, tt, p and w have integer isospins, and the 
isospin symmetry of the strong interaction holds on reactions involving them as well. 
Isospin conservation puts specific restrictions on the structure of the nn interactions 
through meson exchange. 

Meson field theoretic approach 

The isospin invariance of the strong interaction is utilized to obtain an appropriate rep- 
resentation of the elements of a meson exchange theory of the strong interaction. This 
invariance is also taken into account when writing the equations that govern the general 
dynamics under this interaction. In this framework, protons and neutrons are repre- 
sented by Pauli spinors from 5*1/(2)/ fundamental representation, while mesons (the field 
quanta) belong to the triplet representation of su{2) (see e.g. Refs. [7] and [8]): 




(I-l) 
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and 



(1-2) 



V 



which transform under an infinitesimal SU{2)j rotation (parametrized by an isovector e) 
following, 



*jv ^ (1 - • t/2)*jv 



(1-3) 
(1-4) 



T is a vector formed by Pauli matrices (the generators of transformations in SU{2)). 
For a nucleon ^-jv in a pion field 



(ly - mf,)<^f,{x) = .goi75(T • $^(.t))\E'jv(x), (1.5) 

which corresponds to the Dirac equation modified by a canonical transformation ("min- 
imal substitution") of the momentum operator, ly ^ -goiir^{T-^^{x)), in analogy to what 
is done for a charged particle in a electromagnetic field to account for interactions. 75 
ensures parity conservation, and go is a coupling constant to be determined experimen- 
tally. Similarly, as in electromagnetism in the presence of a charged current J""' the 
electromagnetic field A{x) obeys VM(x) = J'='"(x), the pion field in the strong interaction 
obeys, 

- ml)<^^{x) = -,go*iv(a;)i75T*jv(a;), (1.6) 

in which *jv(a;)i75T*jv(x) is identified as the isotopic nucleon current. Eqs. (1.5) and (1.6) 
are covariant under SU{2)i. The homogeneous versions of these equations correspond to 
the Dirac and Klein-Gordon equations for a free nucleon and pion field respectively. In 
a more formal treatment Eqs. (1.5) and (1.6) can be obtained after minimizing the action 
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under the Lagrangian density 



i.e., from 

SS^6 J d'*x£,(*jv,5^«'^,*iv,$.,5^$.) = 0, 
which yields the Euler-Lagrange Equations 



(1-7) 



resuhing in Eq.1.5, and 



resuhing in Eq.1.6. 

The Lagrangian density in the form of Eq.(1.7) is invariant under the transformations 
in Eq.(1.4). This symmetry in conjunction with Euler-Lagrange equations leads to 

in which the quantity in parentheses is identified as the conserved current J" corre- 
sponding to the SU{2)i symmetry. Explicitly, from Eq.(L7), 

= 1/2^-^7^ r^-^ + X 9^$,) . (L8) 

From J" we then obtain the three components of isotopic charge of the system, 

1 = J d'xf = j d'x [l/2*]vT*iv + X ao$.)] , (L9) 

which is a constant of motion. Because the Lagrangian density in Eq.(L7) should satisfy 
the conservation of the electric charge, it should be also symmetric under c/(i) transfor- 
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mations on electrically charged fields, 




(1. 10) 



(I.ll) 



It can be shown that from this symmetry it follows that the electrical charge of the 
system Q is a constant of motion as well. Likewise, the baryonic number B, 

B = 2(Q - I,) 

is a constant of motion in interactions under in Eq.(1.7). 

Scattering matrix The fact that the conservation of b, i, and Q charges in strong in- 
teractions is well established experimentally drew efforts into developing a field theory 
of pions and nucleons based on a Lagrangian such as Eq.(1.7) and on the observed sym- 
metries. The free parameters of such a theory according to Eq.(1.7) would in principle 
be the nucleon and the pion masses, and the pion nucleon coupling constant g which 
are experimental observables. A quantitative analysis is then developed to evaluate the 
phenomenological accuracy of This is usually done by studying scattering processes 
for which the cross section can be calculated from entries of a scattering matrix 8^. 

The scattering matrix is the probability amplitude of a system (generally a mul- 
tiparticle system)in an initial state in which there's no interaction to evolve through 
interaction into a final state also away from the region of interaction, i.e.. 



in which cf), and 0/ are solutions of equations such as the homogeneous versions of Eqs. 
(1.5) and (1.6). In the Heisenberg picture (ses e.g. Ref [9]) Eq.(1.12) can be written 
making explicit the role of a time evolution operator: 



(1-12) 



(1.13) 
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in which the Hamiltonian H is given as a function of the field operators 0(x) and tt{x) = 
dQ(j){x) and T[...] stands for the time ordering of the product of operators in [...], e.g., 



These products come from the expansion of e Because the initial and final states 
in an scattering reaction are generally different, the nontrivial entries of the S matrix 
become the entries of a transition matrix 7 which relate to S through 



An invariant matrix element M is then defined by removing the condition of 4-momentum 
conservation, i.e.. 



If the interaction is weak, the invariant matrix elements can be computed from the 
lowest order terms of an expansion of Eq.(1.13). The terms of such expansion can be 
graphically represented by Feynman Diagrams such as the one shown in Fig.I.l. The 
vertex factors are determined from the interaction terms of the Lagrangian or Hamilto- 
nian, while the propagators or internal lines are obtained from the equations of motion 
derived from the free Lagrangian. e.g., from 

{q^ - m^)G'Jix, x') = S{x - xy 

we obtain the pion propagator. 



T[(j){x)(j){y)] = 0{xo - yo)(j){x)(j){y) + e{yo - Xo)<p{y)(f>{x) 



S = l + i7. 




(1-14) 



in which 



(1-15) 



11 



enters as a factor in calculating M in momentum space, as illustrated in Fig.I.l. Ex- 
ternal lines emerge from taking the T ^ oo limit in defining the scattering matrix as in 
Eq.(1.13). The corresponding factors for external lines are corresponding representa- 
tions in momentum space of solutions for the free equations of motion. 

N igWJp+q)y^r'WJp) 

I 
I 

7t 1 id ..(q^ -m 

I 

N' I 

L> 

Figure I.l: Nucleon Nucleon scattering through the one pion exchange 
mechanism. 

Feynman Diagrams The matrix elements M for nn scattering through the pion 
exchange mechanism are then computed based on an expansion in which the lowest 
order term in g (of second order) is illustrated in Fig. (I.l). This term corresponds to 
the one pion exchange (OPE) mechanism which is assumed in OPE models to suffi- 
ciently well describe the interaction at large distances (r>2fm). From Fig. (I.l), M for 
NN scattering within OPE can be expressed as, 

-zMoj.e =igo^iv(p + g)75r'^iv(p) ^ igo^^Mip' -ihsT'-^'j^ip'), (1.16) 

in which J^_'^^ is the pion's propagator according to Eq. (1.14). 

It is possible to obtain the interaction potential from Feynman diagrams in the non- 
relativistic limit. For instance, a radial potential can be derived from an angular inte- 
gration of the Fourier transform of the s-wave component of this scattering amplitude. 
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It yields (see e.g. [8]), 



Vs{r) = -2f 



r 




(1.17) 



in which = f| {j^f- The first term in Eq.(1.17) corresponds to the Yukawa poten- 
tial which properly describes the long-range attraction. A short-range repulsion is ac- 
counted for by the Dirac delta term in Eq.(1.17). However, because of its nonrelativistic 
character, Eq.(1.17) is only accurate for distances at which the Yukawa potential is ac- 
curate. In such a region, the parameter / and consequently can be set. It is found 
that experimentally 14, i.e., the ttN coupling constant is much larger than 1. The 
use of a traditional expansion of M in orders of go is not well justified in this case; the 
diagrams of higher order in higher are not necessarily suppressed. Such diagrams 
involve multiple pion exchanges and become more significant at smaller distances with 
a contribution to the nn potential behaving as e""'""'^ in which n corresponds to the 
number of pions exchanged. 

Calculations for diagrams with two or more pion exchanges are considerably more 
tedious and complicated and did not correlate well with experimental data. Alterna- 
tively, it was argued that multiple uncorrelated pion exchanges contribute little to the the 
NN force, and that instead additional heavier mesons in place of correlations between 
the exchanging pions would explain the nn interaction at smaller distances. Particu- 
larly, the short range repulsion and a spin orbit force arise naturally by considering one 
vector-meson exchange diagrams. The p and u vector mesons were discovered in the 
early 60s as resonances of li^ and Stt states respectively. Both mesons have spin 1 and 
isospins 1 and respectively. Accordingly, their couplings to nucleon currents are of 
the form 



which are added to the Lagrangian and/or Hamiltonian as interaction terms along with 
the corresponding kinetic energy terms. 

The strong attraction at the mid range of the nn force however requires the inclusion 




13 



of yet another kind of meson. It was named a- meson and should have spin and 
isospin, and a mass in the range of 400-800 MeV. The existence of such a meson has 
not yet been established at least as a field particle in meson field theories. However, 
many other mesons were discovered increasing the density of degrees of freedom of 
the theory and consequently the number of interaction terms in the Lagrangian. Baryon 
resonances of nucleons and mesons such as A-isobars needed to be included as well 
which contributed largely to the previously neglected uncorrelated multipion exchange 
diagrams. The number of free parameters such as coupling constants grew dramatically 
and soon the Isospin invariance together with the more fundamental symmetries (parity 
and Lorentz invariance) did not seem to sufficiently constrain the Lagrangian of the 
theory to a limited form desired for a field theory of fundamental degrees of freedom 
such as in QED (that has only one coupling constant). 

Renormalization and the Landau pole The large number of free parameters meant 
that the developing theory lacked predictive power. Specially in the short range in 
which contributions other than the one meson exchange become relevant. This issues 
and the fact that hadrons were shown to have internal structure restricted the validity 
of the meson field theoretic approach to distances at which deviations from conceptual 
structureless hadrons are negligible. 

Furthermore, limitations to the perturbative approach also arise from divergences in 
meson meson interaction diagrams such as the one shown in Fig. (1. 2a). Such diver- 
gences are dealt through renormalization schemes by isolating a divergent term from a 
convergent part of the diagram and then introducing counter terms into the interaction 
Lagrangian. These counter terms generate extra diagrams that cancel the divergent term 
in the original diagram. 

The renormalized Lagrangian is then used to redefine parameters such as mass, cou- 
pling constant and wave function normalizations. Diagrams are then calculated using 
the redefined parameters in following Feynman rules. A scale dependence for instance 
is introduced in the renormalized or effective coupling constant g. The value of this 
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constant is now going to depend on the four momentum transfered at eacli vertex of a 
diagram. It also depends on its value g measured at some other scale. In the case of QED 
g = e = V47ra, with a « measured at almost zero momentum transfer (Q^ = -g^ = ^i^). 
The renormalized QED coupling constant is shown to deviate little from this value for 
a large range of energies. At the energy of the Z boson mass (90GeV), a « 1/127. These 
small values justify the use of perturbative expansions in powers of a in the study of 
numerous phenomena at many energy scales in QED. Nevertheless, in the asymptotic 
limit {-q^ 00), the value of the running coupling constant is given by: 

which is not very useful at energy scales in which a becomes greater than one, and of 
not use at all once -q- reaches the pole value (see e.g. Ref. [10]). For QED however, 
a{Q'^) increases very slowly from 1/137, and the pole on a{Q'^) is estimated to occur 
at energies-- W^^'^GeV, far beyond the Plank scale lO^'^GeV) and outside the energy 
domain where physical phenomena can be currently studied. Just as in QED, in a meson 
field theoretic approach to the strong interaction, treating divergent diagrams such as 
Fig. (I.2a) requires adding counter terms to the Lagrangian that generate diagrams such 
as Fig. (I.2b), and consequently a redefinition of the coupling constant (see Eq.(1.19)) 
that in the asymptotic limit just as Eq.(1.18), grows with and has a pole at some scale. 
However, the measured coupling constant at low energies(-- lOOsMeV) is already much 
larger than 1 {g'^/Air ^ 14), which leads to a near pole also at Q lOO's MeV, making an 
issue of the convergence of an expansion of the S matrix in series of Feynman diagrams. 



g^ , (1.19) 

Effective field Tlieories Due to the singularities discussed abuve, a meson field 
theoretic approach to the nn force is not well suited for phenomena at energy scales of 
few GeV. At these energies the internal structures of nucleons and mesons play more 
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a) b) 




i/4! ^ 



Figure 1.2: (a) Divergent fermion loop, (b) Counter term 

explicit roles in their interaction that cannot be accounted for by using hadronic degrees 
of freedom. More fundamental variables are sought for a theory of the strong interaction 
that can be extended to all scales where the strong force is dominant. These variables 
were later associated with the internal constituents of hadrons (quarks and gluons). 

A new Lagrangian, the QCD Lagrangian, is constructed from terms involving quark 
and gluonic fields, and from such a Lagrangian a consistent quantitative description 
of strong interaction phenomena is expected to emerge for all energy scales of inter- 
est. It turns out that in working which such a Lagrangian, there is not an analytic way 
of effectively computing observables at low energies. Instead, an effective field the- 
ory is developed constructing a new Lagrangian with the same global properties of the 
QCD Lagrangian, such as its invariance under chiral transformations in the massless 
quark limit. In this limit, these transformations flip the helicity of the quark fields. The 
spontaneous breaking of such symmetry gives rise to a set of massless pseudoscalar 
bosons which are identified with the lightest pseudoscalar mesons. It is also shown that 
the masses of these mesons come from the explicit breaking of the chiral symmetry 
when quark mass terms are added to the fermionic part of the QCD Lagrangian. Hence, 
hadronic degrees of freedom emerge through spontaneous and explicit symmetry break- 
ing of the free quark Lagrangian. 

An effective Lagrangian is accordingly written to group these effective hadronic 
degrees of freedom and observables are computed through the resulting Feynman dia- 
grams, which now are ordered in powers of (f^) where p is some external momentum, 
and --IGeV is known as the chiral scale. A perturbative approach is then valid as 
long as this parameter remains small, i.e., for energies of the order of hundreds MeV. 
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Thus, again, the use of the hadronic degrees of freedom in the analytical computation 
of observables based on a field theory is limited to a certain energy regime. 

While the quark mass term in the effective Lagrangian break chiral symmetry, mass 
differences between quark species breaks su{2) isospin invariance. A larger, S';7(3) flavor 
symmetry, is also broken by these mass differences. This symmetry is suggested by the 
discovery of hadrons with 'strange' decaying properties as it is described in the next 
section. 'Strange' pseudoscalar mesons named kaons compliment pions to form the 
group of eight light pseudoscalar mesons associated with the eight Goldstone bosons 
emerging from the spontaneous breaking of chiral SU{'i)R x S'C/(3)i symmetry. 

Strangeness 

In 1947 a new particle was discovered with the property of decaying into two hadrons 
p and TT at a rate much slower (half life ^ I0^"s) than typical decays through the strong 
interaction (lOaoS). It was named A° and since the net number of baryons or baryon 
number (s) was thought to be a conserved quantity regardless of the nature of an inter- 
action, the A" particle was concluded to be a baryon with a 'strange' decaying property 
accordingly named "strangeness". This property was also observed around the same 
time in the discovery of a meson named K which would decay into two pions, and in 
the discovery of a set of three baryons called Y,- and that 'slowly' decay into a 
nucleon and a pion. 

The half lives of A, s and k are typical in weak interaction processes such as beta 
and pion decays. Evidence of them taking part in strong interaction processes was 
found through the discovery of a resonance (E(1385)) in the reaction A' p A'V+tt". The 
Resonance e(1385) decays into a final Att system with a half live of ^ lo^^s which is 
characteristic of the strong interaction. This reaction and many others involving strange 
particles that followed were found to occur at these high rates. The common property 
of these reactions was that if one strange particle was present in the initial state, one 
strange particle will be present in the final state of the reaction as well. This is not the 
case for the weak decays explored earlier. 
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However, in some cases in which two or more strange particles were present in the 
initial state of a strong reaction, the final state did not necessarily have the same number 
of strange particles. This is the case for instance in the reaction 

K'^K" ^ TT-p. (1.20) 

To infer a conservation law related to strangeness, a new additive quantum number S 
was introduced. As the baryonic number is positive for baryons and negative for an- 
tibaryons, the strangeness quantum number can also take positive or negative values. 
While for nucleons 5 = 0, for K mesons (kaons) S = +i, and for A and s baryons 
5=1 with corresponding antihadrons having opposite strangeness. With this conven- 
tion the empirical evidence shows that reactions under the strong interactions conserve 
strangeness. 

Naturally, to preserve isospin as a good quantum number, values for both S and T 
should be assigned to the discovered strange baryons supported by the empirical ev- 
idence. Table I.l shows these assignments along with corresponding masses(ii/) and 
spin quantum numbers J. It also includes the heavier baryon named "cascade" (s) in 
reference to the two step weak strange decay (s ^ Att ^ Nt^t:) observed in its discovery 
that generated a cascade of particles. This behavior is understood if for s, 5 = -2 which 
is also consistent with S conservation in strong reactions involving s. 

SU{?>) Flavor and the Quark Model 

Table I.l suggests the classification of the shown hadrons in families of isospin (/) 
multiplets of a given strangeness (5). These multiplets coincide with dimensional rep- 
resentations of the 5;7(2) group of unitary transformations. The dimension of the rep- 
resentation would correspond to the number of particles in the multiplet. For instance, 
the family of pions is a three dimensional representation of SU{2), while the dimension 
of the s multiplet is 2. 

The conservation of strangeness hinted to a larger symmetry group for the strong 
interaction. The proposed group was the 5C/(3)-flavor symmetry group, of which SU{2) 
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Table I.l: Hadron Families 



was a subgroup. While the fundamental representation of SU{2) is two dimensional, 
which is realized in the nucleon and cascade duplets in Table I.l, the dimension of 
the fundamental representation for S'C/(3) is 3. And while the two states in the SU{2)i- 
fundamental representation are labeled u (or up for t = +1/2) and d (or down for t ^ 1/2), 
the states in the fundamental representation of su{i)f are labeled u, d and s. 

The baryons in Table I.l can be arranged in a eight-dimensional representation of 
S'C/(3) as represented in the in /a, y-plane (with y = s + b and h = t) Fig.1.3. Because of 
the mass differences between baryons, this symmetry is actually broken. This octet is 
part of the representations from arising from combining three SU{i) fundamental repre- 
sentations, 3«)3«)3=3(g)(6®3*)=10s®8Msffi8AfA©l (see e.g. [11]). Sub-indexes s, A and m 
stand for symmetric, antisymmetric and mixed under the interchange of flavors. 

The J 3/2 baryons were later presumed to form the 10-dimensional representation 
of S'C/(3) allowed in 'i®2>®2>. This picture was completed by the discovery of a baryon 
with5' = -3 {n-)m 1964. 

The pion triplet in Table I.l is built from a combination of the isospin doublet 2, 
and its conjugate representation which explains the presence of antimesons to complete 
the triplet,2®2*=3®l. In SU{i) with the same idea, the mesons should be arranged in 
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Figure 1.3: Baryon octet 
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Figure 1.4: Baryon decuplet 
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multiplets from 3(8)3*=8el. This structure corresponds to a 8-dimensional multiplet and 
a singlet. The mesons 77 and 77' were discovered in 1961, and join the mesons in Table 
I.l to complete the 8 and 1 representations as shown in Fig.1.5. Here, SU{3) is broken 
by the mass difference between the mesons. 

^ - ^ 

/ \ 
/ ^ 
/ __ \ 

71 / „ ^ 7t* 

^ 

\ -1/2 __ +1/2 ' T 

\ / 
\ __ / 

Figure 1.5: Meson nonet 

Despite the elegance of associating the hadron families with representations of a 
symmetry group, there was not a concrete explanation of why only some of the SU{3) 
representations were realized in nature. In particular, why there was not a hadron family 
associated with the fundamental representation of su{3) (i.e. 3) as it was the case for 

SU{2)j. 

Gell-Mann [12] and Zweig [13] suggested in 1964 that this family of particles ex- 
ists. They were named quarks, and it was further proposed that hadrons were actually 
systems of confined quarks from which they attain their intrinsic properties. The par- 
ticles in this fundamental SU{3) triplet are named after their flavors u, d and s, and as 
it's illustrated in Fig.1.6, their baryon numbers are B = 1/3, and B = 1/3 for their cor- 
responding antiparticles in the conjugate representation. Then in accordance with the 
group theory formalism of representations, baryons are now states of a system made up 
of three quarks 3«)3®3, and mesons are states of a system of a quark and and antiquark 
(from the conjugate of the fundamental representation) 3(g)3*. This way the hadrons' 
baryonic numbers result from the sum of their respective constituent quarks' baryonic 
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Figure 1.6: Quarks in SU{'i)j fundamental representation 



The liadrons' spin states should now come from the combination of the spin states 
of their two or three constituent quarks. In the language of representations, they are 
multiplets of dimensions allowed by 1®1 for mesons, and 2®1®2 for baryons, where 
2 is the fundamental representation of S'C/(2)-spin. The allowed multiplets are le3 for 
mesons and 2ms®'^ma®^ for baryons {Ms and Ma means symmetric and antisymmetric 
with respect to the exchange of the spins of the first two quarks) . This means that 
mesons can be grouped in multiplets of spin J = o and J = i, while baryons are grouped 
in two multiplets of spin J = i/2 and one with J = 3/2. 

These 5L/(2)-spin representations are then combined with the 5'C/(3)/ representations 
to fully characterize a hadronic {\h)) state in terms of its constituent quarks' states 
{\qq...)). For instance in the Bras-Kets notation a meson \m) or a baryon \b) with defined 
spin and flavor numbers is expanded in \qq) or \qqq) respectively according to: 

\H J, J^,I, T, Y)) = Cjtj..CJr^,T.,Y.,Y. X 

J^l,J^2,Tl,T2,l'i,Y2 

(1-21) 

and 

(1.22) 
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where 

Qi = q{Jzi,Ti,Yi), 

and the coefficients c are Clebsch-Gordan coefficients of the expansion. The labeHng 
of c^^^ is simpHfied understanding that there is an additional dependency from the la- 
beHng of the hadron's multiplet, i.e., a pair of numbers (p,q) analogous to J in SU{2) 
that define a representation in SU{i)f. The values of the C coefficients for both SU{2)- 
spin and S';7(3)-flavor also depend on the symmetry properties of the spin and flavor 
representations of \h), i.e., the sub indexes S, A, and M. Then, the combined represen- 
tation (SU(2),SU(3)) of which \h) is part of inherits a symmetry subindex as well from 
the product of the symmetries of the spin and flavor representations: For instance, the 

Symmetric: (5,5) {m,m) {A, A) 
Antisymmetric: (5, A) {m,m) (A, 5) 

Mixed: (5,M) (Af,5) (m,m) {m,A) {A,m) 

Table 1.2: Symmetry of (SU(2), SU(3)) representations. 

baryon decuplet is symmetric only in the spin-flavor combination (4,10). In fact, the 
only symmetric combined representations of ((2(g)2(g)2),(3(8)3(g)3)) are in (2,8) and (4,10), 
which correspond to the observed baryon J = i/2 octet and J = 3/2 decuplet. How- 
ever, if as assumed quarks are fermions, according to Pauli's principle they should only 
form antisymmetric states, therefore in principle a combination such as (4,10) should 
not exists for baryons. The contradiction is more explicit in baryons such as the A++ 
resonance. In its ground state |J = 3/2, = +3/2), all its constituent quarks are in the 
same state (|u(J = 1/2, = +1/2))), which is forbidden for identical fermions regardless 
of the exactness of 5C/(3)/. The anomaly was later resolved by introducing the concept 
of color charge which eventually led to the formulation of quantum chromodynamics 
(QCD), the fundamental theory of the strong interaction. 

The Parton Model 

Aside from the conflict between the idea of fermionic quarks and the exclusion prin- 
ciple, at the time there was not yet direct empirical evidence of their existence as real 
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particles. The phenomenological successes of the SU{3)f model in describing the hadron 
spectrum did not require quarks to be real particles. However, if in fact they exist and 
are the building blocks of hadrons, the strong interaction should be more fundamen- 
tally described through a theory of interactions between quarks. In this light, nucleons, 
and hadrons in general are composite systems of more elemental particles. The dis- 
crete nature of their structure was later revealed through experiments in deep inelastic 
scattering. 

From inelastic proton-proton scattering experiments it was observed that most of the 
hadrons emerging from the collision were vastly produced collinear with the collision 
axis. The strong suppression at large angles hints that the collision evolves in a very 
weakly interacting or dilute medium. 

Deep inelastic electron proton scattering experiments however showed that high en- 
ergy electrons have a large probability of scattering off protons with a significant en- 
ergy and momentum transfer. This meant that electrons were being deflected through 
the interaction with very localized concentrations of charge within the proton. In the 



parton model proposed by Feynman in 1969 (see Refs.[14] and [15]), the electrons in- 
teract incoherently with nearly structureless particle-like entities named partons inside 
the hadrons. The incoherence of this interaction is better understood in a reference 
frame where the hadron's longitudinal momentum P ^ oo accordingly named infinite 
momentum frame. Having that in the rest frame, partons inside the hadron interact 




Hadrons 



Figure 1.7: Deep inelastic scattering 
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by exchanging finite energy and momentum, their interaction times are finite as well. 
Then, when boosted to this "infinite momentum" frame the interaction time between 
partons is extremely dilated such that they appear basically free to the electromagnetic 
probe. 




par ton 



Figure 1.8: Parton picture of electron hadron interaction in DIS 



Scaling in DIS 

In the parton model developed by Bjorken and Paschos [16], the cross section for e+p^ 
e' + X is constructed from the cross section of an elementary electron parton elastic 
scattering (see Fig. (1.8)). Assuming that such parton of some specie (flavor, spin,...) i 
is a fermion (gj, the invariant amplitude for this elementary scattering is obtained using 
Feynman rules and for the spin averaged squared amplitude we have, 

in which eQ, is the electric charge of the struck parton g, , and the invariants s, i, and u are 
the Mandelstam variables for this subprocess: 

s=ip + kf i={k- k'f = 

u = {p- k'Y (1.24) 
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such that in the massless limit they satisfy s + £ + u = 0. Then, for the differential cross 
section we have. 



dt 167rs' ' 



t2 



(1.25) 



The invariants s, t, and u defined in Eq. (??) can be related to the Mandelstam 
variables s = {p + kf, t^{k- k'f and u^{P- k'f defined for the e+p^ e + x reaction by 
assuming that the parton's momentum is collinear with the proton's momentum. Then, 
defining ^ such that p = in the infinite momentum frame, we have that, 

s = £,s i=t, (1.26) 

in which ^ can be related to experimental variables. In the massless limit, for the parton 
we have that. 

Ok {p + qf = 2p-q + q'' = 2CP-q + q^ . 

Then with 



")2 



that in the lab reference frame corresponds to x = j^f^^- Equation (1.25) then corre- 
sponds to the differential cross section for an electron scattering off a quasi free parton 
of specie t concentrating x times the momentum of its parent proton (measured in the 
infinite momentum frame). 

Eq.(1.25) can now be written in terms of measurable s, Q'^, and x as. 



1+1-^ (1.28) 



dQ2 \ \ xs 

Thus, Eq.(1.28) is the contribution to e + p^ e + x from the scattering of the electron 
from any single parton of specie i and momentum fraction x. The total contribution of 
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partons of specie i and momentum fraction x is going to be weighed by the probability 
fi{x) of finding one of such parton in the proton. 

To obtain the differential cross section for the e + p e + X reaction, Eq.(1.28) is 
weighed by fi{x), summed over all parton species and integrated over all momentum 
fractions x to yield, 



i 

= jdxj:Mx)Qf-^i^i+{i-9iyy (1.29) 

On the other hand, from general principles the differential cross section of the inclusive 
reaction e + p^ e + x can be expressed as follows, 

d^^J '^"^ [yF.{^,Q)+^l-y- j --^ j , (1.30) 

in which Fj and Fa are inelastic structure functions, and y = In the large s limit 
(s >> Mp, y « ^), comparing Eq.(1.29) and Eq.(1.30) one obtains 

2xF,ix,Q') = F,ix,Q') =Y,Q>M^)- (1-31) 

i 

The first equality is well supported by experiments. It is known as the Callan-Gross 
relation, and is a consequence of the partons having spin i. Then, as it follows from 
Eq.(1.31), the inelastic structure functions, that can be extracted experimentally, are 
expected to become independent of within the parton model. This independence is 
referred as Bjorken scaling. Such prediction was confirmed in 1969 by the SLAC-MIT 
experiment in which the scaling behavior expected from Eq.(1.29) was confirmed for 
iGeV^ < Q2 < sGeW 

The experimental confirmation of the Bjorken scaling in ep eX scattering seemed 
to validate the picture of the nucleon to be a collection of almost free constituents (par- 
tons) as seen by a hard electromagnetic probe. Partons are then identified with the 
quarks introduced in the previous section, and the parton's specie i is associated with 
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Figure 1.9: Bjorken Scaling emerges from the the electron's incoherent scat- 
tering of the proton's partons in DIS 

the quark's flavor. However, the observed scaling contradicts a quantum field theory 
that would describe the interaction between these quarks. Since they do interact to be 
bound in nucleons, even at the smallest distances they are expected to exchange field 
quanta and deviations from a scaling pattern were to be expected instead. 

The above mentioned conflict was resolved in the 1970's by the discovery of field 
theories in which the interaction is characterized by asymptotic freedom (the strength 
of the interaction decreases as increases). Such behavior was discovered as a result 
of the regularization of diverging diagrams that, similar to QED and 0* field theories, in- 
troduces a dependence in the redefinition of the coupling constant. However, unlike 
QED and 0*, the value of the running coupling constant in these new theories can de- 
crease to zero as 00. This feature favors the main assumption of the parton model, 
namely that at very short distances the interaction between the partons in the nucleon 
is negligible. Nevertheless, deviations from scaling in DIS are still expected from a fi- 
nite probability of field quanta emission by the probed quark. They become significant 
for larger Q'^ or for smaller x than those reached by the SLAC-MIT experiment. Such 
violations of Bjorken scaling can be systematically studied through perturbation theory 
thanks to the smallness of the coupling constant. 

Quarks now become the fundamental fermions in a quantum field theory that de- 
scribes the strong force, while hadrons become multiquark systems bounded by their 
interaction under such force. As described in the following section such interaction not 
only binds quarks inside hadrons, but it also solves the apparent contradiction between 
the fermionic nature of quarks and the exclusion principle. 



28 



Color Charge and QCD 

In order to keep quarks as fermions and for them to obey Pauli's statistics, Greenberg 
proposed in 1964 [17] that a group of seemingly identical quarks should differ at least in 
a hidden quantum number for them to form a ground state hadron. This hidden quantum 
number was named 'color' or color charge, and its hidden property means that it cannot 
be observed, or that it cannot be probed in any observed particle. This proposition is 
in agreement with the fact that the SU{3)f fundamental representation is not realized 
in nature, if it's assumed that quarks have a nonzero color charge. Since quarks are 
'colored' objects, they have to be confined within colorless systems such as hadrons. 

In the language of representations, quarks are color non-singlets that combine to 
form the color singlet hadrons. A meson for instance is formed by a quark of color c 
and an antiquark of color c such that the total color of the system is . For a baryon 
where the color charges of its quarks are ci, ca, and cg we have that ci + C2 + cs = 0. 
Then, a quark within the baryon could carry any color charge c provided that the color 
charges of the other two balance it through the above relation. The smaller number of 
independent options for c is 3, thus a quark may be on one of three color states or on a 
linear combination of three states traditionally labeled R, G and B after red, green and blue 
in analogy with theory of colors where the combination of red, green and blue yields 
white. 

The special unitary group SU{3) is then the simpler choice to represent a quark color 
state. These states then belong to the fundamental representation of S';7(3)-color or 3. 
Now, a system of a colored quark and an anticolored antiquark can be represented by 
3(g)3*=le8 as discussed earlier regarding 5'c/(3)-flavor. Likewise, three colored quark 
states are represented by 3(8)3(8)3=1©8©8©10. Unlike S'f7(3)-flavor, only the singlets (1) 
are observed in nature. 

The e + e+ ^ Hadrons reaction 

As opposed to S'[/(3)-flavor, S'[/(3)-color is an exact symmetry. While a quark in a flavor 
state has specific mass and charge, it's color state is independent of any other intrin- 
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Figure 1. 10: Elementary particle picture of e e+ Hadrons 

sic property. The absence of color nonsinglets constrains direct evidence for this color 
degree of freedom hence direct evidence of this symmetry. Indirect evidence however 
is found in experiments in electron positron collisions where the energy from the anni- 
hilation creates multiple hadrons emerging from such collisions. The mechanism that 
creates these hadrons should start with the creation of a quark {q) and an antiquark (g) 
from a photon to which the e, e+ pair fused to. The q and q pair is created with enough 
center of mass (cm.) energy, that it can create more qq pairs from the vacuum ultimately 
hadronizing before reaching the detectors. As it is assumed in the parton model of DIS, 
for e~e+ hadrons the dynamics of the reaction is mainly dictated by the elementary 
subprocess e-e+ qtqi, whosc total cross section in the massless limit differs from that 
of e-e+ i-rfj.+ only by a electric charge factor Q, and by a color charge factor. Summing 
over all the possible color states of q in the colorless qq system, this color charge factor 
is 3, thus, 

f7(e"e+ qiqi) = 3(3iCr(e"e+ fJ-'fJ.'^), (1-32) 
and summing over all possible quark flavors i in order to obtain the total cross section 

for e~e+ — hadrons, 

a{e^e^ ^ hadrons) = y ^ ojer li^ji) 
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= 3^Q?a(e-e+^M-/x+). (1.33) 

i 

Then, the existence of a su{3) color degree of freedom implies that 

ji ^-(^-^^^ hadrons) ^^y^, (1.34) 

The ratio i? then depends on the number of quark flavors that can be produced in the 
elementary pair creation. The number depends on the center of mass energy of the e-e+ 
system. If the center of mass (cm.) energy equals the mass of the lightest bound state 
for a given quark flavor, that flavor enters in the sum in Eq.(1.34). Then, for different 
flavor families we have, 

3[(i)' + (|)V(|)^] =2 for u,d.s 
R=< 2 + 3(1)' =f for u,d,s,c (1-35) 

f+3(|)' =f for u,d,s,c,b 

A comparison of Eq.(1.35) with experimental data is illustrated in Fig.(I.l 1). The agree- 
ment improves as the center of mass energy increases far from the flavor threshold 
resonances regions. A description for such regions would involve more complex mech- 
anisms than the one discussed here. The results displayed in Fig.(I.l 1) strongly suggest 

the existence of a 3-degenerated color degree of freedom which consequently reinforces 
the parton-quark picture of hadronic structure. 

However, the main motivation for the study of hadronic structure, understanding the 
fundamental origin of the strong nuclear force, requires a description of how quarks, 
the elementary constituents of hadrons, interact. As described so far, evidence of the 
existence of these elementary particles is based on the notion of noninteracting mul- 
tiquark configurations, but as mentioned earlier such configurations are possible as an 
asymptotic limit in a field theory describing interactions between quarks. Just as QED 
describes interactions involving electrical charged particles and photons, this field the- 
ory describes interactions of color-charged quarks and the corresponding field quanta. 
Like photons in QED, these field quanta are bosons and they couple to corresponding 
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Figure I.l 1: World data on e+e hadrons cross sections (top) and R (botton) 
as a function of the center of mass energy {^f{s)) of the system [18]. 

color-charged currents, thus the name Quantum Chromodynamics (QCD) for this field 
theory. 

Gauge Field Theory of Quarks and Gluons 

The QCD field quanta are known as gluons, from the expectation that the force emerg- 
ing from quarks emitting/absorbing gluons binds (glues) them to form hadrons. Just as 
pions, the field quanta of the previously studied S't/(2)i,„,p,„ meson field theory, carried 
isospin charge and formed a 3-dimensional representation of SU{2)i,osviu, gluons carry 
color-charge and correspond to a 8-dimensional representation of SU{i)^oiov Unlike pi- 
ons which are pseudo-scalar (spin 0) fields, and like photons, gluons are vector (spin 
1) fields. Taking these considerations into account when writing representations of the 
elements of the quark gluon field theory, we have that for a quark field of a given flavor, 
taking into account color degrees of freedom. 



q{x) = 4>q{x) 
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(1.36) 
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while a gluon field is represented by, 



C^ix) = {Gt{x)} (1.37) 
in which (i=l,...8). Under a global SU{3)coior transformation, 

q{x) — > q'{x) = e 2 q 
q{x) q'{x) = q{x)e~^ , (1.38) 

in which {6*,} is a set of 8 real constant parameters. {A^} are the Gell-Mann matrices, the 
generators of the SU{3) group of transformations that as Pauli matrices in SU{2) satisfy 
the commutation relations 

and the normalization condition, 

tr {\,X,) = 25'' (1.40) 



The Lagrangian density for a noninteracting quark field, 

L ^ q{x) (ly - m) q{x) (1-41) 

is invariant under the transformations in Eq.(1.38). Quark-gluon dynamics arises how- 
ever if we require the Lagrangian of the theory to be invariant under local SU{3) trans- 
formations, U{x), i.e., if a space-time dependency is introduced in the set of parameters 
e in Eq.(1.38). Then, 

q{x) — >■ q'{x) — U{x)q{x) — e 2 q 
q{x) q'{x) = q{x)U-\x) = qe"^^ . (1.42) 

The Lagrangian density L in Eq.(L41) is not invariant under these transformations. This 
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is because of the differential operator v acting on the parameters 0{x) which generates 
extra terms in the transformed Lagrangian. However, if a new Lagrangian is built from 
Eq.(1.41) by replacing v with 

D,^^,~^g^Giix), (1.43) 

thus obtaining, 

Lj = qix){i]p-m)q{x). (1.44) 

Having that under the local SU{3) transformation of Eq.(1.42), 

A'G;(x) -> Uix)X'G'^ix)U-\x) ~ '-[d,U{x)]U-\x), (1.45) 

the new Lagrangian density Lj in Eq.(1.44) is then invariant under local SU{3) trans- 
formations, since Eq.(1.45) guarantees that Dq{x), known as the covariant derivative of 
q{x), transforms by the same rule as q{x). Equation (1.45) is known as a SU{3) gauge 
transformation on G^ix), and (x) are known as gauge fields which then are chosen to 
represent gluons. Thus, Dq introduces a quark-gluon interaction term (q-/\Gq) in a local 
SU{3) invariant Lagrangian. A special case, U{1) gauge transformations correspond to 
the local gauge symmetry of the QED Lagrangian that gives rise to the photon fermion 
interaction'. 

As it stands, Eq.(L44) accounts for the dynamics of the quark field. However just as 
in QED G"" fields by themselves should contribute to the total Lagrangian. They do so 
through the gauge invariant term, 

^^^-T^'^"" (1-46) 



in which. 



F;.-^d,-Gl-d^-Gl-tg 



GlGi. (L47) 



'See e.g. Refs, [9] and [19] for a complete introduction of gauge field theories and their quantization 
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Because of the Abelian nature of U{i), in QED the third term in Eq.(I-47) is zero. In 
QED, Ff,^ is recognized as the electromagnetic tensor, and classically the equations of 
motion obtained from Eq.(1.46) yield the homogeneous or free Maxwell equations for 
electric and magnetic fields. 

For QCD on the other hand, from Eq.(1.39), the commutator term in Eq.(1.47) is not 
zero and f^^ is not a linear function of G„. Consequently, Eq.(1.46) introduces gGGd^G 
and g'^GGGG terms in the total Lagrangian which now takes the following form. 



These three-gluon and four-gluon fields terms give then rise to elementary interactions 
between gluons at the lowest an the next to lowest order expansions in the coupling con- 
stant (a, = 1^) in the perturbative approach to QCD's elementary processes. Fig. (1. 12) 
illustrates the building blocks for such interactions according to Eqs.(1.47) and (1.48). 
In contrast, from a QED Lagrangian, only an analogous to the first vertex in Fig. (1. 12) 
is present. 



Figure 1.12: Elementary interactions in Lqcd- (b) and (c) are introduced by 
£) (Eq.(1.46) from the nonabelian nature of SU{3)^oior 

Asymptotic Freedom and pQCD 

In perturbation theory, the expansion of QCD observables in Feynman diagrams in- 
volves loop contributions such as those to the gluon propagator displayed in Fig. (1.13). 
Factors corresponding to the tree-gluon vertex in Fig. (1. 12) generate the gluon loop di- 
agram correction to the gluons Green's function. As it is the case for QED, for QCD, 
these contributions diverge individually, with such divergences emerging from the in- 
tegration over all momentum space of the loop momentum variable not fixed by 4- 



L = q{x) [lYd, - m) q{x) + gq(x)Y -Gl{x)q{x) - -F^^Ff ■ 



(1.48) 




(b) 



(c) 
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momentum conservation. Thus, as discussed earlier, through renormalization schemes 
these divergences can be absorbed in the redefinition of the parameters of the theory 
such that the methodology of perturbative expansions in Feynman diagrams can be 
safely applied. 





Figure 1.13: Contributions to the gluon propagator. The loop diagrams are 
divergent, and their regularization leads to the dependence in the redefi- 
nition of the coupling constant (Eq.(1.49)). 



A scale dependence is then introduced in the redefinition of the QCD's coupling 
constant. In the large limit it's found that at the one loop correction. 



l + £^(33-2n,)/o.g^ 



(1-49) 



in which a{^i^) is the value of the coupling constant measured at some scale fi^, and iif is 
the number of quark flavors in the theory. With the substitution. 



-12n 



(33-2n/)a.(M=) 



(1.50) 



Eq.(1.49) can be rewritten. 



127r 



(1-51) 



Unlike a for QED in Eq.(1.18), provided that Uf < 16, a, for QCD in Eq.(1.51) decreases 
towards zero as ^ cx). Thus, the interaction mediated by gluon exchanges grows 
weaker as the distance between the interacting particles becomes smaller. This property 
is known as Asymptotic Freedom and its discovery by David Gross, Frank Wilczek [20] 
and David Politzer [21] in the 1970's validated within the framework of a field theory 
the main assumption on which the earlier parton model successfully described phe- 
nomena such as scaling in ep^ eX. Namely, that to a electromagnetic probe with high 
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resolution (large Q^) quarks deep inside the nucleon behave basically as free particles. 

Quantum chromodynamics then provides an important feature in understanding prop- 
erties of hadronic structure. It also provides a methodology for examining corrections 
to results predicted within the parton model approximation. Equation (1.49) is valid 
for momentum energy scales at which < i. For such momentum-energy regime, 
the use of perturbation theory at the elementary particle level is well justified. There- 
fore, corrections to results from the parton model arise from calculable diagrams of 
gluon emissions by the struck parton. For instance, the leading order correction to the 
structure function F2 is of order aa, and it comes from the superposition of the gluon 
emission diagrams in Fig. (1. 14). The quark emerging from the second or third diagrams 




(a) (b) (c) 



Figure 1.14: (a) quark-photon interaction controlling the dependence in ep 
DIS. (b) and (c) gluon emission diagrams introducing a dependence in F2 
and braking the scaling in Eq.(1.31). 

carries a fraction z of the momentum of the initial parton. Unlike the situation in the 
first diagram, this quark can also carry transverse momentum relative to the virtual pho- 
ton's in the infinite momentum frame. This transverse momentum is balanced by the 
transverse momentum of the emitted gluon to recover the zero transversity of the initial 
quark. Therefore, The elementary cross sections of these 7*9 gq subprocesses are de- 
pendent on this z momentum fraction, and on the transverse momentum that the quark 
acquires after emitting the gluon. After integrating over all possible values allowed for 
this transverse momentum, we have that; 

= aQ^d,^P,,{z)log{^^ (1.52) 

where do = Air'^a/s, and Pqq{z) = |^ is called splitting function and represents the prob- 
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ability of a quark emitting a gluon that carries a i - z fraction of the parent quark's 
momentum. 

For the second line in Eq.(1.52) it is used that 

and a ^ cut off is introduced to regularize the divergence when o. Thus, as opposed 
to the parton model approximation, QCD predicted gluon emissions introduced a 
dependence in the elementary parton electron interaction. As a result, the structure 
functions are corrected as follows, 

= T.< (i(-) + Q")) ' (1-53) 

9 

with 

A„,<3^).g,,,(2!)£^!,(„,P„(i). (1.54) 

The parton model limit is met when the above fluctuation is negligible.Then, q{;x) = f^{x) 
and Eq.(1.31) is recovered from the general form. 



^^%^ = E^^9(x,g^). (1.55) 



Experimental data on F2 is shown in Fig. (1.15). Deviations from scaling behavior of Fa 
are observed across values of x being more dramatic at small x at which as increases 
the structure function increases as well. The structure function F2 scales approximately 
in the region 0.2 < x < 0.3, while for larger x it slowly decreases with increasing 
Equation (1.54) can be rewritten to formally exhibit the evolution of the distribution 
functions q{x,Q'^) (see e.g. Ref. [22]), 

' -M^. Q') = ^ f -i{y, Q')P,, i -] ■ (1.56) 



dlogQ^ ' ' 2ti 

Thus, from the preexisting knowledge of an experimentally extracted distribution func- 
tion for some fixed Ql, through Eq.(1.56), q{x,Q'^) can be computed for any large 
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Figure 1.15: behavior of F2. Figure taken from Ref. [23] 

value of Q^. The evolution from this fixed Ql is dominated by the dependence of a,, 
i.e., q{x,Q'^) evolves logarithmically with Q^ and thus, for large the deviation from 
scaling is rather a subtle effect as can be seen in Fig. (1.15). 

Equation (1.56) is further corrected by considering that the struck quark may have 
come from a quark- antiquark splitting of a prexisting gluon in the proton. The con- 
tribution from pair creation to q{x, Q^) depends on a preexisting distribution of gluons 
(g{x,Q^)). Then, Eq.(1.56) becomes 



dlogQ^ 2n y 



(1.51) 



in which the splitting function P^g = " represents the probability for the struck 

quark to carry a z fraction of the momentum of the initial gluon. The function Pqg{z) can 
be obtained from the cross section of the subprocess 7*3 qq in the same way that 
is obtained from the cross section of the subprocess 7*9 gq. 

Solving Eq.(1.57) requires knowing g{x,Q'^). Because gluon sources in the nucleon 
are both quark and preexistent gluons an evolution equation analogous to Eq.(1.57) can 
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be derived through the same methodology for g{x, Q"^) yielding 

in which the splitting functions Pg^iz) and Pgg{z) can be obtained from elementary dia- 
grams of a gluon emitted with a z momentum fraction of parent quark or a parent gluon 
respectively. In all we have for ^ < i the following splitting functions 

4 1 + 

^ ^^^^ 

Pgg{z) = 6(^l-i + ^+.(l-Z)). (1.59) 

The equations (1.57 and 1.58) are known as the DGLAP (Dokshitzer [24],Gribov and 
Lipatov [25], and Altarelli-Parisi [26]) QCD evolution equations, and they describe the 
internal longitudinal momentum structure of the nucleon as it evolves in the large 
region. 

Going back to Eq.(1.29), to account for the evolution of parton distributions with 
the cross section for ep eX is corrected to 

^ = j dx'£q{x,Q-)^-^^2^0^. (1.60) 

The evolution equations show that the factorization implied in the parton model is val- 
idated by the slow evolution of q{x,Q'^) at large In contrast, as seen in Eq.(1.28) the 
elementary subprocess's cross section has a large dependence on Thus, for these 
kinematics the short distance physics of the reaction is mostly contained in this elemen- 
tary subprocess' dynamics while the long distance effects are isolated in factors such as 
the parton distributions q{:x) in Eq.(1.60). 

These factors then describe the dynamics of partons (quarks and gluons) as they 
interact to form hadrons. The full description of these factors however is beyond the 
reach of pQCD, because the behavior of the coupling constant is not longer described 
by Eq.(1.51). How this effective coupling behaves at long distances is not well under- 
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stood; Eq.(1.51) is no longer valid for < a « 200MeV, at which values a« > i, making 
perturbative expansions in Feynman diagrams no longer convergent. 

The strength of the interaction between quarks is believed to grow as the distance be- 
tween them increases which leads ultimately to the confinement of quarks into hadrons. 
Although the dynamics of this confinement cannot be described by pQCD, general fea- 
tures of hadrons' structure are expected to be described through the use of the QCD 
Lagrangian in nonperturbative approaches of which Lattice QCD is the most widely 
accepted. 

In addition to explaining the deviation from the scaling pattern of structure functions 
in DIS, from the behavior of a, dictated by Eq.(1.5 1), QCD also introduces corrections to 
the results in Eq.(1.35) for e, e+ annihilation into hadrons. Accounting for gluon emission 
by one of the quarks of the created quark antiquark pair introduces an additional energy 
dependency in R, 



Again QCD predicts a scaling violation of order logQ'^ for the otherwise independent 
behavior of R. Current experimental data however do not reach a region where this 
difference is observable as seen in Fig.(I.l 1), thus making the effects of gluon emission 
negligible. 

If now within this approximation the reaction ee+ hx (with h a hadron of a given 
specie) is considered, its cross section is also predicted to scale according to 



in which D''^{z) is known as a fragmentation function, and similarly to the parton dis- 
tribution functions in ep eX, it represents the probability that the hadron h carries a 
fraction z = of the energy of the parent quark (g) in the quark antiquark pair from 
which it is produced. Likewise, the scaling predicted by Eq.(1.62) is broken by a logQ'^ 
evolution of the fragmentation functions arising from gluon emission. Another source 




(L61) 



a dz 



[D^,{z) + D'l{z)] 



(1.62) 



41 



of scale breaking arises from the crossing of the charm threshold at which the heavy 
cc system is created almost at rest and subsequently decaying weakly into many low z 
hadrons. 

Drell-Yan processes 

Hadron-hadron interactions can be naturally approached through the same method- 
ology applied to the electroproduction processes. One standard example is the reaction 
■p-p l-l+x in which a lepton-antilepton pair is created at large transverse momentum 
in a proton-proton collision. It is known as the Drell-Yan process and proceeds at the 
parton level as shown in fig. (1. 16). The Drell-Yan cross section within the quark-parton 

P 



P 

Figure 1.16: Drell-Yan process, pp 1^1+ in the QPM 

model (QPM) is, 

^(pp ^ 1-1+ X) = J2j'^'^J dyq{x)q{y)-^{qq ^ /+r), (1.63) 

in which 

^{qq ^ m-) = el^SiQ^ - [xp, + yp.f), (1.64) 

with X and y being the fraction of the respective parent nucleons momenta carried by 
the quark and antiquark participating in the annihilation subprocess, and being the 
invariant mass of the 1^1+ system. 

For large s and we have that = xys. Making this condition explicit in Eq.(1.63) 
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yields 

±^{pp^l-l+X)^^Y.< Jd^l dyq{x)q{y)s(l^xy^Y (1.65) 

q 

which suggest the introduction of a scaling variable w = xy ^ ^ resulting in the scaling 
law, 

Q'^,{pp^l-l+X) = y{w) (1.66) 

This scaling is well satisfied by experimental data from which can be extracted. 
Also, just as it happens in DIS and e e+ hx the scaling predicted in Eq.(1.66) is broken 
by logQ'^ coming from corrections to fig. (1. 16) including quark gluon interactions. 

Cross sections and other observables of numerous processes involving hadrons can 
be studied by making use of parton distribution functions and fragmentation functions. 
The universality of these functions further allows the study of QCD subprocesses con- 
trolling a" order corrections in reactions such as those shown in fig. (1. 17) 




a) b) 

Figure 1.17: Typical inclusive DIS processes involving gluons in the corre- 
sponding parton subprocesses. 

These studies are made possible by the large kinematic variables involved in the 
reactions (s, Q,pt) that allow the use of pQCD in describing the dynamics of the embed- 
ded subprocess. The kinematic ranges in which this description is applied is known as 
the hard kinematic regime, and the controlling subprocesses are known as hard subpro- 
cesses. Consequently, the short distance factors such as the cross sections of the hard 
subprocesses and the long distance factors such as parton distribution functions (/) and 
fragmentation functions (d'') are known as hard factors and soft factors respectively. 
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Asymptotic freedom has allowed the experimental study of QCD at the elementary 
particle level. Because of confinement, in experiments, the reactions studied do not pro- 
duce free quarks or gluons which are the fields associated with the degrees of freedom 
in the QCD interaction Lagrangian. But because at large energy and momentum scales 
the QCD's coupling constant weakens and changes slowly, the elementary particle pro- 
cesses at the core of these reactions dominate the kinematic dependence of the reaction 
and thus the short distance dynamics, making the study of the reaction an almost direct 
study of the hard subprocess. Corrections to this picture are also a test of QCD, since in 
the hard kinematics, such corrections come from the calculable pQCD evolution of the 
soft factors. 

The kind of reactions considered so far are of the form ab {cd...)x, i.e., in exper- 
iments only particles (cd...) in the final state are detected, while x sums up all those 
that are not. The final states of these reactions then are not specified. The cross sec- 
tions given are obtained including all possible final states, thus these kind of reaction 
are known as inclusive or semi-inclusive reactions. The extensive study of inclusive 
reactions has reinforced QCD as the theory of the strong interaction that arises from the 
interactions between quarks and gluons. However, precisely because all the hadronic 
states are not completely resolved a description of the strong force can't be extracted 
through QCD from inclusive processes alone. 

In what follows we focus on a more constrained kind of processes known as exclu- 
sive reactions. Such processes are of the form ab [cde...), thus all the particles in the 
final state are detected. The role of QCD degrees of freedom in hadronic structure and 
in hadronic interactions is more carefully investigated here because of the constrains of 
constructing the final hadron states from the scattered elementary particles. This diffi- 
culty is of course avoided through the parton model in the previously studied inclusive 
reactions. 
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1.3 Hard Exclusive Processes 

In this section we introduce a methodology and a formalism that is closely followed in 
the study of the processes of interest in this dissertation. 

With the purpose of developing a QCD description of interactions between hadrons 
(as bound states of quark systems) and at the same time gaining further insight into 
their quark structure, the studies of exclusive reactions at large momentum transfer 
have attracted particular interest. These hard processes facilitate the study of the short 
distance regime of the strong nuclear force in reactions such as hard proton proton (pp) 
and proton neutron (pn) elastic scattering. The characteristic kinematic region for hard 
exclusive processes is defined by large kinematic variables, i.e. s, —t, —u >> m^, and -, 
and - fixed. 

s 

To illustrate the extent through which the analysis of hard exclusive scattering pro- 
cesses potentially probes the quark-gluon dynamics that underlies hadronic processes 
we first consider the reaction 7*A^ N. Fig.(1.18) illustrates how the reaction may 
proceeds at the constituent quark level. 




Figure 1.18: ep elastic scattering factorization according to Eq.(1.67). 

The large momentum transfer by the virtual photon to the nucleon is distributed 
among the constituents through the hard subprocess Th. Unlike the situation in inclusive 
DIS, here, coherence is required to form the outgoing nucleon, thus all the constituents 
are involved in the hard subprocess. Fig.(1.18) also suggests the factorization of the hard 
scattering amplitude (Th) from the soft factors corresponding to the constituent distri- 
butions of the incoming (ip) and outgoing (ip*) nucleons. These factors are convoluted 
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in the matrix element [27] of tlie diagram in Fig. (1. 18), 

= ^ J[dx] J[dy](p''{yi)TH{xi,y„Q)ip{x,). (1-67) 

This matrix element in the infinite momentum frame corresponds to the nucleons's mag- 
netic form factor which at large -t = and s» ~t can be extracted from experiments 
using the relation 

^(ep->ep) = 2™^^. (1.68) 

In Eq.(1.67), 

[da;] = (5(1 — ^^a;^) JJ^dxi, 

in which Xi (yO is the fraction of the initial (final) nucleon's light cone momentum car- 
ried by its constituent i. The tp{x) {(p'iy)) is the amplitude for finding in the incoming 
(outgoing) nucleon the configuration of near on- shell partons entering (leaving) the hard 
Th blob. 

Quark wave functions of hadrons 

The nucleon shown in Fig.(1.18) is represented by a system of three partons, three 
quarks that in this case is the minimum number of quarks that are needed to reconstruct 
the helicity and isospin of the nucleon. The nucleon in Fig. (1. 18) illustrates what is 
known as the minimal component of the Fock expansion of the nucleon wave function. 
In Fock space the nucleon is expanded in states with a definite number of constituents, 
i.e, 

ip^qqq + qqqg + qqqqq + ... 
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Labeling each component by N (number of constituents) the nucleon wave function in 
terms of its quark and gluon constituents is given by [28] 

\p,h)='Y] [[dx][cPkT] — ]. ipix,, kiT, hi)\xi,hT, h,), (1-69) 

J ^/Xl...XN 

in which 

CJV \ 
j=i / 1=1 

where \xi,kiT, h^) represents an individual constituent state. In the light cone gauge [27], 
the amplitudes f{x) in Eq.(1.67) is related to ipixi^hr, h,) by 

...,a;jv) = [d'^kT]tp{Xi,kiT,hi,fi), (1-70) 



in which it is understood that the additional helicity and flavor expansion coefficients 
such as those in Eq.(1.22) are contained in ip{x). 

Scaling Laws 

The hard subprocess in Fig.(I.18)involves the far off shell (k'^ ^ Q'^ » ni^) propaga- 
tion of the incoming nucleon's constituents into the outgoing constituents that form the 
final nucleon. A typical diagram contributing to the eqqq eqqq subprocess is shown 
in Fig. (1. 19). It is known as a minimally connected diagram and guarantees that the 




X ,p "//^"^ ^ y , (p+q) 



x.p 












(y 3^ j)p +y 



y Jp+q) 



y 3(p+q) 



Figure 1.19: Typical minimally connected diagram contributing to T„ in eN 
elastic scattering. 



momentum q transfered by the electron probe steers all the quarks collinearly in the di- 
rection of the outgoing nucleon. At ^ oo this kind of diagrams dominate in Th since 
they represent the lowest order interaction in among the constituent quarks. The 
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asymptotic dependence of Th can then be obtained from pQCD. Having that the quarks 
in Th interact through vector exchange, each pair of external lines contributes a factor 
Q. This is seen for instance in the vertex element, 

u{h' = h,k + qh^'uih, k)\Q2^^ cx Q (1.71) 

in which u are quark helicity spinors. \f h' ^ h this factor behaves as ^. Each fermion 
propagator contributes a factor of ^ and each vector boson propagator contributes a 
factor - ^, thus resulting in 

Th a{-^fTH{x,y) 

which from Eq.(I.67)and Eq.(1.68) leads to 

(1.72) 

or for fixed cm. angle of ep ep scattering (t/s fixed, t, s oo, and s » -t), 

|(e,^ep)^(41!/(t/.). (1.73) 
at s° 

The contributions to Th from diagrams with additional constituent lines taking part in 
the hard subprocess fall off as powers of ^ faster as compared to the contribution from 
Fig.(1.18) that contains the minimal number of constituents. Hence, the minimal Fock 
component of the nucleon contributes to the dominant term of the scattering amplitude. 

Through the same rules that lead to Eq.(1.73), for the en en scattering reaction we 
have that, ^ ^, while it is known that for e/x ^ e^, ^ ^. Thus in general we have 
that for eH eH, 

sV+-«^"""^ /„(-), (1.74) 
at s 

in which Nh is the number of minimal constituents in particle H. 

A hadron hadron scattering amplitude in the hard kinematic regime can be factorized 
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T 



Figure 1.20: a^a^ elastic scattering controlled by the expansion of con- 
nected constituent diagrams. 

similarly to Eq.(1.67) in the following form: 

'Mab^cd= J Y\. [dXi]ip*(y{Xc)(pl,{Xc)TH{x,,S,t)ipA{Xc)(pB{Xc). (1-75) 

i=a,b,c,d 

Then in the above equation Th is dominated by connected diagrams involving con- 
stituents from the minimal Fock components of the interacting hadrons such as those 
shown in Fig. (1.75) for nucleon nucleon (nn) elastic scattering. Although there are at 
least thousands of such diagrams contributing to Th, they contain the same number of 
fermion and vector boson external and internal lines. Then as before, at fixed angle in 
the asymptotic limit there is a factor per each pair of quark external lines, a ^ per 
each quark internal line, and a ^ factor per each vector gluon internal line. Then for 
NN elastic scattering, Th ^ j^T{x,e^,^), and for the nn elastic scattering amplitude one 
obtains, 

M^^^^^ ~ -M{t/s), (1.76) 

and 



dt 



s 



\M\ 



■j^fNNit/s). (1.77) 



In general, the asymptotic energy dependence of the invariant amplitudes and of the 
cross sections for exclusive processes is given respectively by the forms (see Refs. [29], 
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and [30]), 




(1.78) 



and as in Eq.(1.77), 




(1-79) 



in which uh is the minimal number of constituents of particle H (h = ABCD...) taking 
part in the hard subprocess, while f{t/s) is an angular function that becomes invariant 
under a change of energy scale. The dimensions of the scale invariant f{t/s) are given by 
the number of constituents taking part of the elementary hard subprocesses (for instance, 
fit/s) is dimensionless for efj. e^i) consequently determining the power fall-off of the 
energy (s) dependence of the differential cross section. 

Evidence of a power law fall-off in energy distributions at fixed angle has been 
observed in experiments for several hard processes involving hadrons, with many of 
them closely fitting the scaling behavior predicted by Eq. (1.79) (see e.g. Refs.[31, 32, 
34, 33, 35, 36]). Empirical agreement with Eq. (1.79) however does not necessarily 
implies the perturbative approach used here to arrive at this scaling law. Nonetheless, 
the correlation between the number of constituents and the energy dependence of the 
exclusive reactions evidenced by this result hints to the dominance of minimal Fock 
component of the partonic wave functions of the hadrons participating in the reactions. 

Quark Interchange 

In the constituent picture, the processes through which hadrons in an exclusive reaction 
interact include the mechanisms shown in Fig. (1. 21). Unlike the situation involving 
reactions at low energy and low momentum transfer, in which the reactions proceed 
through the exchange of mesons, as illustrated in Fig. (1. 21) hard exclusive processes 
involving hadrons proceed through the exchange of quarks or gluons. QCD excludes 
the mechanism of one gluon exchange in hadron-hadron exclusive scattering because 
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hadrons are color singlets (of zero color charge), but the exchange of one gluon would 
result in a net transferring of color charge from one hadron to the other producing 
two non-singlet particles in the final state. The latter scenario is not observed in na- 
ture, therefore gluon exchange mechanisms in hadron-hadron scattering involve the ex- 
change of two gluons, as shown in Fig.(I.21a), or more such that it can ensure a net 
transfer of zero color charge between the interacting hadrons. The other mechanisms of 
interaction shown in Fig.(1.21) (b) and (c), correspond to quark- antiquark annihilation 
and quark interchange. The quark- antiquark annihilation mechanisms are considered 
in scattering reactions in which one of the hadrons participating contains antiquarks 
of corresponding quarks from the second hadron. The reaction then proceeds through 
the annihilation of a quark from one hadron with an antiquark from the other hadron. 
Quark interchange mechanisms are considered for instance in elastic scattering reac- 
tions in which both hadrons contain quarks of common flavor. It has been argued that 




Figure 1.21: QCD description of exclusive hadron-hadron scattering through 
(a) gluon exchange, (b) quark- antiquark annihilation and (c) quark inter- 
change. 

the quark interchange mechanism dominates the exclusive scattering of hadrons that 
posses quarks of common flavor among their constituents [109, 108]. 

The dominance of the quark interchange mechanism has been experimentally con- 
cluded when comparing for instance data on large angle scattering for proton proton 
(pp) and proton antiproton (pp) elastic reactions. It was found that at large angle, the pp 
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cross sections largely dominates pp [37]. At 90° cm angle of scattering and lOGeV of 
beam momentum, the cross section of pp elastic scattering is more than 24 times larger 
than that of pp [37]. If gluon exchange mechanisms ( Fig. (1.21 a)) were the dominant 
mechanisms of interaction, the cross sections for both reaction would be expected to be 
similar because these mechanisms do not depend on the quark- antiquark composition of 
the interacting baryons. Because there are no quarks of common flavor between protons 
and antiprotons, quark interchange mechanisms do not contribute in pp elastic scattering 
while they do in pp elastic scattering. Having that the cross section of pp elastic scatter- 
ing is much larger than that of pp, it is concluded that quark interchange dominates pp 
elastic scattering [37]. 



Figure 1.22: Scattering channels in the quark exchange mechanism of 
baryon-baryon (bb) exclusive scattering. 

Quark exchange in baryon baryon exclusive scattering BaBb Be Bo proceeds through 
either of the two channels illustrated in Fig.(I.22)(a) and (b). The scenario (a) corre- 
sponds to the exchanging quarks scattering in the u channel while the residual system 
scatters in the t channel, while in (b) the residual system scatters in the u channel while 
the 'exchanging' quarks scatter in the t channel. As it follows from the previous sec- 
tion, at fixed angle, the energy dependence of these two contributions is the same. The 
angular dependence on the other hand goes as CtF{t/s) for Fig.(I.22)(a), and c^f{u/s) 
for Fig.(I.22)(b). While F comes from the interaction among constituents, the coeffi- 
cients Ct, and C„ come from the quark wave functions of the baryons participating in 
the interaction. Appendix B illustrates how these coefficients are obtained from a given 
flavor-helicity expansion of the quark wave functions of the interacting baryons. In such 
approach, the residual system consists of the two remaining quarks in the minimal Fock 
component of the wave function of the interacting baryon. 




(a) 



(b) 
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1.4 Hard Processes Involving Nuclei 

The descriptions of nuclear processes based on QCD is met with the empirical difficulty 
of quarks and gluons not being experimentally observable particles. Up to certain scale 
they are instead 'hidden' within the hadrons that mediate and take part of the interaction 
binding nucleons in the nucleus. The experimental constrains and the computational 
challenges brought on by confinement made the pursuing of a truly elementary particle 
description of the nucleus very impractical when compared to the alternative effective 
field theoretical frameworks. 

Nonetheless, reactions involving nuclei in the hard kinematic regime (at energy 
scales in which pQCD may be applicable) can probe the quark-gluon dynamics of small 
size nucleon-nucleon (nn) configurations ( small enough such that their overlapping 
dynamics takes the nucleons' internal structure into account). Such configurations are 
expected to be naturally present in dense nuclear matter. 

Small NN configurations can also be probed in light nuclei for instance through ex- 
clusive reactions in which nuclei absorb enough energy to break into nucleons emerg- 
ing with large relative transverse momenta. This kind of reactions is the focus of the 
research work detailed in chapters III and IV and Refs.[98] and [114] in which the ex- 
plicit role of subnucleonic degrees of freedom is investigated in selected hard processes 
involving light nuclei. The following section lists the main features of characteristic 
approaches developed to incorporate QCD in the description of these processes. 

Hard breakup of A^A^ system in nuclei 

In the research work described in this dissertation, we focused on the reaction 7 + A ^ 
{NN)+{A-2) in which (NN) is produced at large center of mass angle. Two-body breakup 
reactions involving nuclei at high momentum and energy transfer play an important role 
in studies of nuclear QCD. The uniqueness of these processes is in the effectiveness by 
which large values of invariant energy are produced at rather moderate values of beam 
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energy. For a photodisintegration process of an nn system we have that, 

s^MN ~ 4m^ + AE^rriM, (1.80) 

in which the produced invariant energy grows with the energy of the probe twice as 
fast as compared, for example, to hard processes involving two protons, in which case 
% + 2E ■ mjv. As it follows from Eq.(1.80), already at photon energies of 2 GeV 
the produced invariant mass on one nucleon, M ^ y^^^f^, exceeds the threshold at which 
deep-inelastic processes become important, M > 2 GeV. 

Combining the above property with a requirement that the momentum transfer in the 
reaction exceeds the masses of the particles involved in the scattering ( 
order for the reaction to reach the hard scattering kinematic regime. In this regime it is 
expected that only the minimal Fock components dominate in the wave function of the 
particles involved in the scattering. Assuming that all the constituents of minimal Fock 
component participate in a hard scattering, it is expected that the energy dependence 
of the reaction follows the constituent counting rule of Eqs.(1.78) and (1.79) [29, 30] . 
These predictions have been confirmed for a wide variety of hard processes involving 
leptons and hadrons (see e.g. Refs.[31, 32, 34, 33, 35, 36]). 

One of the most interesting aspects of the constituent-counting rule is that its ap- 
plication allows us to check the onset of quark degrees of freedom in hard reactions 
involving nuclei [51, 52]. This is essential for probing the quark-gluon structure of nu- 
clei. For example, if quarks are involved in hard photodisintegration of the deuteron 
then according to Eqs.(1.78) and (1.79) one expects that ^ s " [51]. 

During the last decade there were several experiments in which 90" cm. photodis- 
integration of the deuteron had been studied at high photon energies [53, 54, 55, 56, 
57, 58, 59, 60]. These experiments clearly demonstrated the onset of s " scaling for 
the differential cross section at 90" cm., starting at -B^ > i GeV. Also, the polarization 
measurements[57, 60] were generally in agreement with the prediction of the helicity 
conservation - a precursor of the dominance of the mechanism of hard gluon exchange 
involving quarks. 
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Even though two-body scattering experiments demonstrate clearly an onset of quark 
degrees of freedom in the reaction, they do not affirm the onset of the perturbative QCD 
(pQCD) regime. Indeed it has been argued that the validity of constituent-counting 
rule does not necessarily lead to the validity of pQCD(see e.g. Refs.[61, 62]). In 
several measurements in which the constituent quark rule works pQCD still underes- 
timates the observed cross sections sometimes by several orders of magnitude (see e.g. 
Refs.[63, 64]). The latter may indicate a substantial contribution because of nonper- 
turbative effects although one still may expect sizable contributions from pQCD due 
to generally unaccounted hidden color components in the hadronic and nuclear wave 
functions [65]. 

A similar situation also exists for the case of hard photodisintegration involving nu- 
clei. Even though experiments clearly indicate the onset of s " scaling for the cross sec- 
tion of for example pn reactions at 90''c.m., one still expects sizable nonperturbative 
effects. Theoretical methods of calculation of these effects are very restricted. They use 
different approaches to incorporate nonperturbative contributions in the process of hard 
photodisintegration of the deuteron. The reduced nuclear amplitude (RNA) formalism 
includes some of the nonperturbative effects through the nucleon form factors[66, 67], 
while in the quark-gluon string model (QGS)[68] nonperturbative effects are accounted 
for through the reggeization of scattering amplitudes. Also recently, large cm. an- 
gle photodisintegration of the deuteron for photon energies up to 2 GeV was calcu- 
lated within point-form relativistic quantum mechanics approximation[69] in which the 
strength of the reaction was determined by short range properties of the nn interaction 
potential. 

In the QCD hard rescattering model (HRM)[70, 72] it is assumed that the energetic 
photon knocks-out a quark from one nucleon in the deuteron which subsequently ex- 
periences hard rescattering with a quark of the second nucleon. The latter leads to the 
production of two nucleons with large relative momentum. The summation of all the 
relevant rescattering diagrams results in a scattering amplitude in that the hard rescat- 
tering is determined by the large-momentum transfer pn scattering amplitude, which 
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includes noncalculable nonperturbative contributions. Experimental data are used to 
estimate the hard pn scattering amplitude. The HRM allows us to calculate the abso- 
lute cross section of 90° cm. hard photodisintegration of the deuteron without using 
additional adjustable parameters. 

Also, within the QGS[74] approximation and the HRM[43, 75] rather reasonable 
agreement has been obtained for polarization observables[60]. 

Hard photodisintegration of ^He 

Although all the above-mentioned models describe the major characteristics of hard 
photodisintegration of the deuteron they are based on very different approaches in the 
calculation of the nonperturbative parts of the photodisintegration reaction. To inves- 
tigate further the validity of these approaches it was suggested in Ref.[76] to extend 
the studies of high energy two-body photodisintegration to the case of large angle cm. 
breakup of two protons from the -^He target. In this case not only do the predictions 
of the above-described models (RNA, QGS, HRM) for absolute cross section diverge 
significantly, but also the two-proton breakup reaction from ^He provides additional ob- 
servables such as spectator-neutron momentum distributions that can be used to check 
further the validity of the models. 

Detailed analysis of reactions involving hard breakup of both pp and pn pairs from 
the ^He target is presented in chapter III. New insight into the nature of large cm. angle 
scattering is gained through the comparative study of pp and pn breakup processes. One 
important observation is that the relative strength of pp to pn breakup is larger than the 
one observed in low energy reactions. This characteristic is related to the onset of quark 
degrees of freedom in hard breakup reactions in which effectively more charges are 
exchanged between two protons than between proton and neutron. 

Another signature of the HRM is that the shapes of the energy dependencies of s"- 
scaled differential cross sections of pp and pn breakup reactions mirror the shapes of the 
energy dependencies of s"'-scaled differential cross sections of hard elastic pp and pn 
scatterings. 
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Within the HRM one observes also that pp and pn hard breakup processes are sensi- 
tive to different components of the -^He ground state wave function, resuking in different 
spectator-nucleon momentum dependencies for pp and pn hard breakup cross sections. 

Because of the different ground state wave function components involved in pp and 
pn breakup reactions, the HRM also predicts significantly different magnitudes for trans- 
ferred longitudinal polarizations for these two processes. 
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Chapter II 

PROTON NEUTRON ELASTIC SCATTERING 



This chapter looks into an asymmetry in the angular distribution of hard elastic proton- 
neutron scattering with respect to 90" center of mass scattering angle. It will be demon- 
strated that the magnitude of the angular asymmetry is related to the helicity-isospin 
symmetry of the quark wave function of the nucleon. An estimate of the asymme- 
try within the quark-interchange model of hard scattering demonstrates that the quark 
wave function of a nucleon assuming the exact SU(6) symmetry predicts an angular 
asymmetry opposite to that of experimental observations [110]. On the other hand the 
quark wave function derived from the diquark picture of the nucleon produces an asym- 
metry consistent with the data. Comparison with the data allows extracting the relative 
sign and the magnitude of the vector and scalar diquark components of the quark wave 
function of the nucleon. These two quantities are essential in constraining QCD models 
of a nucleon. Overall, it is concluded that the angular asymmetry of a hard elastic scat- 
tering of baryons provides a new venue in probing quark-gluon structure of baryons and 
should be considered as an important observable in constraining the theoretical models. 

For several decades elastic nucleon-nucleon scattering at high momentum transfer 
( t, -u > M% GeV^) has been one of the important testing grounds for QCD dynamics 
of the strong interaction between hadrons. Two major observables considered were the 
energy dependence of the elastic cross section and the polarization properties of the 
reaction. 

Predictions for energy dependence are determined by the underlying dynamics of the 
hard scattering of quark components of the nucleons. One such prediction is derived 
from the quark-counting rule [29, 30] according to which the differential cross section 
of two-body elastic scattering {ab cd) at high momentum transfer behaves like ^ ~ 
g-(na+ni,+nc+nrf)^ whcrc rcprcscnts the number of constituents in particle i (i=a,b,c,d). 

For elastic nn scattering, the quark-counting rule predicts s^^^ scaling which agrees 
reasonably well with experimental measurements (see e.g. Refs.[34, 33, 35, 36]). In 
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addition to energy dependence, the comparison [37] of the cross sections of hard exclu- 
sive scattering of hadrons containing quarks with the same flavor with the scattering of 
hadrons that share no common flavor of quarks demonstrated that the quark-interchange 
represents the dominant mechanism of hard elastic scattering for up to ISR energies (see 
discussion in [38]). 

For polarization observables, the major prediction of the QCD dynamics of hard 
elastic scattering is the conservation of helicities of interacting hadrons. The latter 
prediction stems from the fact that the gluon exchange in massless quark limit conserves 
the helicity of interacting quarks. 

The quark counting rule and helicity conservation however do not describe com- 
pletely the features of hard scattering data. The energy dependence of pp elastic cross 
section scaled by sj^V exhibits an oscillatory behavior which indicates the existence of 
other possibly nonperturbative mechanisms for the scattering[39, 40]. These expecta- 
tions are reinforced also by the observed large asymmetry, A„„ at some hard scattering 
kinematics [41] which indicates an anomalously large contribution from double helicity 
flip processes. These observed discrepancies however do not represent the dominant 
features of the data and overall one can conclude that the bulk of the hard elastic nn 
scattering amplitude is defined by the exchange mechanism of valence quarks which 
interact through the hard gluon exchange (see e.g. Refs.[42, 38]). Quark-interchange 
mechanism also reasonably well describes the 90 cm. hard break-up of two nucleons 
from the deuteron[70, 43]. 

However, the energy dependence of a hard scattering cross section, except for the 
verification of the dominance of the minimal-Fock component of the quark wave func- 
tion of nucleon, provides rather limited information about the symmetry properties of 
the valence quark component of the nucleon wave function. 

In this chapter it is demonstrated that an observable such as the asymmetry of a hard 
elastic proton-neutron scattering with respect to 90° cm. scattering may provide a new 
insight into the helicity-flavor symmetry of the quark wave function of the nucleon. 
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Namely we consider 



cr(6l) - cr(7r - 



(11. 1) 



where a{9) - is the differential cross section of the elastic pn scattering. We will dis- 
cuss this asymmetry in the hard kinematic regime in which the energy dependence of 
the cross section is ^ .s^°. Our working assumption is the dominance of the quark- 
interchange mechanism (QIM) in the nn elastic scattering at these kinematics. 




Figure II. 1 : Typical diagram for quark- interchange mechanism of nn nn 
scattering. 

Within the QIM, the characteristic diagram for pn elastic scattering can be repre- 
sented as it's shown in Fig.II.l. Here, one assumes a factorization of the soft part of 
the reaction in the form of the initial and final state wave functions of nucleons, and of 
the hard part which is characterized by the QIM scattering that proceeds with five hard 
gluon exchanges. This hard factor generates the energy dependence in accordance to 
the quark counting rule. In order to calculate the absolute cross section of the reaction, 
one needs to sum hundreds of diagrams similar to one of Fig.II.l. However, for the 
purpose of estimation of the asymmetry in Eq.(II.l) the important observation is that 
the hard scattering kernel is flavor-blind and conserves helicity. As a result, one expects 
that the angular asymmetry will be generated mainly through the underlying spin-flavor 
symmetry of the quark wave functions of the interacting nucleons. 

The amplitude of the hard elastic a + b ^ c + d scattering of Fig.II.l, within quark- 
interchange approximation, can be presented as follows: 



{cd\T\ab) = l«2,/3;,7;)(^I|a;,/32,72> 



x(a2,/32 , 72:ai/3W( I H \ Oi, A, 71, 02/3272) • (ai,/3i,7i I V'a) (a2, /32, 72 I V't), 
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(11.2) 



in which (at, a'.), (A,/?-) and (7^,7,') describe the spin-flavor quark states before and after 
the hard scattering, H, and 

Cl,,,^(a,A7 I (11.3) 

describes the probability amplitude of finding the a, ^,7 helicity-flavor combination of 
three valence quarks in the nucleon j[42]. 

To be able to calculate the c^^^ factors, one represents the nucleon wave function 
through the helicity-flavor basis of the valence quarks. We use a rather general form 
separating the wave function into two parts characterized by two (e.g. second and third) 
quarks being in spin zero - isosinglet and spin one - isotriplet states as follows: 



,3 — 1 l3 — -1 

^b(^.X^,^.,J-ir^^rf'^...)}. (11.4) 

1 -^.n^^ n,23 -'-'''23 2 ' iV 23 J 

in which j% and /ijv are the isospin component and the helicity of the nucleon. Here, the 
fci's are the light cone momenta of quarks. These momenta is represented by (xi,kn_), in 
which Xi is a light cone momentum fraction of the nucleon carried by the j-quark. We 
define Xj,h and tijs as helicity and isospin wave functions, where j is the spin, h is the 
helicity, / is the isospin and its third component. The Clebsch-Gordan coefficients are 
defined as (ji,TOi; j2,"^2 1 j,m). Here, represents the momentum dependent part of the 
wave function for (/ = 0, J = 0) and (/ = 1, J = 1) two-quark spectator states respectively. 
Since the asymmetry in Eq.(II.l) does not depend on the absolute normalization of 
the cross section, a more relevant quantity for us is the relative strength of these two 
momentum dependent wave functions. For our discussion we introduce a parameter, p: 

(^1.1) /TT c^ 

which characterizes an average relative magnitude of the wave function components 
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corresponding to (/ = 0, J = o) and (/ = i, J = i) quantum numbers of two-quark "specta- 
tor" states. Note that the two extreme values of p define two well know approximations: 
p = 1 corresponds to the exact SU(6) symmetric picture of the nucleon wave function 
and p = will correspond to the contribution of only the scalar diquark configuration 
in the nucleon wave function (see e.g. Ref.[44, 45, 46, 47] in which this component 
is referred as a scalar or good diquark configuration {[qq]) as opposed to a vector or 
bad diquark configuration denoted by {qq)). In further discussions p is kept as a free 
parameter. 

To calculate the scattering amplitude of Eq.(II.2), we assume the conservation of the 
helicities of the quarks participating in the hard scattering. This allows us to approxi- 
mate the hard scattering part of the amplitude, H, in the following form: 

Inserting this expression into Eq.(II.2) for the QIM amplitude, one obtains[42]: 

{cd\T\ ah) = Tri^r'^M"") (II.7) 

with: 

where we sum over the all possible values of /3 and 7. Furthermore, we separate the 
energy dependence from the scattering amplitude as follows: 

{cd\T\ ab) = (II.9) 

and define five independent angular parts of the helicity amplitudes as: 

cj,, = (++ I T{e) I ++); cf>, = {— I T{9) I ++>; 
03 - (+- I ne) I +-); 04 = -(-+ I T{e) I +-); 

05 = (-+ I T{e) I ++). (11.10) 



62 



Here, the "-" sign in the definition of 04 follows from the Jacob-Wick helicity convention[48] 
according to which a (-1) phase is introduced if two quarks that scatter to tt - e^^ angle 
have opposite helicity (see also Ref.[42]). 

Using Eqs. (11.7,11.8) for the non-vanishing helicity amplitudes of Eq.(II.lO) one ob- 
tains: 



for pp ^ pp'. 



{3 + y)F{9) + {3 + y)F{TT~0) 
{2-y)F{e) + il + 2y)Fi7r-e) 

-(1 + 2y)F(e) - (2 - y)F{n - 0) (II. 1 1) 



and for pn ^ pn: 



c^, = {2^y)F{e) + {l + 2y)F{^-e) 
= {2 + y)F{d) + {l + Ay)F{7r-9) 
(1)4 = 2yF{e) + 2yF{TT-e) (11-12) 

with 02 = 05 = due to helicity conservation. Here: 

y = x{x + 1) with X = — ^ (11.13) 

3(1 + 

and F{9) is the angular function. Note that the p = i case reproduces the SU(6) result of 
Refs.[42] and [38]. The results of Eqs.(II.l 1) and (11.12) could be obtained also through 
the formalism of the H-spin introduced in Ref.[38]. In this case, the helicity amplitudes 
are expressed through the average number of quarks to be found in a given helicity-spin 
state. These numbers will be directly defined through the wave function of Eq.(II.4). 

Introducing the symmetric and antisymmetric parts of the angular function F as 
follows: 

siO)^m±Iil^, aiO)^m^Ii^, (11.14) 
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Figure II.2: Asymmetry of pn elastic cross section. Solid dotted line - SU(6), 
with p = i, dashed line diquark-model with p = 0, solid line - fit with p = 0.3. 



and using Eq.(II.12) for the asymmetry as it is defined in Eq.(II.l) one obtains: 



6ai0)s{0){l ~ 2y ~ 3y^) 
a{0y{l - 3y)2 + 3s(6i)2(3 + 6y + 7y^) ' 



(11.15) 



One can make a rather general observation from Eq.(II. 15), that for the SU(6) model, 
(p = 1, y = |) and for any positive function, a{9) at 6* < f , the angular asymmetry has a 
negative sign opposite to the experimental asymmetry (Fig. II. 2). Note that one expects 
a positive a{e) at 6* < f from general grounds based on the expectation that in the hard 
scattering regime the number of t-channel quark scatterings dominates the number of 
It-channel quark scatterings in the forward direction. 

As it follows from Eq.(II.15), a positive asymmetry can be achieved only for i - 
2y - Sy'^ > 0, which according to Eq.(II.13) imposes the following restrictions on p: 
p < 0.49 or p > 2.036. The first condition indicates the preference for the scalar diquark- 
like configurations in the nucleon wave function, while the second one will indicate 
the strong dominance of the vector-diquark component which contradicts empirical 
observations [44, 45, 46]. 

In Fig.II.2 the asymmetry of pn scattering calculated with SU(6) (p = i) and pure 
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scalar-diquark (p o) models are compared with the data. In these estimates we use 
F{9) = C ■ sin-'^{6){i - cos{6))-^ as the dependence of the angular function[49], which is 
consistent with the picture of hard collinear QIM scattering of valence quarks with five 
gluon exchanges, and reproduces reasonably well the main characteristics of the angular 
dependencies of both pp and pn elastic scatterings. Note that using a form of the angular 
function based on nucleon form-factor arguments[38, 42], F « (i cos{9))-^ will result 
in the same angular asymmetry. 

The comparisons show that the nucleon wave function (II.4) with a scalar diquark 
component (p = o) produces the right sign for the angular asymmetry. On the other 
hand, even large errors of the data do not preclude to conclude that the exact SU(6) 
symmetry (p = i) of the quark wave function of nucleon is in qualitative disagreement 
with the experimental asymmetry. 

Using the above defined angular function F{9), Agoo is fitted in Eq.(II.15) to the data 
at -t, -u > 2 GeV^ varying p as a free parameter. We used the Maximal Likelihood 
method of fitting excluding those data points from the data set whose errors are too 
large for meaningful identification of the asymmetry. The best fit is found for 

-0.3 ±0.2. (11.16) 

The nonzero magnitude of p indicates the small but finite relative strength of a bad/vector 
diquark configuration in the nuclear wave function as compared to the scalar diquark 
component. It is intriguing that the obtained magnitude of p is consistent with the 10% 
probability of "bad" diquark configuration discussed in Ref.[46]. 

Another interesting property of Eq.(II.16) is the negative sign of the parameter p. 

Within a qualitative quantum-mechanical picture, the negative sign of p may indicate 
for example the existence of a repulsion in the quark- (vector- diquark) channel as op- 
posed to the attraction in the quark - (scalar-diquark) channel. It is rather surprising that 
both the magnitude and sign agree with the result of the phenomenological interaction 
derived in the one-gluon exchange quark model discussed in Ref.[45]. 

In conclusion, we demonstrated that the angular asymmetry of hard elastic pn scat- 
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taring can be used to probe the symmetry structure of the valence quark wave function 
of the nucleon. We demonstrated that the exact SU(6) symmetry does not reproduce 
the experimental angular asymmetry of hard elastic pn scattering. The use of nucleon 
wave functions consistent with the diquark structure results in an asymmetry in better 
agreement with the empirically observed. The fit to the data indicates 10% probability 
for the existence of bad/vector diquarks in the wave function of nucleons. It also shows 
that the vector and scalar qq components of the wave function may be in opposite phase. 
This will indicate a different dynamics for q - [qq] and q - (qq) interactions. 

The relative magnitude and the sign of the vector (qq) and scalar [qq] components can 
be used to constrain different QCD predictions which require the existence of diquark 
components in the nucleon wave function. These quantities in principle can be checked 
in Lattice calculations. The angular asymmetry studies can be extended also to include 
the scattering of other baryons such as A-isobars (which may have a larger fraction of 
vector diquark component), as well as strange baryons which will allow us to study the 
relative strength of (qq) and [qq] configurations involving strange quarks. 



66 



Chapter III 

HARD BREAKUP OF A NUCLEON NUCLEON SYSTEM IN THE He 

NUCLEUS 

This chapter investigates the large angle photodisintegration of two nucleons from the 
^He nucleus within the framework of the hard rescattering model (HRM). In the HRM 
a quark of one nucleon knocked out by an incoming photon rescatters with a quark of 
the other nucleon leading to the production of two nucleons with large relative momen- 
tum. Assuming the dominance of the quark-interchange mechanism in a hard nucleon- 
nucleon scattering, the HRM allows the expression of the amplitude of a two-nucleon 
break-up reaction through the convolution of photon-quark scattering, nn hard scat- 
tering amplitude and nuclear spectral function which can be calculated using a nonrel- 
ativistic -^He wave function. The photon-quark scattering amplitude can be explicitly 
calculated in the high energy regime, whereas for nn scattering one uses the fit of the 
available experimental data. The HRM predicts several specific features for the hard 
breakup reaction. First, the cross section will approximately scale as s Secondly, the 
s" weighted cross section will have the shape of energy dependence similar to that of s^" 
weighted nn elastic scattering cross section. Also one predicts an enhancement of the 
pp breakup relative to the pn breakup cross section as compared to the results from low 
energy kinematics. Another result is the prediction of different spectator momentum 
dependencies of pp and pn breakup cross sections. This is because of the fact that same- 
helicity pp-component is strongly suppressed in the ground state wave function of -^He. 
Because of this suppression the HRM predicts significantly different asymmetries for 
the cross section of polarization transfer nn breakup reactions for circularly polarized 
photons. For the pp breakup this asymmetry is predicted to be zero while for the pn it is 
close to |. 

This chapter is organized as follows: Section III.l, within the HRM, presents a 
detailed derivation of the differential cross section of the reaction of hard breakup of 
two-nucleons from a ^He target. In Section III. 2 the formulas derived in the previous 
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section are used to calculate the differential cross section of a proton-neutron breakup 
reaction, while in Section III.3 calculations are done for a two-proton breakup reaction. 
Section III.4 considers the relative contribution of two- and three-body processes for 
hard breakup reactions involving A>3 nuclei. In Section III. 5 numerical estimates are 
presented for differential cross sections of pn and pp breakup reactions. In Section III. 6 
the polarization transfer mechanism of the HRM is discussed, and estimates of the 
asymmetry of the cross section with respect to the helicity of the outgoing proton are 
presented. Results are summarized in Section III. 7. 

The details of the derivation of the hard rescattering amplitude are given in Ap- 
pendix C. The quark- interchange contribution to the hard nn elastic scattering am- 
plitude is discussed in Appendix D. Appendix B describes a method for calculating 
quark-charge factors within quark-interchange mechanism of nn hard elastic scatter- 
ing. A complete list of HRM helicity amplitudes for high energy two-nucleon breakup 
is presented in Appendix G for both, deuteron and ^'He photodisintegration. 

III.l Hard Photodisintegration of Two Nucleons from ^He 
Reference frame and kinematics 

We are considering a hard photodisintegration of two nucleons from the -^He target 
through the reaction: 

7 He ^ (NN) + N, , (III. 1 ) 

in which two nucleons (NN) are produced at large angles in the "^-nn" center of mass 
reference frame with momenta comparable to the momentum of the initial photon, 
q i>i GeV/c). The third nucleon, N,, is produced with very small momentum p, < jtiat. 
(Definitions of four-momenta involved in the reaction are given in Fig.III.l.) 

We consider "7-7V7V" in a "g+ = o" reference frame, where the light-cone momenta' 

'The light-cone four- momenta are defined as {p+,p-,p±), where p± = E zLpz- Here the z axis is defined 
in the direction opposite to the incoming photon momentum. 
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Figure III.l: Typical hard rescattering diagram for the nn photodisintegra- 
tion from the -'He target. 

of the photon and the nn pair are defined as follows: 



= {q+,q-,q±) = (o, V^'nn, o), 



M' 



Pnn = {Pnn+:Pnn~,Pnn± 

) = (7^, 0), (III.2) 

where p'^.j^ = p'^-^^ - p^, Ml^ = p^;^,PJv^^,^, and s'^^ = s^n - Ml^. Here the invariants, s^^m 
and defined as follows: 

SJVAT = (g + PnnT = [Pfl + Pf2f 

tNN = {q-pflf = {Pf2-PNNf- (III. 3) 

As it follows from Eq.(III.2) in the limit of the = o" reference frame 

N 

coincides with the center of mass frame of the 7-iViV system. As such it is maximally 
close to the reference frame used for the 7^ ^ p + n reaction in Refs.[70] and [43]. 
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Hard rescattering model 

The hard rescattering model is based on the assumption that in the hard two-nucleon 
photodisintegration reaction, two nucleons with large relative momenta are produced 
because the hard rescattering of a fast quark from one nucleon with a quark from the 
other nucleon. In this scenario the fast quark is knocked out from a low-momentum 
nucleon in the nucleus by an incoming photon. This approach is an alternative to the 
models in which it is assumed that the incoming photon breaks the preexisting two- 
nucleon state which has very large relative momentum in the nucleus. 

The validity of the HRM is derived from the observation that the ground state wave 
functions of light nuclei peak strongly at small momenta of bound nucleons, p^O. Thus, 
diagrams in which an energetic photon interacts with bound nucleons of small momenta 
will strongly dominate the diagrams in which the photon interacts with bound nucleons 
that have relative momenta p>m,^. 

The resulting scenario that the HRM sketches out is as follows (see e.g. Fig.III.l): 
first, the incoming photon will knock out a quark from one of the nucleons in the nucleus 
and then the struck quark that now carries almost the whole momentum of the photon 
will share its momentum with a quark from the other nucleon through the exchanged 
gluon. The resulting two energetic quarks will recombine with the residual quark-gluon 
systems to produce two nucleons with large relative momentum (^q). This recombina- 
tion will contain gluon exchanges and also incalculable nonperturbative interactions. 

Note that for the quark-gluon picture discussed above to be relevant the interme- 
diate masses produced after the photon absorption should exceed the mass scale 
characteristic for deep inelastic scattering, w ^ 2.2 GeV. Using the relation mi„i « 
i/mlr + 2E^mM, from the requirement that > W one obtains the condition E^>2 GeV. 
Additionally, to ensure the validity of quark degrees of freedom in the final state rescat- 
tering, one requires k^^i > 1 GeV/c for the relative momentum fcrei, of two outgoing nucle- 
ons. All of these imposes a restriction on the incoming photon energy, E.,>2 GeV, and 
for transferred momenta -t, -u>2 GeV^ Note that provided a smooth transition from 
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hadronic to quark- gluon degrees of freedom in nuclei one expects the validity of the 
HRM to extend to even lower values of (>i GeV). This expectation was confirmed 
in recent measurements of angular dependencies of the -/d pn cross section for a wide 
range of incoming photon energies[59]. 

To calculate the differential cross section of the hard photodisintegration reaction of 
Eq.(ni.l) within the HRM one needs to evaluate the sum of hard rescattering diagrams 
similar to the one presented in Fig.III.l. We start with analyzing the scattering ampli- 
tude corresponding to the diagrams of Fig.III.l. Using Feynman rules and applying the 
light-cone wave function reduction described in Appendix C, we obtain 



X 



2 



dxi d^kijL dx2 d^fcax *3He "{a,P±,Ps)da d'^p^ nu A\ 

^ ' ^""^^ 2(2^^ 2(2^0^ Jl^) ~2(2^"^^^^^^^^^' ^^^^-^^ 

where the (pi/ < — >p2f) part accounts for the diagram in Fig.III.l(b). Here the four- 
momenta, pu, P2i, Ps, fci, ^2, r, pif and P2f are defined in Fig.III. 1 . Note that ki and fcz define 
the four-momenta of residual quark-gluon system of the nucleons without specifying 
their actual composition. We also define xi, x[, X2 and x'^ as the light-cone momentum 
fractions of initial and final nucleons carried by their respective residual quark-gluon 
systems: xu2) = hB>±. and .t'„, = . For the 'He wave function, a = is the 

•' ^ ' Pl(2)i+ '-'■^> Pl(2)/+ PNN + 

light-cone momentum fraction of the nn pair carried by one of the nucleons in the pair, 
and is their relative transverse momentum. The scattering process in Eq.(III.4) can 
be described through the combination of the following blocks: (a) *3He^"'^'''^"(a,px,Ps), 
is the light-cone ^He-wave function that describes a transition of the -'He nucleus with 
helicity A 4 into three nucleons with A^ , Aa^, and A, helicities, respectively, (b) The 
term in describes the "knocking out" of an Tyn-helicity quark from a An-helicity 
nucleon by an incoming photon with helicity A^. Subsequently, the knocked-out quark 
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exchanges a gluon, ([-igTfY]), with a quark from the second nucleon producing a final 
772/-helicity quark that combines into the nucleon "2/" with helicity A2/. (c) The term in 
{...}2 describes the emerging T^az-helicity quark from the Aa^ -helicity nucleon then then 
exchanges a gluon, ([-igTfi"]), with the knocked-out quark and produces a final 771/- 
helicity quark that combines into the nucleon "1/" with helicity Ai/. d) The propagator 
of the exchanged gluon is C^ir) = with polarization matrix, d"" (fixed by the light- 
cone gauge), and r = (^2 - ^2 + - iPi - ki + q), with l = {p^f - P2,). In Eq.(III.4) the V^'" 
represents everywhere an r/-helicity single quark wave function of a A-helicity nucleon 
as defined in Eq.(C.14) and is the quark spinor defined in the helicity basis. 

The denominator of the struck quark's propagator can be represented as follows: 

(pii - ki + qY - TO^ + = (1 - Xi)s'^^(ac - a + ie), (III. 5) 

where 



etc = 1 



ml{l - .Ti) + rriixi + (fci - XiPi)^ 



(III.6) 



xi(i-xi) 

Here ml and m]^ « m|r are defined in Eqs.(C.8) and (C.l 1), and m, represents the current 
quark mass of the knocked out quark. In what follows, we use the fact that the ^He wave 
function strongly peaks at a = i, which corresponds to the kinematic situation in which 
two constituent nucleons have equal share of the nn pair's light-cone momentum. Thus 
one expects that the integral in Eq.(III.4) is dominated by the value of the integrand at 
\. This allows us to perform a-integration in Eq. (111.4) through the pole of the 
denominator (III. 5) at a = a^. I.e. keeping only the i7r(5( ) part of the relation 

1 . 1 

-i'K5{a — ckc) + P , 



ac — a + IE 



and later replacing by \. Using this relation to estimate the propagator of the struck 
quark at its on-mass shell value (a = a J allows to write, 

-Hi + - E "c(^'i - ^1 + 9)^c(Pi -k,+q) 
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Then for the scattering amplitude of Eq.(III.4) one obtains 



(Ai/,A2/,As I Ai I Xj,Xa) 



O^^Iy HQ.ee,^7 ]u,Jp^.-h) 1^ 

\ L — Xi [L — X2) 



^ w — w " ^"'^ ^ '^'^ ■ ^ ^ ^ 

Next, we evaluate the matrix element of the photon-quark interaction using on-mass 
shell spinors for the struck quark. Taking into account the fact that {pu - k)+ > Ifcx^m,, 
for this matrix element we obtain 

uc{pu - fci + q)[~iQ.eei^i^]u,,^^(j)u - k-,) = ieQa^2E^^{X^)5^-<^5^-'^^^ (III.8) 

where E, = {i- a)(i - xO^^ and = (i - (i - a)(i - x^))^^. 



Further explicit calculations of Eq.(III.7) require the knowledge of quark wave func- 
tions of the nucleon. Also, one needs to sum over the multitude of the amplitudes rep- 
resenting different topologies of quark knock-out rescattering and recombinations into 
two final nucleon states. 

This difficulty can be circumvented by again using that the -^He wave function strongly 
peaks at a = i . We evaluate Eq.(III.7) setting everywhere = I- Such approximation 
significantly simplifies further derivations. As it follows from Eq.(III.6) the a^ = ^ con- 
dition restricts the values of xi of the recoil quark-gluon system to xi 4^, thereby 
ensuring that the quark-interchange happens for the valence quarks with x^ = i - x^ ^ i. 
The latter allows us to simplify Eq.(III.8) setting Ej = E2 ^ Using these approxi- 

mations and substituting Eq.(III.8) into Eq.(III.7) one obtains 



(Ai/,A2/A<, I M, I A^,Aa) = i[X^]e ^ / 



X 



(Vlf,V2f),{V2i),{>^li,^2i) ' 



12^ 



73 



1 - x\ 



{Plf 



'2i 



dxi cPki± dx2 (Pk2± 



me 




(2^)2 



{p,f^P2f). (III.9) 



Xi 2(27r)3 X2 2(27r)3 



JQIM 



Note that due to the 5 factors in Eq.(ni.8) the helicity of the knocked out quark in 
Eq.(III.9) is equal to the helicity of incoming photon, that is r]u = \^,. 

To proceed, we observe that the kernel, [. . .]qim representing the quark-interchange 
mechanism (QIM) of the rescattering in Eq.(ni.9) can be identified with the quark- 
interchange contribution in the nn scattering amplitude (see Appendix D). Such iden- 
tification can be done by observing that in the chosen reference frame, q+ = 0, and 
the quark wave function of the nucleon depends on the quark's light-cone momen- 
tum fraction and transverse momentum only, which are the same in both Eqs.(ni.9) 
and (D.4). For our derivation we also use the above-discussed observation that the 
a = I condition ensures that the quark-interchange happens for the valence quarks 
with .T, = 1 xi - 1. This justifies our next assumption, that valence quarks carry the 
helicity of their parent nucleon (i.e. 7]u = A,). The last assumption allows us to perform 
the summation of Eq. (III. 9) over the helicities of the exchanged quarks ('n2i,vif,V2f) and 
to use Eq.(D.4) to express the QIM part in Eq.(III.9) through the corresponding QIM 
amplitude of a^a^ scattering. Summing for all possible topologies of quark-interchange 
diagrams we arrive at 



where nucleon momenta pi and P2 have half of their cm. momentum fractions and is 
their relative transverse momentum with respect to the direction of the photon momen- 
tum [see Eq.(III.16)]. Here, for example, Qf • (Aa/; | T^™is,P) \ Ai;A2) represents the 



(Ai/,A2/, A, I M I X^,Xa) = ie[X^] x 



+ 
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quark-interchange amplitude of nn interaction weighted with the charge of those inter- 
changing quarks Qf that are struck from a nucleon N by the incoming photon. The sum 
(J2) can be performed within the quark- interchange model of nn interaction which 
allows us to represent the nn scattering amplitude as follows [38]: 

{a'b'\T2T\ab)^^{a'b'\ [U,+rvr,]F,As,t)\ab) (III.ll) 

where /, and are the identity and Pauli matrices defined in the SU(2) flavor (isospin) 
space of the interchanged quarks. The kernel Fij{s,t) describes an interchange of i and 
j quarks. 

Using Eq.(III.ll) one can calculate the quark-charge weighted QIM amplitude, Qi 

{a'b'\T^^^^\ab), as follows: 

^Qf(a'6'|rQ™|a6) = ^{a'b'\ ^ [IJ, + f^r,m)F,A^,t)\ab) 

= Q^.(a'&'|TQjr|afe), (III.12) 

where are the charge factor that are explicitly calculated using the method described 
in Appendix C. These factors can be expressed through the combinations of valence 
quark charges Qi of nucleon N and the number of quark interchanges for each flavor of 
quark, Nq^, necessary to produce a given helicity nn amplitude, as follows. 



N^u + Nad + N^d 

Next we discuss the light-cone wave function of ^He that enters in Eq.(III.lO). The 
important result that allows us to evaluate the wave function is the observation that two 
nucleons that interact with the photon share equally the nn pair's cm. momentum 
(pnn), i.e., a = 5. If we constrain the third nucleon's light-cone momentum fraction 
a, = « 1 and transverse momentum < m,^, then the momenta of all 

the nucleons in the nucleus are nonrelativistic. In this case one can use the calculation of 



'The additional assumption of helicity conservation allows us to express the kernel in the form[38] 
Fi j {s,t) = i [lilj + (Ti^aj]Fij (s, t), where li and operate in the SU(2) helicity (iJ-spin) space of exchanged 
{i, j) quarks[38]. However for our discussion the assumption of helicity conservation is not required. 
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triangle diagrams, which provides the normalization of nuclear wave functions based on 
baryonic number conservation to relate LC and nonrelativistic nuclear wave functions 
as follows[77, 78] 



(III. 14) 

where for *3He,NR we can use known nonrelativistic ^He wave functions (see e.g., [80]). 

Substituting Eqs.(III.12) and (III.14) into Eq.(III.lO) for the two-nucleon photodis- 
integration amplitude we obtain 



1/ I ^JVJV (SiVJVj^Jv) I -^71 •^2i)^3He.NR(^'l' '^T'-P2: '^aiSPsj As)..„jv /„ N 



(III. 15) 



where in the Lab frame of the ^He nucleus, defining the z direction along the direction 
of QLab one obtains 

^ g2-P2. . „ ^ ~ -Pax 
"» = 17I7T; Pi+P2 = -Ps, (III.16) 

with all the momenta defined in the Lab frame. 

Equation(III.15) allows us to calculate the unpolarized differential cross section of 
two nucleon breakup in the form 

da IMP ^jjj ^rj^ 



dtd'pJ{2E,{2nY) 167r(.s - Ml){s^t, - 

where s = (fc^ +pa)^ and 



|Mp = i.i \{Xif,X2f,Xs\M\X,,X^)\\ (10.18) 

Al/,A2/,As,A-j,,A^ 
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As follows from Eq.(ni.l5) the knowledge of quark- interchange helicity amplitudes 
of NN elastic scattering will allow us to calculate the differential cross section of hard 
NN breakup reaction without introducing any adjustable parameter. 



Because the assumption of a, = 1 plays a major role in the above derivations we 
attempt now to estimate the theoretical error introduced by this approximation. This 
approximation by its nature is a "peaking" approximation that is used in loop calcu- 
lations involving Feynman diagrams (one such example is the calculation of radiative 
effects in electroproduction processes; see, e.g., [79]). One way to estimate the accuracy 
of the approximation is to identify the main dependence of the integrand in Eq.(III.7) 
on which can be evaluated exactly and compare with its evaluation at = 5. Using 
Eq.(III.8) as well as Eq.(III.6) that allows us to relate ^ to and assuming that the 
quark wave functions of nucleons at ^ are less sensitive to a^, one arrives at 



This ratio depends on the kinematics of the spectator nucleon, and for the case of < 
100 MeV/c, R{ps) w 1.1, which corresponds to --20% of uncertainty in the cross section 
of the reaction calculated with the q^c = 5 approximation. The uncertainty increases 
with an increase of the momentum of the spectator nucleon. This can be understood 
qualitatively because, for large center of mass momenta of the nn pair, the a = ^ peak 
of the nuclear wave function is less pronounced. 

Quark-interchange and hard A^A^ elastic scattering amplitudes 

The possibility of using NN elastic scattering data to calculate the cross section in 
Eqs.(III.17) and (III. 18) is derived from the assumption that the quark-interchange 
mechanism provides the bulk of the NN elastic scattering strength at high energies and 
large cm. angles. This is a rather well-justified assumption. Experiments on exclu- 
sive large -t two-body reactions [82] demonstrated clearly the dominance of the quark- 
interchange mechanism for the scattering of hadrons that share common quark flavors. 



R{ps) 



3He 



(III. 19) 



/ 
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The analysis of these experiments indicate that contributions from competing mecha- 
nisms such as pure gluon exchange or quark- antiquark annihilation are on the level of 
few percent. This fact justifies our next approximation, to substitute quark- interchange 
A'^A'^ amplitudes in Eq.(III.15) with actual nn helicity amplitudes as follows: 



< 4 




,+ > = 




< -i 




, - > = 




< -i 


I^QIM, 

' ' Kjvat I 


, - > = 


^2 


< 4 


I^QIM, 


> = 


4>3 


< -i 


KATJV 1 


- + > = 


-04 



(III.20) 



All other helicity combinations can be related to the above amplitudes through the parity 
and time -reversal symmetry. The minus sign in the last equation above is due to the 
Jackob-Wick phase factor (see, e.g., Ref.[42]), according to which one gains a phase 
factor of ( i) if two quarks that scatter by n-e angle in cm. have opposite helicities [48] . 
Note that ^/s are normalized in such a way that the cross section for nn scattering is 
defined as 



da' 



dt 



1677 S{S - 4771%) 2 



+ 4105^. 



(III.21) 



Because in the hard breakup regime the momentum transfer -t^ > m^, one can 
factorize the helicity nn amplitudes from Eq.(III.15) at s^n and values defined as 
follows: 



Snn 

t]V 



(q+PNNf = (Pfl +P/2)^ 

(P/2 -PiViv/2) = — + — —■ 



(III.22) 



Using this factorization in Eq.(III.15) for the spin averaged square of the breakup am- 
plitude one obtains 



|Mf = ^^^^ J {2Ql\<l>,\'So + Qlm' + \H')Si2- 
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[{Q^'h + QF'<p4r + (Q~^04 + Q~^03)')] ^34} , 



(III.23) 



where Qp = Qp' + Q%'' and ^iz, 534, and So are partially integrated nuclear spectral func- 
tions: 



Si2{ti,t2,a,p^) = Nnn ^ 



3He,NR 



(27r)2 



, (III.24) 



At— — Aq — — -H- A;^ — — T 



3He,NR 



(pi , Ai , 1 1 ; p2 , A2 , ^2 ; , A3)miv 



(27r)2 



(III.25) 



and 



'S'o — S12 + iSa. 



(III.26) 



In the above equations ti and t2 are the isospin projections of nucleons in ^He. The wave 
function is normalized to | for proton and 1 for neutron. The normalization constants, 
A^jvAT renormalize the wave function to one pp and two np effective pairs in the wave 
function with Npp = i and Np„ ^ 4. 

Equations.(III.17) and (III.23) together with Eqs.(III.20),(III.24),and (III.25) allow 
us to calculate the differential cross section of both pp and pn breakup reactions off the 
^He target. Notice that, on the qualitative level, as it follows from Eqs.(III.17) and 
(III.23) in the limit of s > Mi^^ and snn > m%, the HRM predicts an s " invariant 
energy dependence of the differential cross section provided that the NN cross section 
scales as s However the numerical calculations of Eq. (III. 23) require a knowledge of 
the NN helicity amplitudes at high energy and momentum transfers. Our strategy is to 
use Eq.(III.21) to express nn breakup reactions directly through the differential cross 
section of pn and pp elastic scatterings rather than to use helicity amplitudes explicitly. 

III.2 Hard breakup of proton and neutron from ^He. 

We consider now the reaction 



He ^ (pn) 



(III.27) 
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in which one proton is very energetic and produced at large cm. angles with the neu- 
tron, while the second proton emerges with low momentum < lOO MeV/c In this case 
the hard rescattering happens in the pn channel. Using the « ^4 relation for hard pn 
scattering amplitude (see, e.g., Refs [42, 38, 49]) for breakup amplitude of Eq.(III.23) 
one obtains 

IMf = 1 ^2|<^^|25„ + + + (|<^3|2 + \mSs,} , (III.28) 

where Q'-p^ ^Q'^p + Qp can be calculated using Eqs.(III.12) and (III.13). On the basis of 
the SU(6) flavor-spin symmetry of nucleon wave functions, for the helicity amplitudes 
of Eq.(III.20) using the method described in Appendix B one obtains 

Q7 = I- (111.29) 

We can further simplify Eq.(III.28) noticing that for the pn pair in ^He one has S'12 ~ 
5*34 w ^. This is due to the fact that in the dominant S state two protons have opposite 
spins and therefore the probability of finding one proton with a helicity opposite to that 
of the neutron is equal to the other proton having the same helicity as the neutron's. 
Using this relation and Eq.(III.21) for the pn breakup reaction one obtains 

I M I = 167rsjvjv(siviv - 4tojv) — —■ (lll.JUj 

Inserting it in Eq.(III.17) for the differential cross section one obtains 



dt^ ^'•"^ 2 isr,N-plr,)'i-^-Mi^J dt 



(III.31) 



where a = ^ and is the differential cross section of hard pn scattering evaluated 

at values of snn and defined in Eq.(III.22). The spectral function Sl'^ is defined in 
Eq.(III.26) and corresponds to: 



1 

5r(«,p.) = 4 



/ *3He,NR(^i''^i' 2':P-2,>^2,--j^]PsA'j.)mi 
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(III.32) 



III.3 Hard breakup of two protons from ^He 

We now consider the reaction 

j+'llQ->{pp)+n, (III.33) 

in wliich two protons are produced at large cm. angles while the neutron emerges as a 
spectator with small momentum (p, < lOO MeV/c). 

The relation between and 5*34 is very different from that in the pn case. As a 
result of the fact that two protons cannot have the same helicity, in the S state one has 
that 5i2 < 534- The estimates of the spectral functions based on the realistic ^He wave 
function[80] gives |^ ^ 10 *. Therefore one can neglect the S12 term in Eq.(III.23). 
The next observation is that for pp scattering the helicity amplitudes (^3 and ^4 have 
opposite signs because of the Pauli principle (see, e.g., Refs.[48, 42]). Using the above 
observations and neglecting the helicity-nonconserving amplitude 05 for the pp-breakup 
amplitude we obtain 

= ^^^^i {2(Q-,|03| " Ql\cb,\rs,,} . (III.34) 

The charge factor depends on the helicity amplitude it couples; therefore one esti- 
mates it for the combination of (Q^l^sl - QfI^Ija])- Using SU(6) symmetry for the distri- 
bution of given helicity-flavor valence quarks in the proton and through the approach 
described in Appendix B we obtain 

(q'fI^sI - Q'^rn) = Q'/m - m) (111.35) 

with 

QW = I (III.36) 

It is worth noticing that because of explicit consideration of quark degrees of freedom 
the effective charge involved in the breakup is larger for the case of two protons than 
for proton and neutron [see Eq.(III.29)]. This is characteristic of the HRM model in 
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which a photon couples to a quark and more charges are exchanged in the pp case than 
in the pn case. This is rather opposite to the scattering picture considered based on 
hadronic degrees of freedom in which case the photon will couple to an exchanged 
meson and pp contribution will be significantly suppressed because no charged mesons 
can be exchanged within the pp pair. 

To be able to estimate the cross section of the pp breakup reaction through the elastic 
pp scattering cross section we introduce a parameter 



^2 ^ ^ _ ^ 



(III.37) 



which allows to express the differential cross section of the reaction (III.33) in the fol- 
lowing form: 



dt 



(III.38) 



where we also introduced a factor /3, 



l'^3|-|04| 



(III.39) 



which accounts for the suppression from the cancellation between ^3 and ^4 helicity 
amplitudes of elastic pp scattering.' The spectral function SH in Eq.(III.38) is expressed 
through the ^He wave function according to Eq.(III.25) as follows: 



Sllia.Ps) 



E E 



(III.40) 



III.4 Two- and three-body processes in nn breakup reactions 



For a two-body hard nn breakup mechanism to be observed it must dominate the three- 
body /two-step processes. This is especially important for breakup processes, (III.33) 



'This cancellation was overlooked in earlier estimates of the cross section of pp photodisintegration from ^He 
target (see, e.g., [76]). 
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Ps p n 

Figure III. 2: Diagram corresponding to three-body processes in which the 
hard breakup of the pn pair is followed by a soft charge-exchange rescattering 
of the neutron off the spectator proton. 

because according to Eqs.(III.38 and III.39) the two-body contribution is suppressed 
because of a cancellation between (/)3 and (j)4, helicity amplitudes. 

At low to intermediate range energies (e^ ^ 200 MeV) it is rather well established 
that the pp breakup reaction proceeds overwhelmingly through a two-step (three-body) 
process[83, 84, 85, 86, 87, 88] in which the initial breakup of the pn pair (dominated by 
7r± exchange) is followed by a charge-interchange final state interaction of the neutron 
with the spectator proton. Other two-step processes include the excitation of intermedi- 
ate A isobars in the pn system with the subsequent rescattering off the spectator neutron, 
which produces two final protons. 

The dominance of three-body processes at low energies is related mainly to the fact 
that the two-body pp breakup is negligible because of the impossibility of charged-pion 
exchanges between two protons that absorb an incoming photon. 

At high energy kinematics within the HRM the interaction between protons is car- 
ried out by exchanged quarks because of which the relative strength of pp breakup is 
larger. 

To estimate the strength of three-body processes at high energy kinematics, one 
needs to calculate the contribution of diagrams similar to Fig. III. 2. Because the charge- 
exchange rescattering at the final stage of the process in Fig.III.2 takes place at proton 
momenta p'^. > m^, one can apply an eikonal approximation[91, 92] to estimate its con- 
tribution. 
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For > 2 GeV assuming that the HRM is valid for the first {pn breakup) stage 
of the reaction, for the ampfitude of three-body /two- step process within the eikonal 
approximation[91, 92] one obtains 

1\/r ^Qp.pni'^'Try rphard /-i \ „ 

/ ^3HeNR(Pi'^i'P2'^2;P'' ~ ^_P±^_h±^ (III.41) 

where we suppressed heficity indices for simpHcity and choose the isospin projections, 
ti = -t2 ^ \, corresponding to the initial pn pair that interacts with the photon. Here 
^p^^pn(^Af) is the hard elastic pn scattering amplitude and T^^^:^^^ represents the amplitude 
of the soft charge-exchange pn scattering. Because of the pion-exchange nature of the 
latter it is rather well established that this amplitude is real and can be represented as 
^^AeA\ where A and b are approximately constants [89]. 

Two main observations follow from Eq.(III.41) and the above-mentioned property 
of the charge-exchange amplitude: First, three-body and two-body amplitudes [see, 
e.g., Eq.(III.15)] will not interfere, since one is real and the other is imaginary. The 
fact that these two amplitudes differ by order of i follows from the general structure of 
rescattering amplitudes (see, e.g., Ref.[92]). Equation (III. 15) corresponds to a single 
rescattering amplitude, while Eq.(III.41) to a double rescattering amplitude. 
Second, because of the energy dependence of the charge-exchange scattering amplitude 
at small angles, the three-body contribution will scale like s^^^ as compared to the two- 
body breakup contribution. 

Using Eq.(III.41) one can estimate the magnitude of the contribution of three-body 
processes in the pp breakup cross section as follows: 

J 7^He^(pp)n , 'y^He^(pn)p qpnp ( \ 

"''-'three-body "'"two-body KPa J fill A1\ 

dt^ ^ dt^ SE"{Ps) ^ ^ ^ 

where S'o"(Po) is defined in Eq.(III.32) and for S^^^ips) based on Eq.(III.41) one obtains 
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Here both spectral functions are defined at a = i. 

Using Eqs. (111.42) and (III.43) and the parametrization of T^ll''_i,^^ from Ref.[89] one 
can estimate the relative contribution of three-body processes numerically. Note that 
this contribution is maximal at = i and increases with an increase of the momentum 
of p,. However because of the charge-exchange nature of the second rescattering, this 
contribution decreases linearly with an increase of s.' 

III.5 Numerical Estimates 

For numerical estimates we consider the center of mass reference frame of ^nn system, 
for which according to Eq.(III.3) one obtains 



tN,N = ~ 



{snn _ _ 4m^cos(^c.m.)) + <, (III.44) 



where m^^ = p^j^pnn,,, and pj^jv = Psbs " ^'^ Using the above equation we obtain for 
tjv [Eq.(III.22)], which defines the effective momentum transfer in the NN scattering 
amplitude, the following relation: 

^ _ (.jvjv-Ag^) ^ _ ^,^^_4^^,„,(0^^^)) + „ Mhi. (III.45) 

One can also calculate the effective cm. angle that enters in the NN scattering amplitude 
as follows: 



005(6*^ ) ^ 1 (^^^ ~ ^nn) (y/g^ - V Snn - 4m^cos(6>c.m.)) ^ "^m^ - M^^ (HI 46) 

2(sjvjv - 4m|r) y/SNN 2(sjvjv - 4:171%) ' 

The above equations define the kinematics of hard nn rescattering. 

Energy dependence and the magnitude of the cross sections 

We are interested in energy dependences of the hard breakup reactions of Eqs.(III.27) 
and (III.33) at fixed and large angle production of two fast nucleons in the "-f-NN'' cen- 
ter of mass reference frame. Particularly interesting is the case of 6'c.m.=90° for which 



'Notice that for the case of diagonal pn — > pn rescattering Tpn^np{k±) = scrtote 2 ' 
probabiHty of the rescattering is energy independent. 
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Figure III. 3: . Energy dependence of s" weighted differential cross sections 
at 90° cm. angle scattering in "7-iViV"system. In these calculations one inte- 
grated over the spectator nucleon momenta in the range of O-iOO MeV/c. 

as it follows from Eq.(III.46) cos{e^,-„ ) = 0.5. This means that the cross sections of hard 
breakup reactions at these kinematics will be defined by the nn elastic scattering at 
^c.m. = 60°. In Fig.III.3 the E., and s dependencies of the s^^V weighted differential cross 
sections are presented for the cases of the pp and p7i breakup reactions. In the calculation 
we integrated over the spectator nucleon's momentum in the range of (0-100) MeV/c 
and over the whole range of its solid angle. Also for the parameter C in Eq.(III.37) we 
used C = i, consistent with an estimate obtained within the quark-interchange model of 
pp scattering (see, e.g., [42, 38]). The estimation of the factor /3, which takes into account 
the cancellation between ^3 and 04 helicity amplitudes in Eq.(III.38) requires the knowl- 
edge of the angular dependence for helicity amplitudes. For this we used the helicity 
amplitudes calculated within thequark-interchange model[38, 42] with phenomenolog- 
ical angular dependencies estimated using F(6'c.m.) - s.nHe.r.)(i-cos(9^r.)y' (see, e.g., Refs. 
[38, 49]) which describes reasonably well the data at hard scattering kinematics. 

Several features of the HRM calculations are worth discussing in Fig.III.3: First, 
the breakup cross sections in average scale like s^^^. Note that the absolute (nonscaled) 
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Figure III.4: Invariant energy dependence of s^" weiglited differential cross 
sections of elastic vv scattering at 6'c.m. = 90° and e^.ra. = 60°. are fitted using the 
parametrization described in Appendix F of the available world data [90]. 



values of the cross sections drop by five orders of magnitude in the 2-8 GeV of photon 
energy range. Next, the shapes of the weighted differential cross sections reflect the 
shapes of the s^" weighted differential cross sections of w andpn scattering at 6'c.m. = 60°. 
[see Figs.(III.4) and (III. 5)]. It is worth noting that as follows from Figs. (III. 4) and 
(III. 5) the fits used in the calculation of pp and -pn breakup reactions contain uncertain- 
ties on the level of 10% for -p-p breakup (for sjvjv > 24 GeV^) and up to 30% for pn breakup 
reactions. Consequently, one can conclude that the calculated shape of the energy de- 
pendence of the pn breakup reaction in Fig.(III.3) does not have much predictive power. 
However, for the pp breakup the calculated shape, for up to satat < 24 GeV^ is not ob- 
scured by the uncertainty of the pp data and can be considered as a prediction of the 
HRM. 

Analysis of the first experimental data on pp photodisintegration of the ^He nucleus 
at high momentum transfer from Jefferson Lab [94] has shown excellent agreement 
with the HRM predictions. At d^,^ = 90° and s >12GeV^ the cross section scales with 
s^^ according to counting rules, and as shown in Figs. III. 6 and III. 7, the experimental 
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energy distribution of the breakup process at this kinematics is remarkably similar but 
for a power of s, to the experimental energy distribution of pp elastic scattering at = 
60°. 

It is worth mentioning that considered features of the HRM are insensitive to the 
choice of the above-discussed parameters of c and p, because they only define the ab- 
solute magnitude of the pp breakup cross section. 

pn pn 
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Figure III. 5: Invariant energy dependence of s"' weighted differential cross 
sections of elastic pn scattering at Ocm. = 90° and 6'c.m. = 60". Curves are fitted 
using the parametrization described in Appendix F of the available world 
data [90]. 

The next feature of the calculations in Fig. III. 3 is the magnitude of the pn and pp 
breakup cross sections. The pn breakup cross section [Eq.(III.31)] does not contain any 
free parameter, and similar to the HRM prediction for the breakup of the deuteron[70], 
it is expressed through the rather well-defined quantities. For the estimate of pp breakup, 
however, one needs to know the relative strength of the 03 and 04 amplitudes as com- 
pared to 01 as well as the extent of their cancellation at kinematics of s^n and defined 
in Eqs.(III.22) and (III. 45). Our calculation, determined by phenomenologically justi- 
fied estimates of factors C and I3 in Eq.(III.38) results in the pp breakup cross section 
which is about ten times smaller than the cross section for the pn breakup. This result 
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Figure III.6: Cross section data for pp breakup in ^He photodisintegration at 
90° cm. angle of the j+pp system [94]. The ^ lle{'j,pp)n events were selected 
with p„ < 100 MeV. (Dashed) QGS prediction times 1/5. (Long dashed) RNA 
prediction times 1/200. (Shaded) HRM prediction. 
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Figure III.7: pp elastic scattering data at 60° of center of mass angle of scat- 
tering. These data is used in Eq.(III.38) to calculate the HRM cross section 
of pp breakup shown in Fig. III. 6. 
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indicates however, an increase of pp breakup cross section relative to the pri breakup 
cross section as compared to the results from the low energy breakup reactions. As it 
was mentioned in Sec.III.4, at low energies 200 MeV) the cross section of pp pho- 
todisintegration from ^He is significantly smaller (by almost two orders of magnitude 
according to Ref.[83]) than the pn-breakup cross section. 

Note that the factors C and p introduce an additional uncertainty in estimating the 
magnitude of the pp-breakup cross section. While the factor C can be evaluated in the 
quark-interchange model thus staying within the framework of the considered model, 
the factor is not constrained by the theoretical framework of the model. The latter 
is sensitive to the angular dependence, F{e^,^,), of the helicity amplitudes. To esti- 
mate the uncertainty associated with F{e^,^) we varied it around the form, F{6c.m.) = 
—^73 — ^-i — jj, — in such a way that the results were still in agreement with angular 
distribution of pp scattering at u > i GeV^ We found that this variation changes the 
HRM prediction for the magnitude of pp-breakup cross section as much as 40%. 
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Figure III. 8: Dependence of the ratio of the cross section of three-body /two- 
step process discussed in Sec.V to the cross section of the two-body pp 
breakup at £'7 = 2 GeV on (a) transverse momentum of the spectator neu- 
tron pst at as = 1 and on (b) a, at p^t = 0. 



Because pp breakup cross section is still by a factor of 10 smaller than the pn cross 
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section, one needs to estimate the contribution due to three-body processes in which 
hard pn breakup is followed by soft charge-exchange rescattering. The estimate based 
on Eq. (111.42) is given in Fig. III. 8 where the ratio of three-body to two-body breakup 
cross sections is evaluated for different values of a, and transverse momentum of the 
spectator neutron, ps±. 

Because of the eikonal nature of the second rescattering in three-body processes, 
one expects the cross section to be maximal at a, = i. As Fig. III. 8(a) shows in this 
case, the three-body process is a correction to the two-body breakup process, ^2% for 
Ps± = 100 MeV/c and - 17% for p,^ = 200 MeV/c. Then, starting at Ps± i 300 MeV/c 
the three-body process dominates the two-body contribution. The latter can be verified 
by observing an onset of s"^^ scaling at large (> 300 MeV/c) transverse momenta of the 
spectator neutron in the case of hard pp breakup reactions. Figure 111.8(b) shows also 
that the three-body contribution will be always small for psi_ w MeV/c, and for a wide 
range of a^, which again reflects the eikonal nature of the second order rescattering in 
which case the recoiling of the spectator nucleons happens predominantly at - 90° (see, 
e.g., [93]). Note that one expects the above estimate of the three-body contribution to 
contain an uncertainty of 10-15%, representing the general level of accuracy of eikonal 
approximations. 

On the basis of Fig.III.8 one can expect that overall, for small values of <100- 
150 MeV/c in the high energy limit > 2 GeV) one expects two-body breakup mech- 
anisms to dominate for both pp and pn production reactions. 

Spectator nucleon momentum dependence 

The presence of a spectator nucleon in the hard two-nucleon breakup reaction from ^He 
gives us an additional degree of freedom in checking the mechanism of the photodisin- 
tegration. As foflows from Eqs.(III.31) and (III.32) and Eqs.(III.38) and (III.40) the pn 
and pp breakup cross sections within the HRM are sensitive to different components of 
the nuclear spectral function. This is a result of the fact that the pp component with the 
same helicities for both protons is suppressed in the ground state wave function of the 

91 



^He target. Thus one expects rather different spectator-momentum dependencies for pp 
and pn breakup cross sections. 
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Figure III. 9: Dependence of the s" weighted 90° cm. breakup differential 
cross section on the light-cone momentum fraction of of spectator nucleon, 
a,, calculated at = 4 GeV and = o. (a) The solid line is for pp breakup 
reactions, and the dashed line is for pn breakup reactions. Calculations are 
normalized to the cross sections at a, = i. (b) Ratio of the pn to pp breakup 
cross sections normalized to their values at a, = i. 



The quantity that we consider for numerical estimates is not the momentum of the 
spectator but rather the momentum fraction of the target carried by the spectator nu- 
cleon, as. This quantity is Lorentz invariant with respect to boosts in the q direction, 
which allows us to specify it in the Lab frame as follows: 



_ Es- Ps.z 



i2f 



(III.47) 



where a^ = ^g^^ for i = A, if, 2/ and z axis in the Lab frame is defined parallel to 
the momentum of incoming photon q. Note that the photon does not contribute to the 
above equation because a, = O. In definition of a, we use a normalization such that for 
a stationary spectator i. The a, dependencies of the differential cross sections for 
pp and pn breakup reactions normalized to their values at = i are given in Fig.III.9(a). 
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One feature of a, dependence is the asymmetry of the cross section around a, = i with 
cross sections dominating at a, > i. This property can be understood from the fact that 
the momentum fraction of the nn pair that breaks up is defined through a, as follows: 

o^NN = 3 — ag. (III. 48) 

The latter quantity defines the invariant energy of the NN pair as follows 

Snn = M'^j^ + E^rUnaNN- (III.49) 

Because the cross section within the HRM is proportional to s^]^',' it will be enhanced 
at small values of snn that will correspond to smaller values of unn or larger values of 

The difference of the cross sections because of the different composition of the 
nuclear spectral functions entering the pp and p?i breakup reactions can be seen in 
Fig. III. 9(b) in which case one calculates the ratio of pn to pp breakup cross sections 
normalized to their values at a, = i. The drop of the ratio in Fig. III. 9. (b) at values close 
to a, = 1 is the result of the suppression of the same-helicity two-proton component 
in the ground state nuclear wave function at small momenta. In this case the spectral 
function is sensitive to the higher angular momentum components of the ground state 
nuclear wave function. This yields a wider momentum distribution for the pp spectral 
function as compared to that for pn because no same-helicity state suppression exits 
for the latter. The estimates indicate that differences in a, dependencies of pp and pn 
breakup cross sections are rather large and can play an additional role in checking the 
validity of the HRM. 

III.6 Polarization transfer of the hard rescattering mechanism 

One of the unique properties of the hard rescattering mechanism of two-nucleon breakup 
is that the helicity of the nucleon from which a quark is struck is predominantly defined 

'An additional negative power of invariant energy is provided by factor in the differential cross section 

of the reaction [see Eqs.(III.3 1,111.38)]. 
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by the helicity of the incoming photon = [see Eq.(ni.l5)]. This property is a re- 
suk of the fact that in the massless quark limit the helicity of the struck quark equals the 
helicity of the photon, 77^ = A^, and assuming that at large x the quark carries almost all 
the helicity of the parent nucleon one obtains Ai; « = A^.. 

Because within the HRM, the energetic struck quark shares its momentum with 
a quark of the other nucleon through a hard gluon exchange, it will retain its initial 
helicity when it merges into the final outgoing nucleon. It will also have 
which allows us to conclude that the final outgoing nucleon will acquire the large part 
of struck quark's (as well as the photon's) helicity. This mechanism will result in a large 
(photon) polarization transfer for the hard two-nucleon breakup reactions. 

An observable that is sensitive to polarization transfer processes is the quantity C^/, 
which for a circularly polarized photon measures the asymmetry of the hard breakup 
reaction with respect to the helicity of the outgoing proton. 

A large value of C,j was predicted within the HRM for the hard breakup of the 
deuteron in Ref.[43] that was observed in the recent experiment of Ref.[60]. 

For the case of the ^He target an additional experimental observation will be a com- 
parison of c^. asymmetries for -p-p and -pn breakup channels. For the =*He target we define 
c^. as follows: 



E {|(+,A2^,A. |M|+,A^)|'-|(-,A2„A, |M|+,A.4>n 



E |(Ai^.,A2;,A, |M|+,A^)|' 



(III.50) 



Using Eq.(III.15) and the definitions of Eq.(III.20) for C,, one obtains 



(III.51) 



'+ 



where 



1 1 



d^P2,± 

(27r)2 



(III.52) 
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and s+ = S++ + s+- 



As follows from Eqs.(ni.51) and (III. 52) one predicts significandy different magni- 
tudes for c,, for pp and p?i breakup cases. 

For the pp breakup, S++ < s+- due to the smallness of the nuclear wave function 
component containing two protons in the same helicity state. As a result one expects 



ml + ml 

while for the pn break up case S++ « s+- then one obtains 

where in the last part of the equation we assumed that |(/)3| = = ^l^il. 
III.7 Summary 

The hard rescattering mechanism of a two-nucleon breakup from the -^He nucleus at 
large cm. angles is derived from the assumption of the dominance of quark-gluon 
degrees of freedom in the hard scattering process involving two nucleons. The model 
explicitly assumes that the photodisintegration process proceeds through the knock- 
out of a quark from one nucleon with a subsequent rescattering of that quark with a 
quark from the second nucleon. While photon-quark scattering is calculated explicitly, 
the sum of all possible quark rescatterings is related to the hard elastic nn scattering 
amplitude. Such a relation is found assuming that quark- interchange amplitudes provide 
the dominant contributions to the hard elastic nn scattering. 

The model allows one to calculate the cross sections of the hard breakup of pn and pp 
pairs from ^He expressing them through the amplitudes of elastic pn and pp scatterings, 
respectively. 

Several results of the HRM are worth mentioning: First, the HRM predicts an ap- 
proximate s " scaling consistent with the predictions of the quark-counting rule. How- 
ever, the model by itself is nonperturbative because the bulk of the incalculable part of 
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the scattering amplitude is hidden in the amplitude of the nn scattering that is taken 
from the experiment. 

Second, because the hard nn scattering amplitude enters into the final amplitude of 
the photodisintegration reaction, the shape of the energy dependence of the s^^ weighted 
breakup cross section reflects the shape of the s^o weighted nn elastic scattering cross 
section. Because of a better accuracy of pp elastic scattering data for s,^,^ < 24 GeV^ we 
are able to predict a specific shape for the energy dependence of the hard pp breakup 
cross section at photon energies up to < 5 GeV. This prediction and the overall 
compliance with counting rules have been confirmed by a recent experiment on =^He 
photodisintegration carried out at Jefferson Lab by the Hall A collaboration [94]. 

Another observation is that, when s " scaling is established, the HRM predicts an 
increase of the strength of the pp breakup cross section relative to the pn breakup as 
compared to the low energy results. This is the result of the feature that within the 
quark-interchange mechanism of nn scattering one has more charges flowing between 
nucleons in the pp pair than in the pn pair. This situation is opposite in the low energy 
regime when no charged meson exchanges exist for the pp pair. Even though the large 
charge factor is involved in the pp breakup its cross section is still by a factor of ten 
smaller than the cross section of the pn breakup. Within the HRM, this is due to cancel- 
lation between the helicity conserving amplitudes ^3 and ^4, which have opposite signs 
for the pp scattering. 

Because of the smallness of the pp breakup cross section, within the eikonal approx- 
imation, we estimated the possible contribution of three-body /two- step processes in 
which the initial two-body hard pn breakup is followed by the charge-exchange rescat- 
tering of an energetic neutron off the spectator proton. We found that this contribution 
has energy dependence and is a small correction for spectator nucleon momenta 
<150 MeV/c. However, the three-body /two- step process will dominate the hard pp 
breakup contribution at large transverse momenta of the spectator nucleon starting at 
p,_L > 350 MeV/c. 

The next result of the HRM is the prediction of different spectator-momentum de- 
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pendencies of breakup cross section for the pp and pn pairs. This result follows from the 
fact that the ground state wave function of =^He containing two protons with the same 
helicity is significantly suppressed as compared to the same component in the pn pair. 
Because of this, the pp spectral function is sensitive to the higher angular momentum 
components of the nuclear ground state wave function. These components generate 
wider momentum distribution as compared to say the S component of the wave func- 
tion. As a result the cross section of the pp breakup reaction exhibits wider momentum 
distribution as compared to the pn cross section. Additionally because of the strong s 
dependence of the reaction, the cross section exhibits an asymmetry in the light-cone 
momentum distribution of the spectator nucleon, favoring larger values of a,. 

The final result of the HRM is the strong difference in prediction of the polarization 
transfer asymmetry for pp and pn breakup reactions for circularly polarized photons. 
Because of the suppression of the same helicity pp components in the ^He ground state 
wave function, the dominant helicity conserving component will not contribute to the 
polarization transfer process involving two protons. Because of this effect, the HRM 
predicts longitudinal polarization transfer C^,, for the pp breakup to be close to zero. 
Because no such suppression exists for the pn breakup, the HRM predicts a rather large 
magnitude for C^, « |. 

Even though the HRM model does not contain free parameters, for numerical esti- 
mates we use the magnitude of elastic nn cross sections as well as some properties of 
the NN helicity amplitudes. This introduces certain error in our prediction of the magni- 
tudes of the breakup cross sections. For the pn breakup this error is mainly related to the 
uncertainty in the magnitude of the absolute cross section of hard elastic pn scattering 
which is on the level 30%. For the pp breakup the main source of the uncertainty is the 
magnitude of the cancellation between ^3 and (j)^, which is sensitive to the angular dis- 
tribution of helicity amplitudes. The uncertainty that results from the angular function 
is on the level of 40 %. These uncertainties should be considered on top of the theoret- 
ical uncertainties that the HRM contains due to approximations such as estimating the 
scattering amplitude at maximal value of the nuclear wave function at a i. The latter 
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may introduce an uncertainty as much as 20% in the breakup cross section. 

In conclusion, having that the HRM energy distribution for pp breakup shows good 
agreement with recent experimental data, further experimental evaluation of all the 
above-mentioned set predictions will contribute in verifying the validity of the hard 
rescattering model. Deeper insight is also expected from progress in extracting the 
helicity amplitudes of the hard nn scattering which will allow much improvement re- 
garding the accuracy of the HRM predictions. 

Mounting evidence favoring a HRM picture of A^A^ breakup reactions motivates the 
study of additional processes within the HRM framework. One such extension is dis- 
cussed in the following chapter in which a HRM description of deuteron breakup into 
AA-isobars is presented in contrast with a picture of this breakup channels emerging 
from AA components of the deuteron. It also provides an scenario alternative to that 
of a deuteron transitioning from a baryon-baryon to a six quark system as the relative 
transverse momentum of the outgoing AA system is asymptotically increased. 
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Chapter IV 

HARD BREAKUP OF THE DEUTERON INTO TWO A-ISOBARS 



In this chapter, high energy photodisintegration of the deuteron into two A-isobars at 
large center of mass angles is studied within the QCD hard rescattering model (HRM). 
According to the HRM, the process develops in three main steps: the photon knocks the 
quark from one of the nucleons in the deuteron; the struck quark rescatters off a quark 
from the other nucleon sharing the high energy of the photon; then the energetic quarks 
recombine into two outgoing baryons which have large transverse momenta. Within the 
HRM, the cross section is expressed through the amplitude of pn AA scattering which 
are evaluated on the basis of the quark-interchange model of hard hadronic scattering. 
Calculations show that the angular distribution and the strength of the photodisintegra- 
tion is mainly determined by the properties of the pn A A scattering. 

Through the HRM, the cross section of the deuteron breakup to A++A" is predicted 
to be 4-5 times larger than that of the breakup to the A+A" channel. Also, the angular 
distributions for these two channels are markedly different. These can be compared 
with the predictions derived from the assumption that two hard A-isobars are the result 
of the disintegration of the preexisting AA components of the deuteron wave function. 
In this case, one expects the angular distributions and cross sections of the breakup in 
both A++A- and A+A° channels to be similar. 

IV. 1 Introduction 

Experiments on large center of mass angle breakup of the deuteron into pn channel and 
of a pp system in -^He photodisintegration confirmed the prediction of quark-counting 
rule[51] according to which the energy dependence of the differential cross section at 
large cm. scattering angles scales as ^ s However, calculations of the absolute 
cross sections require a more detailed understanding of the dynamics of these processes. 
The considered theoretical models addressing this task can be grouped by two distinctly 
different underlying assumptions made in the calculations [76]. The first assumes that 
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the large cm. angle nucleons are produced through the interaction of the incoming 
photon with a pre-existing hard two nucleon system in the nucleus[66, 67, 96]. The 
second approach assumes that the two high momentum nucleons are produced through 
a hard rescattering at the final state of the reaction[68, 70, 71, 75, 97] as described in 
the previous chapter. 

In the hard rescattering model (HRM)[70] in particular, by explicitly introducing 
quark degrees of freedom, a parameter-free cross section has been obtained for hard 
photodisintegration of the deuteron at 90*' cm. angle [70, 71]. Also, the HRM predic- 
tion of the hard breakup of two protons from the -^He nucleus described in the previous 
chapter agreed reasonably well with the recent experimental data[94]. 

In this chapter, the HRM approach is extended to calculate hard breakup of the 
deuteron into two-A-isobars produced at large angles in the 7-d center of mass reference 
frame. Within the HRM the relative strength of 7^ A++A and A+A" cross 
sections is calculated as they compare with the -/d pn cross section. 

The investigation of the production of two energetic A-isobars from the deuteron has 
an important significance in probing possible non-nucleonic components in the deuteron 
wave function (see e.g. Refs.[99, 100, 101, 78, 102]). The study of the deuteron pho- 
todisintegration into AA-isobars channels was proposed as a venue for investigating the 
evolution of a nucleon-nucleon system into a six quark system. The onset of a six quark 
picture of the deuteron could then be marked by a large increase of the jd AA cross 
section. The latter prediction assumes that such cross section is small for a nucleon 
dominated deuteron wave function because of its suppressed AA components. In con- 
trast, for a six quark deuteron, nn and AA components contribute with comparable 
strength to the deuteron wave function (roughly 10% and 8%, respectively) while more 
than 80% is contributed by CC (hidden color) components for which unlike iv or A, C 
has a color charge. 

High energy jd AA with AA emerging at large transverse momentum is thought to 
probe the onset of hidden color components in the deuteron. Assuming that the high 
AA system was created in the initial state of the interaction, in the asymptotic limit we 
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have that 

Jt Jt ' 

Under this same assumption we also have that, 

Jt ~ Jt ' 

since both AA channels in general contribute with the same strength to the spin-isospin 
wave function of the deuteron. 

In contrast, throughout this chapter, the role of hard rescattering in these processes 
will be assessed. It will allow us also to explore another venue for checking the basic 
mechanism of high momentum transfer breakup of nuclei into two baryons. As it will 
be shown in the sections that follow, the HRM description of these reactions results in 
distinct predictions for angular distribution of the A-isobar pair at large cm. production 
angle as well as their relative strength compared with the production of the pn pair at 
the same kinematics. 

Despite experimental challenges associated with the investigation of two A-isobar 
breakup of the deuteron [104], there are ongoing efforts in performing such experiments 
at the Jefferson Lab [105, 106] which will provide the opportunity to asses the validity 
of the HRM scenario by testing its predictions discussed in this chapter. 

IV.2 Hard Rescattering Model 

We consider the photoproduction of two baryons, Bi and B2, in the reaction, 

j + d^B^ + B2 (IV.l) 

in which the baryons are produced at large angles in the 7 center of mass reference 
frame. 

According to the HRM, the large angle breakup of the NN system proceeds through 
the knock-out of a valence quark from one of the nucleons with subsequent hard rescat- 
tering of the struck-quark with a valence quark of the second nucleon. The two quarks 
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then recombine with the spectator systems of nucleons forming two emerging baryons 
with large transverse momenta. The hard rescattering provides the mechanism of shar- 
ing the photon's energy among two final baryons. 

The invariant amplitude of the photodisintegration Eq.(IV.l) is calculated by apply- 
ing Feynman diagram rules to diagrams similar to Fig.IV.l. During the calculation we 
introduce undetermined quark wave functions of baryons to account for the transition 
of the initial nucleons to the quark- spectator systems, and also for the recombination of 
the final state quarks with these spectator systems into the final two baryon system. 




Figure IV. 1 : Deuteron photodisintegration according to the HRM 

Fig.IV.l displays the chosen independent momenta for three loop integration in- 
volved in the invariant amplitude. As it was the case for ^He photodisintegration in the 
previous chapter, two major approximations simplify further calculations. First, using 
the fact that the struck quark is very energetic we treat it on its mass shell. Then the 
struck quark's propagator is evaluated at it's pole value at such magnitudes of nucleon 
momenta that maximize the deuteron wave function. These approximations allow us to 
factorize the invariant amplitude into three distinguished parts. The first, representing 
the transition amplitude of the deuteron into the (pn) system, which can be evaluated 
using a realistic deuteron wave function. The second is the amplitude of photon-quark 
interaction, and the third term represents the hard rescattering of the struck quark with 
recombination into a two large transverse momentum baryonic system. Combined with 
the initial state nucleon wave functions, the rescattering part is expressed through the 
quark- interchange (QI) amplitude of pn b^B-, scattering. A detailed derivation is given 
in section III.l in conjunction with Appendices C and D (see also Ref.[70]). After the 
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above mentioned factorization is made, the overall invariant amplitude of -fd -> B1B2 
reaction can be expressed as follows: 

(Aij, A2/ I M I \^,\d) = ie[X^] X 

+ tSt^^^/'^i/ I T(^^^B,s.).(^,i^) I A,,;A,)M/^''(p,„A,,;p2„A,)|^J(IV.2) 

where A^,Ad, \if and A2/ are the helicities of the photon, deuteron and the two outgoing 
baryons respectively. Here 'i'^''ipii,Xu;p2iA2i) is the Ad-helicity light-cone deuteron wave 
function defined in the g+ = reference frame. The initial light-cone momenta of the 
nucleons in the deuteron are pu ^ {a,, - ^,pu± = -p±) and p,,, = {a^, = |,p2>± = p±) with Ai, 
and Aai being their helicities respectively. The ^ factor with s' = s - comes from the 
energetic propagator of the struck quark before its rescattering. The squares of the total 
invariant energy as well as the momentum transfer are defined as follows: 

s = {q+p,r = {p,f+P2fr = 2E';'-M, + M', 

t = {p,,-qf = {p,,~p,f (IV.3) 

where q, pa, p^- and pa/ are the four-momenta of the photon, deuteron and two outgoing 
baryons respectively. The lab energy of the photon is defined by £'^"\ and is the mass 
of the deuteron. The transfer momentum, <jv in the rescattering amplitude in Eq.(IV.2) 
is defined as: 

= ip.f Pu - qf = {P2S - P.) - {P2S -^r = l + ^-^, (IV.4) 

where the approximation in the right hand side follows from the assumption that the 
magnitudes of light-cone momentum fractions of bound nucleons dominating in the 
scattering amplitude are au ^a2i^\, and that the transverse momenta of these nucleons 
are negligible as compared to the momentum transfer in the reaction, < |ijv|, |mjv|. 
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In Eq.(IV.2) the following expression 



Qt{^2f]^lf I T^Jn-^BtB2),i{s^'^N) \ Ki] ^2i) (IV.5) 

represents the quark-charge weighted QI amplitude of pn B1B2 hard exclusive scat- 
tering. The factor Qi corresponds to the charge (in e units) of the quark that interacts 
with the incoming photon. In a further approximation we factorize the hard rescattering 
amplitude from the integral since the momentum transfer entering in T^pn^BiB2)As,tN) 
significantly exceeds the Fermi momentum of the nucleon in the deuteron. Also, after 
calculating the overall quark-charge factors, the QI scattering amplitudes are identified 
with the NN B1B2 helicity amplitudes as follows: 

(A2/; Ai, I T«/_^^sJs,f^) I Ai.; X2.) = (IV.6) 

where 9^^ is the effective center of mass angle defined for given s and In- 
The differential cross section for unpolarized scattering is obtained through: 

d(J^d^BxB2 _ 1 1 \'\/r\2 



(IV.7) 



mia^B,B2^\\ E l(A,„A,, |M|A„A,)r, (IV.8) 



dt 167r(s-A/|)' '^''^^i^^ 

where 

LTvrP, „ „ = 

32 



with the invariant amplitude square averaged by the number of helicity states of the 
deuteron and photon. 

IV.3 Cross section of the breakup reaction 

We derive the amplitude of the breakup of the deuteron into the pn pair from Eq.(IV.2) by 
introducing the independent helicity amplitudes of pn elastic scattering Eq.(B.2) and by 
separating the quark-charge factors into and Q^^ which correspond to the scattering 
of the photon off the quark of the first and the second nucleons in the deuteron. Then, 
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for Eq.(IV.8) one obtains: 



11 

2 3 2? 

+ Ss4{\Q 



iVi, 



} 



+ 25o|(Q^^+Q^^)</)5r 



(IV.9) 



where the light-cone spectral functions of the deuteron are defined as follows: 



5*34 

So 



E E 
E E 

A=-l (Ai = -A2 = -i) 
'S'l2 + '5'34- 



(27r)2 



(IV. 10) 



Eq.(IV.9) can be further simplified if we assume (see e.g. [49]) that 03 ~ <p4, as well as 

5*12 ~ Ss4 = ^, which results in: 



= ll^Qi„„^m' + \^2\' + m' + m' + \^,\'+Mcf>s\']. (iv.ii) 

Using the expression of the differential cross section of elastic pn scattering: 

J-NN^NN ( ON ^ 1 11 

J = lkM7^)>f + + 1*1" + IM'^Mm. (IV.12) 

and the relation between the light-cone and non-relativistic deuteron wave functions [77, 
70, 107, 81] at small internal momenta: *d(Q!,p±) = (27r)i*d,]vfl(p)Vm]v in Eq.(IV.9), for 
the differential cross section on obtains from Eq.(IV.7): 

It " 7' dt ^^^-^^^ 
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where we neglected the difference between im% and m^. Here the averaged non-relativistic 
spectral function of the deuteron is defined as follows: 



^ A = l 2 
So,NR = oYl 



3 

A=-l Ai,A 



(IV. 14) 



where ^'d jvA is the non relativistic deuteron wave function, which can be calculated using 
realistic nn interaction potentials. 

The quark-charge factor, Qp^p,, = |[70] accounts for the amount of the effective 
charge exchanged between the proton and the neutron in the rescattering. It is esti- 
mated by counting all the possible quark-exchanges within the pn pair weighted with 
the charge of one of the exchanged quarks (for more details see Appendix B). The 
result in Eq.(IV.13) is remarkably simple and contains no free parameters. It can be 
evaluated using the experimental values of the differential cross section of the elastic 
pn scattering, '''"""^"rff''^""' ' • The angle e^^^ entering in the pn pn cross section is the 
center of mass angle of the scattering corresponding to the nn elastic reaction at s and 
tjv. It is related to e^,^„ of the pn photodisintegration by (See Appendix E): 



It is worth mentioning that as it follows from the equation above, 9c.„^. = 90° photodisin- 
tegration will correspond to the = 60" hard pn elastic rescattering at the final state of 
the reaction. 

IV.4 Cross section of the -fd AA breakup reaction 

We use an approach similar to that in Sec. IV. 3 to derive the invariant amplitude of the 
7(i AA reactions. In this case Eq.(IV.2) requires an input of the helicity amplitudes of 
the corresponding pn AA scattering. One has a total 32 independent helicity ampli- 
tudes for this scattering. To simplify further our derivations, we will restrict ourselves 
by considering only the seven helicity conserving amplitudes given in Eq.(B.3). Using 
these amplitudes in Eq.(IV.2) and separating the quark-charge factors into Q'^^ and (5^^ 
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similar to Eq.(IV.9) one obtains 



|M|2 



1 1 



5i 



7d^AA 2 3 2s 



\Q 



^i^.^ + g^=<^3r 



(IV. 16) 



where 5i2 and S34, are defined in Eq.(IV.lO). Similar to the previous section, we simplify 
further the above expression assuming that all helicity conserving amplitudes are of the 
same order of magnitude. Assuming also that S12 ~ 5*34 « ^, we obtain 

where gF AA = g^' + g^" = i is obtained by using the same approach as for the case of 
the pn breakup in Sec.IV.3. Using now the expression of the differential cross section of 
pn AA scattering, 

dt 



Idtt (s — 4m^j 2 



106^ + 107^ + 108^ + 109^] (IV.18) 



as well as the relation between light-cone and non relativistic deuteron wave function 
discussed in Sec.IV.3, from Eq.(IV.7) we obtain the following expression for the differ- 
ential cross section of the jd -> AA scattering: 



dt 



dt 



So.. 



(IV. 19) 



where So,nb. is given in Eq.(IV.14). The effective cm. angle entering in the argument 
of the differential cross section of pn AA reaction can be calculated by using Eqs. 
(IV.3) and (IV.4) (see also Appendix E) to obtain 



2^ 



Am%cos6 



Ami 



(IV.20) 



As it follows from Eq.(IV.19), provided there are enough experimental data on high 
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momentum transfer pn AA differential cross sections, the AA cross section can 
be computed without introducing an adjustable free parameter. However, there are no 
experimental data on hard exclusive pn AA reactions with sufficient accuracy that 
would allow us to make quantitative estimates based on Eq.(IV.19). Instead, in the next 
section we will attempt to make quantitative predictions based on the quark-interchange 
framework of hard scattering. 

IV.5 Estimates of the relative strength of the AA breakup reactions. 

The results presented in this section are calculated considering the experimental ob- 
servation [37] that the quark- interchange [108] represents the dominant mechanism of 
hard exclusive scattering of baryons that carry valence quarks with common flavor. The 
quark-interchange mechanism however will not allow us to calculate the absolute cross 
sections. Instead, we expect that its predictions will be more reliable for the ratios of 
the differential cross sections for different exclusive channels. 

As an illustration of the reliability of calculations of cross section ratios in the QI 
model, in Fig.IV.2 compares the QI predictions for the ratios of pn to pp differential 
cross sections at 90" cm. scattering. Here, we compare predictions based on SU(6)[42, 
38] and diquark (See chapter II and Ref.[110]) symmetry approaches for the valence 
quark wave function of the nucleons. As the comparison shows, one achieves a rather 
reasonable agreement with the data without any additional normalization parameter. On 
the basis of this agreement, we now estimate the ratio of the differential cross sections 
of 7rf AA to the 7d pn cross sections. We use both SU(6) and diquark- symmetry 
quark wave functions of the nucleon and A-isobars (see Appendix A) in the calculation 
of the pn AA amplitudes (see Appendix B). 

To calculate the photodisintegration amplitudes we go back to Eqs.(IV.9) and (IV. 16) 
and evaluate the quark-charge factors using SU(6) or diquark symmetries of the valence 
quark wave functions of baryons. For this we separate the t and u channels in the helicity 
amplitudes: 

Ms,o:,^,) = 0^(5,0^™.) + (IV.21) 
108 



t 

ao.8 

0.6 



0.4 

0.2 




SU(6), p=l 
Diquark, p=:0 



"A 



5 10 15 20 25 30 

s(GeVO 

Figure IV.2: (Color online) Ratio of the pn pn to pp pp elastic differential 
cross sections as a function of s at = 90". 



and then treat the charge factors for the given nucleon N as: 



Q^'cl^i = Ql 



(IV.22) 



This yields the following expression for the photodisintegration amplitude of Eq.(IV.2) 



E7b[Q*"^'^: + Q""^^r].,J*^^(;'i,A.;p,A,)|^|. (IV.23) 



7(i — > j9n scattering 



For the -fd pn amplitude, the charge factors calculated for the helicity conserving 
amplitudes according to the QI framework yield for both SU(6) and diquark models 
(see Appendix B) 



i3 - "^J 



2QJ'^^ = 2Q^„„ (IV.24) 
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with Qf,p,^ = \ and independent of j. Using tliese relations in Eq.(IV.22), from Eqs.(IV.23) 
and (IV.9) one obtains 



6i + 03 



(IV.25) 



where the different predictions of SU(6) and diquark models follow from the different 
predictions for the pn -pn helicity conserving amplitudes given in Eq.(II.12). 

7(i — > A+A^ scattering 

The calculation for the 7rf A+A'' amplitude yields the same quark-charge factors as 
for the 7cf -> -pn reactions in Eq.(IV.24). Using the helicity amplitudes of the jm A+A" 
scattering from Eq.(B.12) and the expressions for the photodisintegration amplitudes 
from Eqs.(IV.23,IV.16) one obtains 



2s' 



-5.3 



+ 20--(/.n + 



+ 2(?!)3 - 04 



+ 



+ 2C 



+ 20^ - (i>l 



(IV.26) 



where the different predictions of SU(6) and diquark models follow from the different 
predictions for the pn A+A" helicity conserving amplitudes given in Eq.(B.12). 

7(i — > A++A scattering 

For the charge factors in the 7d A++A scattering within the quark-interchange ap- 
proximation from Appendix B we obtain: 



Q 



tNi 



'— - Qfma - 3. 



(IV.27) 
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Inserting these charge factors in Eqs.(IV.23,IV.16) one obtains for the photodisintegra- 
tion amplitude: 

+ ^34 [(203 - 04)' + (204 - 03)' + 5108^ } . (IV.28) 

where predictions for the helicity conserving amplitudes of pn A++A- are given in 
Eq.(B.13). 

Numerical Estimates 

Using Eqs.(IV.25), (IV.26) and (IV.28) with the baryonic helicity amplitudes calculated 
in Appendix B we estimate the ratio R{e,,^,) of the -fd AA to 7^ pn differential cross 
sections at given s and 0^.^,, angle. For simplicity we consider the kinematics in which 
s » 4mi, which allows to approximate both Eqs.(IV.15) and (IV.20) to, 

cosez,^ « ^ + ""f— . (IV.29) 

Before considering any specific model for angular distribution, one can make two gen- 
eral statements about the properties of the photodisintegration amplitude. First, that 
from the absence of the u channel scattering in the pn A++A- helicity amplitudes (see 
Eq.(B.13)), one observes that R{e^,^,) can not be a uniform function of e^,^. Second, that 
independent of the choice of SU(6) or diquark models, the jd A++A- cross section is 
always larger than the cross section of the jd A+A- reaction. 

We quantify the above observations by parameterizing the angular function f{6'^,^), 
which enters in Eqs.(II.12,B.12,B.13), in the following form[49]: 

sin[0)^{l — cos{6)Y 

known to describe reasonably well the elastic pp and pn scattering cross sections. 

Magnitudes of the ratio R at 9c.m = 90" are given in Table IV. 1, while the angular 
dependencies (solid curves for diquark model and dashed curves for SU(6) model) are 
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Figure IV.3: (a)Ratio of the jd AA to •yd pn differential cross sections 
and (b) ratio of the jd A++A- to yd A+A" differential cross sections as a 
function of e^^ . 



presented in Fig.IV.3(a). They clearly show strong angular anisotropy and the excess 
(by a factor of 4-5) of the A++A- breakup cross section relative to the cross section of the 
A+A" breakup (Fig.IV.3(b)). These results show that the ratio of the jd AA to yd pn 
cross sections is very sensitive to the choice of SU(6) or diquark models of the wave 
functions. However, because of the absence of isosinglet two-quark state in the A wave 
functions, the p parameter dependence that characterizes the choice of SU(6) or diquark 
models in the baryons wave functions is factorized and enters only in the normalization 
factor of the pn AA helicity amplitudes. As a result, the ratio of the ^d A++A to 
yd A+A" cross scctions (Fig.IV.3b) is independent of the choice between SU(6) and 
diquark models for the baryons wave functions. 

Finally, it is worth discussing how these results compare with the predictions of 
models in which the production of two A's is a result of the breakup of the pre-existing 
AA component of the deuteron wave function. In this case, the final state interaction 
is dominated by soft scattering of two A's in the final state which will induce similar 
angular distributions for both A++A and A+A" channels (see e.g. [113, 91]). As a result, 
we expect essentially the same angular distribution for both A++A~ and A+A" production 
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channels. Also, because of the deuteron being an isosinglet, the probabilities of finding 
preexisting A++A and A+A" are equal. For coherent hard breakup of the preexisting 
A's we will obtain the same cross section for both the A++A and the A+A° channels. 





i?(90°) 


-yd^ BB 


SU(6) 


Diquark 


-fd A+A" 


0.47 


0.11 


A++A- 


2.01 


0.47 



Table IV. 1 : Strength of AA channels relative to p?2 in deuteron photodisinte- 
gration at e^,„, = 90°. 

One interesting scenario for probing the preexisting A's in the deuteron is using the 
decomposition of the deuteron wave function, in the chiral symmetry restored limit, 
into the nucleonic and non-nucleonic components in the following form[99, 100, 101]: 

*T=0,S=1 = (^)**iViV + (^)**AA + (^)*^'CC, (IV.31) 

9 45 5 

where -^cc represents the hidden color component of r = o and 5 = i six-quark con- 
figuration. Since A++A and A+A" components enter with equal probability in the total 
isospin T = configuration, one expects close 0.8) strengths for deuteron breakup to 
A++A or A+A" channels as compared to the strength of the deuteron breakup into the 
pn pair. The latter result should be compared with the similar ratios presented in Table 
IV. 1 from HRM, and with the HRM angular distributions in FigIV.3. 

In contrast to Eq. (IV.31), the pn component of the non relativistic deuteron wave 
function largely dominates over other baryon baryon components. From early works 
investigating deuteron composition by means of phenomenological and one pion ex- 
change potentials (see e.g, Refs.[l 18, 119, 120, 121, 122]), nn* components have been 
estimated to contribute overall of the order of 1% with a^a^*(1440) ^ 0.0i% , while AA 
contributions are estimated in the range 0.01-3%. These two kinds of baryonic states 
can be comparable contributions to the deuteron, thus it is expected that at low ener- 
gies the channels of deuteron breakup into nn* and AA may be comparable as well 
which, in addition to the fact that iv*(i440) can decay into a nucleon and a pion as well 
into a A-isobar and a pion, would make it difficult to differentiate experimentally. Note 
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that deuteron breakup into nn* will not interfere with the amplitude of AA production 
at large center of mass angles, since the decay products of the produced resonances 
occupy distinctly different phase spaces in the final state of the reaction. 

With pn in the S wave contributing 90% to the deuteron's wave function, the scenario 
of deuteron breakup into the pn channel at low relative transverse momentum clearly 
dominates the other baryon-baryon channels emerging from correspondent baryon- 
baryon components of the deuteron. The situation at kinematics in which the hard 
rescattering approach applies on the other hand can potentially be very different. As 
seen in Fig.IV.3(a), in the HRM the quark wave functions of the baryons play a crucial 
role in the predictions of the strenght of the deuteron breakup to AA channels relative 
to the pn channel. 

A similar result is expected if these channels are compared with the deuteron breakup 
into the nn* channels. Quark-wave functions of iv*(i440), if expanded in the form II.4, 
will introduce a different parameter p* (from the fact that N* is a radial excitation of 
n) playing the same role as the p used for nucleons. Then, the relative strength of the 
deuteron breakup into nn* (luo) channel relative to the pn channel will be determined 
by two parameters, p and p*. Without certainty on which of the considered quark wave- 
function approaches is more accurate, it cannot be concluded that a specific channel 
dominates for instance at 90° cm. in order to compare with the prediction at the onset 
of hidden color components. 

IV.6 Summary 

The hard rescattering model of large cm. angle photodisintegration of a two-nucleon 
system was extended to account for the production of two A-isobars. The HRM allows 
to express the cross section of 7^ pn and jd AA reactions through the large cm. 
angle differential cross section of pn — > pn and pn AA scattering amplitudes. 

Because of lack of experimental information on pn AA scattering, the quark- 
interchange model was further applied to calculate the strength of the AA cross 
section relative to the cross section of jd pn breakup reaction. We predicted a sig- 
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nificantly larger strength for the channel of breakup as compared to the A+A" 

channel which is related to the relative strength of the pn A++A and pn A+A" scat- 
terings. Because of the different angular dependences of these hadronic amplitudes, we 
also predicted a significant difference between the angular dependences of photopro- 
duction cross sections in A++A- and A+A" channels. 

These results can be compared with the prediction of the models in which two A's 
are produced due to the coherent breakup of the AA component of the deuteron wave 
function. In this case one expects essentially similar angular distributions and strengths 
for the A++A- and A+A" breakup channels. 
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Chapter V 
CONCLUSIONS 



The preceding chapters presented quantitative approaches that invoke QCD degrees of 
freedom in describing reaction mechanisms in nuclear and baryonic interactions. Such 
approaches were developed specifically for nucleon-nucleon elastic scattering, and for 
photodisintegration of a nn system in light nuclei. Both cases were analyzed in the hard 
kinematic regime in which it is expected that hard subprocesses among the baryons' 
constituents control the reactions thus enabling the study of the quark dynamics of the 
interaction. 

A hard scattering model of high energy pn elastic scattering was discussed in chap- 
ter II. World data evidenced the energy dependency of this reaction to be dictated by 
the constituent counting rule for a wide range of center of mass angles of scattering 
around 90°. The angular dependence on the other hand is also affected by the helicity- 
flavor constituent structure of the interacting nucleons. The structure of the nucleon was 
modeled in a quark-diquark picture of the nucleons, and it was shown that such picture 
better describes an angular asymmetry observed in the experimental data. Using these 
quark wave functions of nucleons in this quark diquark picture, such asymmetry was 
obtained as a function of a parameter (p) within the quark interchange (Qi) model. The 
parameter (p^) measures the average strength of the quark- vector diquark relative to the 
quark-scalar diquark components of the nucleon's wave function. By fitting the asym- 
metry in this model to the experimental asymmetry, it was found that 

p = -0.3 ±0.2, 

i.e., scalar diquarks contribute 90% on average to the nucleon helicity isospin structure. 
Hence, through this model, the experimental data disfavors the traditional SU{6) three 
quark structure of the nucleon in which both quark-diquark components contribute with 
the same strength, i.e., p = i, and that produces an opposite asymmetry. Also, having 
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that p < 0, indicates that the vector diquark has a negative phase relative to the scalar 
diquark component... 

In chapter III, the hard rescattering model (HRM) of nn breakup was developed for 
pp and pn breakup channels in -'He photodisintegration at large energy and momentum 
transfer. The analysis resulted in quantitative predictions for the differential cross sec- 
tions of both channels. A general feature of these cross sections is the dependence of 
the energy distributions that in the HRM arises from understanding the photodisintegra- 
tion process in terms of quark degrees of freedom which results in such energy depen- 
dence behaving according to quark counting rules (See Eq. 1.79). The latter prediction 
has been verified experimentally for deuteron breakup and recently for pp breakup in 
'^He at 90° cm. of the 7 nn system at beam energies larger than 2GeV^ The onset 
of QCD degrees of freedom marks also the dominance of a one step/two body process 
picture of the reaction. In a hadronic description, this process is suppressed from not 
having charged meson exchanges between the two protons. For pp break up, these ex- 
changes are only possible if the neutron participates. The reaction is then dominated 
by a three-body /two- step process. Such process is much largely suppressed at the onset 
of a explicit quark description of the reaction such as the HRM; its energy dependence 
falls off faster than s~i2. Calculations for additional observables were also carried out 
within the HRM resulting in further numerical predictions on spectator nucleon's mo- 
mentum distributions (See Fig.III.9) and the polarization transfer asymmetry C^,. In 
particular for pp it was found that C,, is almost canceled, while for pn breakup, C^- « |. 
The experimental evaluation of these predictions will properly constrain the validity of 
the HRM approach for these reactions. 

The hard rescattering mechanism was discussed in chapter IV as it applies to double 
A-isobars production in deutron photodisintegration. These studies were motivated by 
the focus on identifying a transition of the deuteron from a baryon-baryon system to a 
six quarks system. It is believed that such transition is signaled by a sudden increase of 
the double A-isobars production in relation to the deuteron breakup in the pn channel. 
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The latter is believed to dominate at small to moderate transverse momenta from the 
dominance of the pn component of the deuteron over the negligible AA components, and 
from the assumption that the baryons emerge from a large relative transverse momenta 
component of the deuteron wave function. 

Such assumption was contrasted in this chapter in which the breakup reaction was 
studied instead within the framework of the QCD hard rescattering model (HRM). The 
baryons emerging with large transverse momenta are produced in a pn rescattering pro- 
cess triggered by the incident photon interacting with a low relative momentum pn com- 
ponent of the deuteron wave function. This rescattering was modeled within the quark 
interchange mechanism (QIM), and within the HRM, angular distributions for the two 
AA channels were obtained and compare to each other and to the breakup to pn channel. 
From such comparisons illustrated in Fig.IV.3, it was conclude that in the rescatter- 
ing picture each A A channel has a distinct angular distribution with A++A dominating 
A+A . This is in clear contrast with what is expected if the two A-isobars emerge from 
a AS component of the deuteron in which case both AA channels should have the same 
strength. 

The results concluded above were arrived to from considering the studied reactions 
to be controlled by elementary particle subprocesses. Such an assumption was justi- 
fied by the kinematic characteristic of hard processes, i.e., the energy and momentum 
invariants were much larger than the masses of the interacting baryons which facili- 
tated for instance a short distance-long distance factorization of scattering amplitudes, 
at a quark-baryon level as it is discussed in chapter II, and at the quark-baryon-nuclear 
level which corresponds to the hard rescattering model developed in chapters III and 
IV. This methodology finds a broad range of applications in many more reactions of in- 
terest. Some of them are natural extensions of the processes studied here. For instance 
a program to pursue, as it was pointed in chapter II, is the further development of the 
studies presented there as it concerns to the possible strange production channels in ex- 
clusive eN or NN scattering, which would make full use of the SU{6) symmetry of quark 
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and diquark states in constructing quark wave-functions of hadrons providing also a 
potential probe of strange quark distributions in nucleons. This program can be as well 
continued into studying strange production channels off nuclear targets in breakup re- 
actions, and further probing the role of hard rescattering mechanisms analogous to such 
studied in chapters III and IV. The experimental difficulties brought on by the decays 
of the produced hadrons in these processes can soon be overcome by updated capabili- 
ties such as the 12GeV upgrade program at Jefferson Lab, which is also instrumental in 
assessing most of the results presented in this dissertation. 

Another venue of interest, and also a mayor motivation driving JLAB's 12GeV up- 
grade, is the study of near threshold J/V- production reactions off proton and nuclear 
targets. These reactions are relevant in identifying the leading gluon exchange mech- 
anisms between hadrons, as well as in extracting gluon contributions to the structure 
functions of nucleons at large x. Because of the large mass of the c quark(-1.2GeV), 
the cc fluctuation that will evolve into a J/ijj meson for instance in 7^ j/ij^p has a small 
transverse size, and since the threshold invariants are large as well, pQCD approaches to 
these processes may apply ( see e.g. Refs. [115], and [116]), and quantitative descrip- 
tions can be constructed through factorization methods analogous to those developed 
throughout this dissertation. Such is the case also for j/ip production off a deuterium 
target which study is of interest in probing the hidden color components of the nuclear 
wave function; in a pQCD picture of the j/^ppn reaction's hard subprocess, the 
cc fluctuation of the photon can interact by single gluon exchanges with non-singlet 
quark clusters within the nucleus [117]. This mechanism is expected to give a signif- 
icant contribution in this reaction. Also,in analogy to what was done in chapter IV, in 
an alternative mechanism, this reaction proceeds through cc fluctuation scattering off 
one of the nucleons, then evolving into a hadronic state and rescattering with the sec- 
ond nucleon into an outgoing nucleon and a J/i' meson. Thus, a variation of the hard 
rescattering model can as well be developed for obtaining qualitative and quantitative 
predictions for these heavy quark production reactions. 

In summary, these phenomena and several more await for a coordinated program 
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focusing experimental and theoretical developments into identifying explicit signatures 
of quark-gluon dynamics of hadron interactions that shed light into building a compact 
QCD description of the strong force. It is one's expectation for the methods utilized 
and the results reported in this dissertation to further encourage research efforts in such 
program in which this methodology can find applicability. 
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APPENDICES 
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A Baryonic Wavefunctions 



Nucleon and A-isobar wave functions are built from the minimal Fock component of 
the corresponding wave function, initially assuming a quark-diquark expansion where 
individual quark wave functions correspond to SU(6) eigenstates. A single quark state 
of spin S=l/2 and isospin 1=1/2 is joined by a diquark state with S=0 and 1=0, or S=l and 
1=1. The former corresponding to what is known as a scalar diquark [qq], and the latter 
known as a vector diquark qq. The baryonic wave functions are then expanded in q[qq] 
and q(qq) states with the proper Clebsh-Gordan coefficients and normalization, and 
introducing a parameter p that determines the relative amplitude of the vector diquark 
with respect to the scalar diquark sector of the expansion: 

ip^"^'^ cti q[qq]+ px qqq (A.l) 

q[qq] and q(qq) expansions vary according to the baryon's isospin and helicity. For 
instance, for A-isobar wave functions q[qq]=0 since the total isospin components won't 
add up to 1=3/2. A nucleon wave function is expanded as follows, 

V'-"''^" E ^^M.^C'"A,.„/'''tof'' (A.2) 

i3,h=l/2,-l/2 t3,A=l,0,-l 

where i%, i^,t^, are the third components of the isospin of nucleon, single quark and 
diquark respectively, and h^, h and A are the corresponding helicities. Cf^ .3^ and C^;, 
are the Clebsh-Gordan coefficients for isospin and helicity expansions in q(qq) states 
respectively. 

Since [qq] and (qq) are singletxisosinglet and tripletxisotriplet representations of qq 
states in order to ensure a symmetric wave function, Eq.(A.2) can be expanded in qqq 
states. Doing so for the case of a proton with positive helicity yields, 

N(p)|p(l/2)) = ^l\u{+)u{+)d{-)) - ^—l\u{+)d{+)u{-)) - p\d{+)u{+)u{-)) 
+ ^-^H+)d{-)u{+)) - ^-^H-^)u{~)d{+)) ~ p\d{+)u{^)u{+)) 
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+2p\d{-)u{+)u{+)) - p\u{^)u{+)d{+)) - p\u{~)d{+)u{+)) (A.3) 



where, u and (+) represent q with i'=l/2 and h=l/2 , and d and (-) represent q with 
i^=-l/2 and h=-l/2 respectively, and 

N{p) = \ ^ (A.4) 

All other nucleon wave functions can also be obtained by properly applying spin or 
isospin ladder operators starting from Eq.(A.3), e.g., 

|n(l/2))=r_|p(l/2)) (A.5) 

where 

T-ki9293) = 1(^-91)9293) + 191(^-92)93) + |9i92(T-93)) (A.6) 

and r- is constructed from Pauli matrices acting on the quark isospin state. Similarly, 
all A-isobar wave functions can be constructed through ladder operators in 1=3/2 and 
J=3/2 representations starting with A++ with helicity hA=3/2, 

|A++(3/2)) = \u{+)u{+)u{+)) (A.7) 

then for instance, to obtain A++ with helicity hA=l/2, 



|A++(l/2)) = i=a-|A++(3/2)) (A.8) 



or in general, 
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Likewise, 



with a and r acting on the three-quark components of the expansion as indicated in 
Eq.(A.6) on the helicity and isospin states respectively of the single quark states. 

A-isobar wave functions 

Either through the process described by Eqs.(A.9) and (A. 10), or by explicitly ex- 
panding each A-isobar state in q(qq) states and then qqq states through the use of 
Clebsh-Gordan coefficients, the following A-isobar helicity-isospin wave functions are 
obtained. 



|A++(3/2)) 
|A++(l/2)) 
|A++(-l/2)) 
|A++(-3/2)> 



\u{+)u{+)u{+)) 

-^\u[^)u{+)u{+)) + l=\u{+)u[-)u{+)) + -^\u{+)u[+)u{-)) 
l^\u{-)u{-)u{+)) + l^\u{-)u{+)u{-)) + l^\u{+)u{-)u{-)) 

\u{-)u{-)u{-)) 



|A+(3/2)) 
|A+(l/2)) 



|A+(-l/2)) = 



|A+(-3/2)) = 



l=\d{+)u{+)u{+)) + l^\u{+)d{+)u{+)) + l^\u{+)u{+)d{+)) 
h^d{-)u{+)u{+)) + h^d{+)u{-)u{+)) + 

+ ^\u{-)d{+)u{+)) + \\u{+)d{-)u{+)) + l\u{+)d{+)u{-)) 
+ l\u{-)u{+)d{+)) + ^-\u{+)u{-)d{+)) + ^-\u{+)u{+)d{-)) 
^-\d{-)u{-)u{+)) + l\d{-)u{+)u{-)) + l\d{+)u{~)u{-)) 
+ l\u{-)d{-)u{+)) + + l^\u{+)d{-)u{-)) 

+ ^-\u{-)u{-)d{+)) + ^-\u{-)u{+)d{-)) + ^-\u{+)u{-)d{-)) 
l^\d{-)u{-)u{-)) + l^\u{-)d{^)u{-)) + l^\u{-)u{-)d{-)) 



(A.11) 
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A° and A" wavefunctions are obtained by the replacement u^d in A+ and A++ wave- 
functions respectively. 
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B Helicity Amplitudes in the Quark Interchange Model 



B.l Baryon-Baryon Scattering Helicity Amplitudes 

We are using helicity states to label the entries of the photodisintegration and the baryon- 
baryon scattering matrices. The number of independent helicity amplitudes for a given 
ab cd processes can be expressed through the total spin of the scattering particles as 
follows[lll, 112]: 



where is the total spin of particle i and for the photon we replace (si+i) by 2. The factor 
i follows from the constraint due to the parity conservation. For elastic scattering, there 
is a further reduction in N due to time reversal invariance, and if the scattering particles 
are identical, or lie in the same isospin multiplet, the number of independent helicity 
amplitudes is reduced further[lll, 112]. For the pn elastic scattering case, out of the 
possible 16 helicity amplitudes only five are independent[112] for which we use the 
following notations: 



For the pn -> AA scattering amplitude, we have from Eq.(B.l), iV=(2)(2)(4)(4)/2=32 
independent helicity amplitudes. We use the following notations for the helicity con- 



N=-- {2s, + l){2s, + 1)(25, + l)(25rf + 1) 



(B.l) 




(B.2) 
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serving independent amplitudes of pn A A scattering: 

01 

03 
04 
06 
07 
08 
09, 

(B.3) 

which are consistent with the definitions in Eq. (B.2). 

B.1.1 Helicity Amplitudes in the Quark-Interchange Model 
Quark Interchange model 

Following the approach presented for example in Refs.[108, 42, 38], the scattering am- 
plitude for a process ab cd, in which a, b, c and d are baryons, is obtained from, 

{cd\T\ ab) ^ ^{iPl \ a'2,p[,j[){iPl \ a[,l3'^,i^) 

x(Q;2,^2>72>a'i/^hi I H \ ai, ft, 71, 02/3372) • (ai,ft,7i | Va) (0:2, ft, 72 I ipb), (B.4) 

where (ai,a[), (ft, ft') and (7^7;) describe the spin-flavor quark states before and after the 
hard scattering, H, and 

Ci,,,, = {a,f5,j\^,) (B.5) 

describes the probability amplitude of finding an a,/3,7 helicity-flavor combination of 
three valence quarks in the baryon j. These coefficients are obtained from the expansion 
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of the baryon's spin-isospin wave function in three-quark valence states as follows: 

1 1 I I , 1 , . 1 , 

><ix^.X^,^.,J-irl'^r^\,_^)}. (B.6) 

' -^-5 2 ' 23 -^'''23 2 ' iV 23 J 

The indexes 1 and 23 label the quark and the diquark states. The first term corresponds 
to quarks 2 and 3 being in a helicity zero isosinglet state, while the second term corre- 
sponds to quarks 2 and 3 in helicity 1-isotriplet states. Where x and r represent helicity 
and isospin states with helicity h and isospin projection respectively. For the wave 
functions of A-isobars a = o and p = i, while for nucleon wave functions a = i and the 
parameter p characterizes the average strength of the isotriplet diquark radial state rela- 
tive to that of the isosinglet state. Two extreme values of p = i and p = o correspond to 
the realization of the SU(6) and good diquark symmetries in the wave function. 
Using Eq.(B.6) in Eq.(B.4) for the hadronic scattering amplitude one obtains: 



where 

which accounts for all possible interchanges of a and a' quarks leaving /? and 7 quarks 
unchanged. In the QI model the interchanging quarks conserve their corresponding 
helicities and flavors, this is accounted for in the matrix elements of A in Eq.(B.7.), 

Eq.(B.7) has two terms, first (referred as a t term) in which four quarks scatter at 
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angle and two (interchanging) quarks scatter at tt 6*^^ and the second (referred as 
a u term) in which two interchanging quarks scatter at e^^^ , while four spectator quarks 
scatter at tt - 6*^^ . 

B.1.2 Helicity Amplitudes in the Quark Interchange Model 

Through the above procedure using Eq.(B.7) for the helicity amplitudes of pn scattering 
one obtains: 



(2 - + (1 + 2y)/(7r - C„J (B.IO) 

'3&.) = (2 + y)/&.) + (l+4y)/(7r-C„J 

2ym^„.) + 22//(^-C™.) 

0, 



were. 



y = 



2 p f 2 p 



31 + p^ \ 3 1 + p2 
For pn A+A° scattering amplitudes we obtain: 



^^aa(2/(C™.)-/(^-^.~™.)) 
^A^AA(4/(C„J + /(^-^r'™.)) 
^A^AA(/(C^.)) + /(7r-C™,) 



1 + • (B.ll) 



"^6 



" - ^""(2/(^?-„.)-/(7r-C..,)) 



(2/(C™.)-/('r-e:„.)) 



3V3 

A^AA 

^A^aa/(^?."^.) 

^A^AA/(^-e,^™.), 



(B.12) 
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and similarly for the amplitudes of the pn A++A- scattering, QI model gives: 



— NaA r/nN \ 



cm. 



N . 
cm. / 



0, 



(B.13) 



For both sets of equations in (B.12) and (B.13), we have 

iVAA^^, (B.14) 
1 + 

which shows that the strength of the two AA channels relative to each other is indepen- 
dent of the value of p. This is not the case for their strengths relative to the pn channel; 
from Eqs. (B.ll) we see that the p dependence of the helicity amplitudes in pn pn 
cannot be factorized. 



B.1.3 Quark-Charge Factors 

In the hard rescattering model, photodisintegration amplitudes are expressed in terms 
of hadronic scattering amplitudes weighted by the charges of struck quarks, Eq.(IV.5). 
We further split the amplitude of Eq.(IV.5) into t and u channel scatterings: 
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where Q*^"^ is the charge of the quark, struck by the incoming photon from the nucleon 
N with further e^^ or n-e^^ scattering. The helicity amplitudes are also split into t and 
u parts 

= cJ(C™.) + c„/(7r-C„.), (B.16) 

with 0* and (/>" corresponding to the 6i^,„ or tt - e^^ scattering terms in Eq.(B.7). 

Using the above definitions and Eqs.(B.7,B-8,B.ll,B.12,B.13) the charge factors Q* 
and Q" are calculated using the following relations: 



rbd 



where summation is understood for repeated a indices, Q(a) is the charge in e units of a 
quark a and the index j labels the process ah cd. 
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C Calculation of the ^He(7, nn)n scattering amplitude 



Applying Feynman diagram rules for the scattering amplitude corresponding to the di- 
agram of Fig.(ni.l)(a) one obtains 

{XfiAf2,K I M I A^, Aa) = 

(iVl) : / ^!£i^i^^4l^^L±^z^(fc,) . . . HgTf^,^ - - - ''^^ + 



en.) : / 



{Pnn - P2iY - ml, + ie pli - ml, + ie {2ttY 
[iP2^ - - (Pu - h) - {q - l)Y + ie' 

(C.l) 



where the momenta involved above are defined in Fig.III.l. Note that the terms above 
are grouped according to their momenta. As such they do not represent the correct 
sequence of the scattering presented in Fig.III.l. To indicate this we separated the 
disconnected terms by "• • •". 

The covariant vertex function, r^He describes the transition of the -^He nucleus to a 
three-nucleon system. The vertex function describes a transition of a nucleon to one- 
quark and a residual spectator quark-gluon system with total momentum fc^, (i = 1,2). 
The function S{k) describes the propagation of the off-mass shell quark-gluon spectator 
system of the nucleon. As is shown below, this nonperturbative function can be included 
in the definition of a nonperturbative single quark wave function of the nucleon. 

Using the reference frame and the kinematic conditions described in Sec. III. 1 we 
now elaborate each labeled term of Eq.(C.l) separately. 

(^He)-term. Using the light-cone representation of four-momenta and introducing the 

light-cone momentum fraction of the nn pair carried by the nucleon 2i as q = 

one represents the nucleon propagators as well as the momentum integration d^p^i in the 
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following form; 



2 2- / P^iJ- \ 
P2i- '^N = 01- PNN+[P2i- ) + »e 

apNN+ 

{Pnn -P2if +ie^ Pnn+{^ ~ a){ — — - P2i-) ~ {m% + p^i^) + ie 

Pnn+ 

d'^P2i = PNN+^dadp2i-d^P2i±. (C.2) 



Using these relations in Eq.(C.l) we can integrate over dpii- taking the residue at the 
pole of the (2z)-nucleon propagator, i.e., 



- = -2m[...]\ , . (C.3) 



After this integration one can use the following relations in Eq.(C.l): 

j4i + = ^ Ux.^,{p2i)U\2,iP2i) 



iPNM - P2.f - ml = (1 - a)(M^^ - ^N +Pli± ^ 

a[l — a) 



M 



2 _ "'w+P2i^ 
NN a(l-a) 



pInn -~pl2i + mjv = V'uAi.(pii)wAi.(pii) H 7. 7"^- (C.4) 

Zpf^i^j^^ 

Furthermore we use the condition p^jv+ > K^^atat - '"l^i^-!')^ ) to neglect the second term 
of the right-hand part of the third equation in Eq.(C.4) . This relation is justified for 
the high energy kinematics described in Sec.III.l as well as from the fact that in the 
discussed model the scattering amplitude is defined at a « i. 
Introducing the light-cone wave function of ^He [77, 78, 81] 

,T,A^,Ai,A,,A,, s UxAPnN -p)Ux^{p)uxAPs)^m^ 



M 



NN a(l-a) 



and collecting all the terms of Eq.(C.4) in the (=*He:) part of Eq.(C.l) one obtains 

{'He:)= 2^ / Ma,i(Pi)ma.,(P2)— .3 - (CO) 

^ — ' J I — a a z(27rj^ 
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(Nl:). To evaluate this term in Eq.(C.l) we first introduce 



Pii+ (1 - oi)pnn+ ' 
_ _ I- a -yN 

Pif+ 1 - a 

where a' = Furthermore we perform the fci integration such that it puts the spec- 
tator system of the A^i nucleon at its on-mass shell. This will results in 

where Vs(fc) represents the nucleon's spectator wave function with mass m„ and spin s. 
Note that in the definition of tp, one assumes an integration over all the internal momenta 
of the spectator system. Using Eq.(C.8) for the (Nl) term one obtains 

Now we evaluate the propagator of the off- shell quark with the momentum, pu - h. 
This yields: 



1 



^P^^ ™9 (l-xi)(m^i ) 



7+ 



2(1 - a){l - Xi)pnn+ ' 

where the effective off-shell mass of the nucleon is defined as 

M^^q;(1 - a) - 777^(1 -a)-p\ 



(C.IO) 



(C.ll) 



As it follows from Eq.(C.lO) at the high energy limit, p%j^^ > m%, one can neglect the 
second term of the RHS (off-shell) part of the equation if (i - a){i - x^) ^ i. As is shown 
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in Sec. III. 1 [see discussion before Eq.(III.9)], the essential values that contribute in the 
scattering amplitude correspond to a « 1 and (i - xi) ^ i. Therefore the second term in 
the right-hand side part of Eq.(C.lO) can be neglected. Using the closure relation for 
the on-shell spinors for Eq.(C.lO) one obtains 

j/u ™9 _ E^i, U^^,{Pl^ - kl)u^,^{pu - fcl) (C 12) 



(Pu Kij {l-x,){m%, ^r^^^^ ) 

Similar considerations yield the following expression for the propagator of the quark 
entering the wave function of the final nucleon "i/": 



where x[ is defined in Eq.(C.7). 

By inserting Eqs.(C.12) and (C.13) into Eq.(C.9) and defining quark wave function 
of the nucleon as 

,T,A,^/„ ^ /, N , U^ip-k)^Pl{k)TMU^^{p) (C]4) 

for the (ivi :) term we obtain 

(^10 2^ J {l^x[) ^vifiPif-ki)--- 



(N2:). This term can be evaluated following similar considerations used above in the 
evaluation of the (Nl:) term. Introducing light-cone momentum fraction of the specta- 
tor system of the second nucleon as 

P2i+ apNN+ 

(C.16) 



k2+ 


^2+ 


P2i+ 


apNN+ ' 




a 




= —X2 


P2f+ 


a' 
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for the (N2:) term we obtain 



(iV2 :) 



1^ I n— ;m ""^^-Mi 

V2f •V2i •S2 



(1 - X',) 

*^^"''^'(P2», 2:2,^21) dX2 (Pk 



(1-2^2) 



X2 2(27r) 



(C.17) 



Collecting the expressions of Eqs.(C.6), (C.15) and (C.17) in Eq.(C.l) and rearrang- 
ing terms to express the sequence of the scattering, we obtain the expression of the 
scattering amplitude presented in Eq.(III.4). 
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D Calculation of the nucleon-nucleon scattering amplitude 



In this section we consider a hard nn elastic scattering model in which two nucleons 
interact through the(QIM). The typical diagram for such scattering is presented in Fig. 
D.l. Applying Feynman diagram rules for these diagrams one obtains 

a) 




b) 




Figure D.l: Quark interchange contribution to nucleon-nucleon scattering 



rpQlM __ 

^ NN — 



{pif - fci)2 - ?7i2 + ie 
{N2) : / —^^^lJ!^I-^!l—^S{k2) 



{P2f - ^2)^ - m'^g + ie 



-iPlf ^P2f), 



(D.l) 
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where definitions of the momenta are given in Fig.(D.l). The procedure of reducing the 
above amplitude is similar to the one used in the previous section. First we estimate 
the propagators of each nucleon's spectator system at their pole values, ki^. = ""l^Jl'^ 
(i = 1, 2) by performing the fc,, integration, which yields 

f[...]s{h)dh_^~'^^y2i'^{h)^i{h)\ . (D.2) 

Furthermore, because v\ » m% one can apply similar to Eqs.(C.12) and (C.13), ap- 
proximations for propagators of interchanging quarks leaving and entering the corre- 
sponding nucleons. Then using the definition of single quark wave function according 
to Eq.(C.14) for the (iVi) and {N2) terms, one obtains similar expressions that can be 
presented in the following form: 

[Nl:) ^ J (l~x') ^mfiPif - ki) ■ ■ ■ 

• • • ^^KAPu - fci) —2(2^- (^-3) 

The {N2 :) term is obtained from the above equation by replacing 1^2. Regrouping (Nl) 
and (N2) terms given by Eq.(D.3) into Eq.(D.l), for the amplitude of nucleon-nucleon 
scattering in QIM we obtain 



VliV2iVlfV2f 



X 



^. .dxi (Pki± dx2 d^k'zi 



xi 2(27r)3 X2 2(27r)3 



(D.4) 



Note that in the above expression we redefined the initial momentum of "Afi" nucleon 
to pi^ to emphasize its difference from which enters in the photodisintegration am- 
plitude. In the latter case p^i is not independent and it is defined by the momenta of the 
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two remaining nucleons in the ^He nucleus. 
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E Kinematic Relations in -fNN bb 

For the process 7 + nn B1 + B2 we define, 

S = {q+PNNf = (Plf +P2ff 

t = (Plf - lY = {P2f - PmnT (E.l) 

In our approximation, pjvat = = 2^2; • Then, 

t = m%,^ + M^^ - 2p2fPNN (E.2) 

The two emerging baryons are producing from the rescattering of the nucleon absorbing 
the incoming photon and the second nucleon in the nn system. For this rescatering we 
define. 



tN = {Plf - Pli - qf = {P2f - P2if 

I Pnn\ 

= [P2f — I 



2 



2 2 4 



2 / 

f2 

(E.3) 



where the last equality is obtained using Eq.(E.2). 
For the case in which mg^ = ms, = mg we have that. 



2 s-M, 

t — m„ „ _ 

2yi 



2 



- a/s - Am%cos9^,,n^ , (E.4) 



where 6*^ is the angle of scattering in the center of mass reference frame of the ^nn 
system. Then from Eq.(E.3), 



tj^ — nig 



2 Ml^ s - Ml^ 



4yi 



- - 4^m%cos9^,^^ (E.5) 
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We also have that, 



+ ml + ml + \^J{s - Ami) (s - '^mDcosO'^^^^, 



(E.6) 



where 6*^^ is the angle of scattering in the center of mass reference frame of two nucle- 
ons scattering into two baryons at a center of mass energy of E^,^ = ^. 9^.^. and can 
be related to each other through Eqs. (E.5) and (E.6), 



^(s - Am%) (s - 4m|) 



■ 



1\J. A 



2Vi 



— \/ s — Am%cos9i. 



.(E.7) 



Then for instance if the final state baryons emerge from 7 + nn Bi+ B2 at 9^,,^, = 90° 
then in the corresponding N + N ^ Bi + B2 process. 



9c.m. = 90°) = arcos 



1 s — Ami, 



2 y s - 4m| 



(E.8) 
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F Parametrization of nn Elastic Scattering Experimental Cross Sections 



Unpolarized nucleon-nucleon elastic scattering data is used to parametrize ^ and 
such that they can be used in calculations of photodisintegration of nucleon pairs ac- 
cording to the hard rescattering model described above. For proton-proton, elastic scat- 
tering center of mass energy distributions are fitted for different center of mass angles 
of scattering through the following parametrization: 



= 77^) e'^^^^Ris^e^,,^) (F.l) 



dt ' ' ' KlOGeV^ 



in units of in which s is input in GeV^. With R{s,9^^ introducing an oscillation 
factor in a fixed 9^^ energy distribution, 



R{s, e,,^,) = i?„ ( 1 + as''cos{ujln{ln{^)) + 6} + ^a'^s"" ) , (F.2) 



in which the parameters adjusted for different center of mass scattering angles are given 
in Table F. 1 . 

^(s,6'„„) is then obtained by interpolating at the corresponding center of mass en- 
ergy s and center of mass angle of scattering e^,^, through the fits of energy distributions 
at the selected angles in Table F. 1 . 

An alternative parametrization follows the form. 



-^(S, ^'c.m.) = (Si^^'c.m.) ^'^ R{s,0c.m.)^^a„s", (F.3) 

11=0 

in which the parameters on R{s, 6'„„) and 7 are now fixed to R^ = 4.5 x 10*, a = 0.O8, k = 0.5, 
uj = = 0.1, and 6 = -2, and s is replaced by s + 2GeV^ in R's argument for s > 2{)GeV^. 
The polynomial expansion coefficients a„ are given in Table F.2. For s < 20GeV'^ at e'^ ^ s 
between 55° and 90°. For s > 20GeV'^ and n > 0, a„ = while ^ 1. Then the polynomial 
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Fitting parameters 






am 


ln{sp) 


a 






A 


5 


90° 


0.055 


0.015 


3.22 


0.120 


0.3 


52.36 


0.10 


-1.57 


85° 


0.065 


0.042 


3.10 


0.100 


0.3 


62.83 


0.10 


-4.71 


80° 


0.077 


0.050 


3.05 


0.090 


0.3 


62.83 


0.10 


-4.71 


75° 


0.100 


0.060 


3.02 


0.090 


0.3 


62.83 


0.10 


-4.71 


70° 


0.120 


0.060 


3.01 


0.150 


0.3 


57.12 


0.09 


-1.57 


65° 


0.220 


0.067 


3.01 


0.080 


0.3 


54.16 


0.09 


3.14 


60° 


0.430 


0.065 


3.20 


0.080 


0.3 


55.12 


0.09 


-5.50 


55° 


0.700 


0.060 


3.18 


0.085 


0.3 


62.83 


0.09 


-3.14 


50° 


2.100 


0.033 


3.50 


0.200 


0.3 


31.42 


0.09 


-2.86 


45° 


5.200 


0.033 


3.55 


0.200 


0.3 


31.42 


0.09 


-3.46 


40° 


20.00 


0.020 


3.88 


0.200 


0.3 


31.42 


0.09 


2.42 


35° 


70.00 


0.023 


3.88 


0.280 


0.3 


28.56 


0.09 


1.85 


30° 


400.0 


0.027 


3.88 


0.280 


0.3 


19.63 


0.09 


0.78 


25° 


4000 


0.035 


3.88 


0.330 


0.3 


16.53 


0.09 


0.78 


20° 


15000 


0.030 


4.00 


0.200 


0.3 


16.53 


0.09 


0.78 



Table F.l: i?(s, 9^.^.) (Eq.(F.2)) parameters to be used in Eq.(F.l) for elastic 
scattering 

factor is interpolated in 9^.,,^, to find ^(s,6'„„) for s < 20GeV'^. 





Fitting parameters 


n 

"cm 








as 


a4 


90.0° 


-5.6155 


1.9297 


-0.21216 


1.0073 X 10-2 


-1.7088 X 10--* 


82.5° 


-9.2744 


3.2989 


-0.39409 


2.0202 X 10-2 


-3.6891 X 10-* 


77.5° 


-13.967 


5.0478 


-0.62805 


3.3308 X 10-2 


-6.2775 X 10-* 


72.5° 


-10.865 


3.8837 


-0.47725 


2.5176 X 10-2 


-4.7326 X 10-* 


67.5° 


-6.2905 


2.2328 


-0.27017 


1.4286 X 10-2 


-2.7048 X 10-* 


62.5° 


-8.7180 


3.2655 


-0.42545 


2.3441 X 10-2 


-4.4937 X 10-* 


57.5° 


-10.841 


4.0907 


-0.54142 


3.0252 X 10-2 


-5.8992 X 10-* 



Table F.2: i?(s, 6'c.,„.) (Eq.(F.2)) parameters to be used in Eq.(F.3) for pp elastic 
scattering 

Figures (F.l) and (F.2) show angular distribution and energy distribution fits respec- 
tively on world data. 

Because of fewer experimental data, for proton-neutron, is parametrized by 
fitting angular distributions at fixed incident momentum, through polynomials on cos9^,,^ 
so that. 
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(S, 0^rn) = {sinO^rn) ^ ^ c" {cOsO , (F.4) 



in which the coefficients c" of the sum term are given in Table F.3 for selected values of 
incident momenta. 





Fitting parameters 


PLab{GeV/c) 


Co 


Cl 




cz 




3 


1.7700 X 10-1 


3.1775 X 10-1 


-2.8763 X 10-1 


-2.8763 X 10-1 


3.6249 X 10-1 


4 


2.5404 X 10-2 


4.4345 X 10-2 


-2.1687 X 10-2 


-5.8543 X 10-2 


1.5246 X 10-2 


5 


4.3600 X 10-3 


7.5730 X 10-3 


2.7010 X 10-3 


-2.8781 X 10-3 


-2.4397 X 10-3 


6 


1.2508 X 10-3 


2.5914 X 10-3 


1.6694 X 10-3 


-2.8558 X 10-3 


-2.8406 X 10-3 


7 


3.0720 X 10-^ 


5.9495 X 10-* 


8.7070 X 10-* 


-5.0500 X 10-* 


-1.2625 X 10-3 


8 


6.6280 X 10-5 


1.4567 X 10-* 


3.4920 X 10-* 


3.0956 X 10-5 


-2.8817 X 10-* 


9 


1.9080 X 10-5 


3.5100 X 10-5 


1.3483 X 10-* 


1.3052 X 10-* 


1.7023 X 10-* 


10 


6.9002 X 10-«* 


2.0643 X 10-5 


8.7424 X 10-5 


3.1467 X 10-5 


-4.9238 X 10-5 


12 


2.5381 X 10-'^ 


1.1734 X 10-^ 


4.5911 X 10-" 


3.1104 X 10-5 


3.3758 X 10-5 



Table F.3: Expansion coefficients of the parametrization for pn elastic scat- 
tering of Eq.(F.4). 
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0.1 




Figure F.2: Fits from Eq.(F.l)(solid line) and Eq.(F. 3) (dotted line) to pp en- 
ergy distribution data at 6^,^, = 90° 
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G Helicity Amplitudes of Photodisintegration 



G.l Helicity Amplitudes in Deuteron Breakup 

Using the notation described in appendix A for helicity amplitudes in pn elastic scat- 
tering and working out Eq.(IV.2) for Bi=p and B2=n we obtain the following scattering 
amplitudes for 



i +1|M|+,A, 



1 1 

2'"2 



1 1 

2'"2 



|M|+,Ad 



|M|+,Ad 



-i +i|M|+,A, 



^,-^|M|-,A, 



-^,+^|M|-,A, 



1 1 

-2'+2 



1 1 



|M|-,Ad 



-tt|M|-,A, 



2' 2 



+ (Q^^</)4+Q^^</>3) 



*d'(pi,+;p2,+ 



*d'(pi'+;?'2'-) 



(2^)2 



*d'(Pl'+:P2' + 

y *d'(pi'+;?'2'- 
^d'bl'-:?'^,- 

*d'(pi'-;?'2'+) 



+ (Q^^</)4+Q^^</>3) 



*d'(pi,-;p2,- 



*d'(pi,-;p2,- 



27r)2 



(2^)2 



(2^)2 



27r)2 



(2^)2 



(2^)2 



(2^)2 



27r)2 



(27r)2 



(2^)2 



(2^)2 



27r)2 



(2^)2 



(G.l) 
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in which the charge operators Q^" are defined such that 



(G.2) 



While for 7 + d ^ A + A, considering only helicity conserving amplitudes, we have 



that 



B 
B 
B 
B 
B 
B 
B 
B 
B 
B 
B 
B 
B 
B 



(27r)2 



Q 



JVi 



'8 ^^^^09) / ^^/{Pl,+;P2,~) 



fp^ 

(27r)2 



d^p 



^^/{Pl,-,P2,+) 



d'^p± 



Q 



(27r)2 

dW 
{27ry 

(fp 



1^1, 



(27r)2 

dW 
(27r)2 



(G.3) 



G.2 Helicity Amplitudes of Two-Nucleon Break-Up Reactions off ^He Target 



Replacing QIM amplitudes in Eq.(ni.l5) by nn helicity amplitudes of Eq.(ni.20) and 
using the antisymmetry of the ground state wave function with respect to the exchange 
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of quantum numbers of any two nucleons, one obtains the following expressions for the 
helicity amplitudes of two nucleon breakup reactions off the -^He nucleus, (Ai/,A2/,A, | 
M I A^, Am): 

For a positive helicity photon, 

^ + i +i AJM|+,A^ = b(^{Q''^+Q''-)^,J^^-{p„+;p,,+;p^,X^)^ 



+ (g^^<^3 + g^^<^4) / K^iPi,+;P2,-;Ps,Xs)j^ 



(27r)^ 

^,-^,AJM|+,A.4^ = B(^(O'^^+Q^^)</'2y*He(Pi'+;P2,+;P.,A.)0i 

+ iQ^^+Q"^)^,! KUPr,+;P2,~;Ps.Xs)^ 

i +i AJM|+,A^^ = b(^{Q''^ +Q''-)<l>sJ ^'^,^{p^,+;p2,+;p.,Xs)^ 

d^P± 



(G.4) 



and for a negative helicity photon. 



^ \a,|M|-,A^\ = -i?f(Q^i+Q^^)0i /vI/;.^(pi,-;p2,-;P.,A,)'^'^ 



i +1 AJM|-,A.4^ = -i?(^(Q^i+Q^^)03y'vi.^-(p,,-;p„-;p,,A,)|^ 



(O'^^03 + Q~^04) J ^^iip^,-;p2.+;Ps,X,) 



--,-r-,A,|M|-,A.4^ = -i?(^(Q^^+Q^^)02y *He(pi.-;p2'-;?'-A3)^ 

i,-^,A.|M|-,A^^ = i3(^(g^i+g^^)03y'*^-(p„-;p2,-;p.,A,)|^ 



(G.5) 
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where B = i^^^B^. Because the scattering process is considered in the center 
of mass reference frame in which the z direction is chosen opposite to the momentum of 
the incoming photon, the bound nucleon helicity states correspond to the nucleon spin 
projections ^ for positive and ^ for negative helicities. 
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